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PREFACE 


Over 15 years ago I began the preparation of AGARDograpb Number 186, published in 1974, the precursor with the 
same tide as this publication. The first publication provided a summary of the state of the art in the application of lasers in the 
field of flow diagnostics. Example applications were provided from a number of NATO laboratories, including many from the 
author's laboratory at the US Air Force Arnold Engineering Development Center. 

The present project began as simply an update and expansion of the earlier publication. I quickly came to the realization 
that the field has expanded and matured so much that a virtually complete rewrite was necessary. The resulting volume contains 
more new material than repeated and updated material. New hardware and technology have evolved, including detector arrays, 
small inexpensive personal computers, powerful software, a wide array of spectroscopic techniques, and many other new 
techniques. 

Literature surveys and visits to various laboratories produced several thousand excellent publications of techniques and 
applications. One of the difficulties in compiling this work was the need to limit its length and to choose only representative 
applications from the many fine examples I saw during the laboratory visits. I have added one entirely new chapter to the 
volume — laser spectroscopy — and retained the remaining four which cover (1) a brief introduction to coherent optics and 
lasers, (2) conventional flow visualization techniques that can benefit from the use of lasers as a source of light, (3) flow 
visualization and particle field holography, and (4) laser velocimetry. 

Major differences in the subject of laser flow diagnostics today are that it is a mature technology, being applied routinely in 
hundreds of laboratories, and many varieties of laser diagnostics systems are commercially available. This does not mean that 
the field has stagnated or even reached a plateau in producing new techniques. For example, laser fluorescence and molecular 
velocimetry made possible by excimer lasers is an exploding technology that will have critical application in high speed flows 
owing to revived interest in hypersonic aircraft. 

I fully anticipate the necessity to rewrite this again before the year 2000. 


II y a quinze ans, j'ai entame les premiers travaux en vue de la redaction de l'AGARDographie No. 186. Celle-ci, qui 
annonfait la presente publication du meme nom, est panic en 1974 et representait, a cette epoque, la synthese des 
connaissances dans le domaine des applications laser au diagnostic de I'ecoulement. II comprenait des exemples d'application 
foumis par un certain nombre de laboratoires de I’OTAN, parmi lesquels le laboratoire dont dependait 1’auteur, a savoir 1’US 
Air Force Arnold Engineering Development Center. 

Le present document fut a I’origine une simple mise a jour destinee a apporter un complement d'informations a 
l’AGARDographie 186. Or, j’ai condu tres rapidement que les connaissances dans ce domaine avaient considerablement 
evolues, atteignant un degre de sophistication tel que je serais amene a entreprendre une refonte quasi-totale du document 
original. L’AGARDographie qui en resulte comprend une part plus importante d’elements nouveaux que d'elements repris a 
jour. Du nouveau materiel a fait son apparition, grace aux technlogies nouvelles; notamment: des dispositifs de detection, des 
petits ordinateurs personnels de cout modique, des logiciels performants, ainsi que toute une gamme de techniques 
spectroscopiques et d’autres techniques nouvelles. Les visites que nous avons pu effectuer dans differents laboratoires, ainsi 
que la synthese de la litterature que nous avons faite nous ont permis de disposer de plusieurs milliers de publications d’un tres 
haut niveau concemant les technologies en question. 

Cependant, l’une des difficultes que j’ai rencontree lors de l'elaboration du present document a ete de limiter son volume 
en ne choisissant que des applications representatives parmi les nombreux exemples presentes, au demeurant tous d'une 
qualite remarquable. Tai ajoute un nouveau chapitre au present volume, concemant la spectroscopie laser, qui vient completer 
les quatre chapitres du document original: (1) introduction 4 I'optique coherente et aux lasers; (2) techniques classiques de 
visualisation de I’ecoulement susceptibles de beneficier de I’emploi de lasers comme source iumineuse; (3) visualisation de 
I’ecoulement et holographie du champ des particules; (4) velocimetrie laser. 

AujourdTiui, le diagnostic de I’ecoulement au laser est une technologie dont 1'efficacite n'est plus 4 demontrer a laquelle 
des centaines de laboratoires font appel tous les jours. D’autre part, est apparu sur le marche un grand choix de systemes de 
diagnostic au laser. 

Cela ne veut pas dire que les travaux de recherche dans ce domaine stagnent pour autant, le seuil d'innovation etant loin 
d'etre atteinL A litre d'exemple, nous citerons la fluorescence laser et la velocimetrie moleculaire, tributaires toutes les deux de 
la technique des lasers 4 excimeres; il s'agit 14 d'une technologie en plein essor, qui est appelee 4 jouer un role determinant dans 
le domaine de I'ecoulement 4 grande vitesse, suite au regain d’interet pour les avions hypersoniques. 

En conclusion, il me semble pertinent de signaler qu'une nouvelle refonte document s'imposera avant I'an 2000. 
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LASER APPLICATIONS IN FLOW DIAGNOSTICS 

by 

James D. Trollnger 
Chief SclenClet 

SPECTRON DEVELOPMENT LABORATORIES, INC. 
Costa Hess, California, U.s.A. 


SUMMARY 

This la an expansion and update of AGARDograph No. 186, Laser Instrumentation for 
Ploy Field Diagnostics which was published In 1974. A brief Introduction to coherent 
optIce and lasers Is given to provide a review and common language to the reader. The 
material la presented at a level to be understood by engineers who are working in the 
aerospace field and who are Interested in the experlaental application of lasers in 
research and developaent. An eaphasls la placed on applications In aerodynamics, and 
exaaples have been selected froa a number of laboratories In the NATO countries visited 
by the author while preparing the aaterlal. The applications Include holography, laser 
veloclaetry, spectroscopy, and the use of lasers simply as efficient light sources in 
conventional measurement methods. 


CHAPTER 1 

INTRODUCTION TO LASERS AND LIGHT 


1.0 BACKGROUND 

An elegant and accurate theoretical description of coherent light existed many 
years before such light was available to experimentalists, 2 ,3 especially those who were 
Involved in practical applications. The advent of the laser In the early sixties 
transformed an essentially academic theory into a powerful and practical tool, which 
could be easily applied to direct the path to a vast number of applications. It was, 
perhaps, this coablnetlon of the availability of a powerful new experimental tool, and a 
highly refined theory with which to describe It, that led to a virtual explosion of the 
field of epplled optics. The growth was accelerated not only In experimental fields but 
also In theoretical optica which began to feed on the ever Increasing number of 
experimental anomalies. The large nuober of practical applications, moreover, provided 
both interest and support for a significant increase In the number of people working and 
contributing to the field. In this chapter we will provide a brief summary of coherent 
optica and lasers for the purpose of providing the reader with the necessary language and 
basics to understand later chapters. 

1•1 Description of Light 

1.1.1 Haves and Wavefronts 

Light is a transport phenomenon for transporting electric and magnetic energy from 
one point to another. The physics of electric and magnetic fields is summarized by 
Maxwell's equations, which predict the existence of the energy transport by means of 
waves. To every point In space one con assign a value to electric and/or magnetic field 
strength. The waves are simply spatial and temporal variations In the field that 
propagate at Che speed of light. 

The general solution of the wave equation provides the values of the electric and 
magnetic field vectors In space and time. The two vectors constitute a tensor for each 
point In space. Fortunately, for the discussion here, we can, by certain specialization 
and approximation, actually describe light with simple scalar quantities. The most 
elementary solutions to the wave equation are sinusoids, and these describe laser light 
with sufficient accuracy to be tractable to many practical cases. When a coherent light 
wave passes through a region in space, the electric and magnetic fields vary sinusoidally 
with time at any point In the region, and they vary sinusoidally along any straight line 
through the region. 

Wavefronts are defined as surfaces of equal phase in the region where the 
electromagnetic field exists. The surfaces are normal to the direction of travel for the 
waves of Interest here (transverse waves). The wavefronts are separated by one 
wavelength and they travel at the velocity of light. When the wavefronts are planes, the 
light is said to be collimated and when they are spherical, the light is said to be 
converging or diverging. If the wsvefrontB are other than plane or spherical, they must 
be described mathematically with a aeries of sinusoids. 

1.1.2 Properties 

Light can be described bv vectors and scalars that can be manipulated and used in 
diagnostics. 

1.1.2.1 Polarization - For the cases to be considered here, the electric and magnetic 
vectors areperpendicular to each other, they are perpendicular to the direction of 
travel of the wave, and their magnitudes are proportional to each other (transverse 
waves). This allows description of the wave in terms of either vector. Even with these 
constraints, the vectors can still rotate about the direction of travel. The state of 
polarization of a wave describes the direction of the electric field vector. Linearly, 
polarized light has a single, fixed electric field direction. Elllptlcally, polarized 
light has a rotating electric field vector. Unpolarlzed light is a blend of many 
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possible polarization states. This property of light Is Important because It affects the 
way two light waves interfere with each other as well as the way light reflects from and 
transmits through materials. 

1.1.2.2 Amplitude and Intensity - At any point in space, the amplitudes of the electric 
and magnetic fields In an electromagnetic wave vary at an extremely high frequency 
(typically 10 1 * to 10*5 Hertz). A linearly polarized, plane wave that has a single 
wavelength, and which travels in a direction k parallel to the x-z plane in a uniform 
medium, can be described mathematically as 

u(x,z,t) - 

( 1 . 1 ) 

u Q (r)exp[l(nk.r - «t + e Q )] 

where a c is the phase of the wave at the origin at t - o; k is 2i/x; x is the vacuum 
wavelength; u 0 is the wave amplitude maximum; and n, the refractive index of the medium, 
is the ratio of the velocity of light in a vacuum to that in the medium (see Figure 1.1). 
Equation 1.1 is an elementary solution (or eigenfunction) of Maxwell's equations, and 
describes a wave which is said to be coherent with itself and with all other waves that 
have the same wavelength. 


x 



Figure 1.1 Propagation of a plane electromagnetic wave. 

The wave propagates in a direction normal to 
the wavefronts in the direction of It. 

In principle, one can measure the electric or magnetic field in the wave. However, 
the variation is much too fast for ordinary instrumentation. The energy flux In such a 
wave is proportional to the square of the electric field, and, therefore, it varies at 
the same frequency. A quantity which is measurable by normal sensors or photographic 
material is the time averaged energy flux, which will hereafter be called the "Intensity" 
of the wave. When averaged over times which are long compared to the period of the wave, 
it is given by 


I(x,z,t) - uu* (1.2) 

This equation will be used many times in the following chapters to determine the 
observable in an Instrument. In each case, the complex amplitude at a point is 
determined by vector addition of the asiplitudes of various waves contributing to the 
total field at the point. Then the observable is computed by using Equation 1.2. 

1.1.2.3 Coherence - When two electromagnetic waves are mixed, their electric and 
magnetic vectors add. If the waves have the same wavelength, then the difference in 
phase of the two waves will be constant at any location in space. When this is true, the 
intensity of the wave is maxinun when the two waves are in phase, and minimum when they 
are ISO degrees out of phase. Also the Intensity is maximized when the polarization 
direction of the two waves are the same. If the two polarization vectors are normal to 
each other, then Interference does not occur and the Intensity la simply the sum of the 
intensities of the two individual waves. 


EXAMPLE: The maximum net intensity, in a region where two equal amplitude, coherent, 
polarization aligned waves are mixed, is four times the intensity in one of the waves (in 
phase case), while the minimum intensity is zero (180 degrees out of phase). 


No real light wave comprises a single, exactly decerminable wavelength. The 
frequency content of all real light waves is broadened, by the many broadening mechanisms 
existing in the wave producing phenomena, and because all waves have finite extent in 
space. The practical effects of broadening can be understood by examining the appearance 
of composite waves that comprise two or three frequencies. Figure 1.2 illustrates the 
appearance of waves comprised of one, two, three, and many frequencies. The single 
frequency wavefront as already described is simply a sinusoid characterized by a wave¬ 
length and velocity. When two waves are mixed, the resulting amplitude varies along a 
straight line as the two waves move In and out of phase with each other. (This is 
derived In Chapter 2.) If the waves are equal in amplitude, then the net amplitude 
varies from a maximum to zero along the straight line, or at a point as a function of 
time. By knowing the wavelength difference, ana by determining how many cycles must pass 
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before Che two mvci cove back to the ease relative phase, it la easy to compute a beat 
frequency 
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(d) Many frequencies 

Figure 1.2 Appearance of an electromagnetic wave comprising 
several frequency components. 

If such a wave is split into two equal intensity parts, and If the two parts are 
remixed after having traveled different distances, then the resulting intensity depends 
upon the relative phase of the beat waves. If the beat waves are 180 degrees out of 
phase, then the resulting Interference is effectively reduced since the amplitudes of the 
two waves differ at all times, making it impossible for either of the two waves to cancel 
the other. As more frequencies are introduced to construct the wave, the distance over 
which the two waves can efficiently Interfere becomes smaller and smaller. The distance 
of efficient Interference is called the coherence length of the wave. If a wave la split 
into two waves, and if the two waves are transmitted over different distances before 
being remixed, the Interference (fringes) will be observed to deteriorate as the 
difference In distance increases. That distance difference, at which interference ceases 
to be observable, is called the temporal coherence length of the light wave. 

A similar effect occurs when one attempts to interfere waves extracted from 
different places on the wavefront; that is, from two places taken from a plane which is 
normal to the direction of travel. As the distance is increased, the Interference 
deteriorates until it is no longer observable. This distance is called the spatial 
coherence interval. 

1 .2 Interference 

1.2.1 Plane Waves 

We now examine the interference of two coherent plane waves. Cast Equation 1.1 
into Cartesian coordinates and examine the effects in the xy plane, which is normal to z, 
the direction of travel. In the following, we carry only the c and z coordinates since 
dependence on x and y is the same. 

u(x,z,t) - u c (x,z,t) expi[kn(zcosa + xslna) • «t + * c ] (1.4) 

The intensity is proportional to 

I(x,z,t) - uu* * [u Q (x,z,t)]^ 0.5) 

We will consider cases in which I Is a slowly varying function of time, where I Is 
time averaged over times large compared to the period of the wave but small compared to 
measurable time. Consider the mixing of two waves of like polarization traveling along 
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two different directions. The rigorous solution for the resulting intensity is not 
simple, but an approximate solution for small angles between the two directions of travel 
can be found simply by adding the complex amplitudes and squaring. 

I(x,z,t) - (u, + U 2 )(uf + uj) - |u 01 | 2 + ! « 02 1 2 + u 1 u 2 + u * u 2 (K6) 

Combining Equation 1.4 and 1.6 leads to 

Kx.x.t) - u Q f + u Q | 

+ 2u 01 UqjCoi [zn(k,cos 0) - k 2 C08a 2 ) (1.7) 

♦ xn(k,sinaj - k 2 slna 2 ) - tmt + 4*] 

Equation 1.7 describes s set of alternate bright and dark fringes that are not stationary 
In tine. We shall now consider soae special cases and useful examples of Equation 1.7. 

Case 1 - Sane Wavelengths and Phase at the Origin 

k, - k 2 - k. An - 4* - 0, 

2 2 

I (x, z) - u Q1 + Uq 2 + 2u 01 u 02 cosk[zn(cosa 1 - coso 2 ) 0 . 8 ) 

+ xn(sin»| - sina 2 )] 

The resulting set of fringes is illustrated in Figure 1.3. They comprise sheets of 
light and dark in space. The equation for these planes of equal Intensity is found by 
setting the argument in Equation 1.6 to a constant. Thus, for n - 1, 

(coa<*. - cosao) 

X " (sin., - »rn.“) 8 * conilt - (K9 > 


y 



Figure 1.3 Interference of two plane waves. Equally spaced fringe planes 

bisect the angle between the two beams. The fringes are stationary 
for two beams of the same wavelength. But at different wavelengths 
fringes move normally to the planes. 

For two particular planes the fringes are especially simple. 

(a) One Beau Normal to the Plane of Observation (a, - 0) 

I(x,z) - u 0 2 + Uq 2 +2u 01 u 02 cosk[x(1 - cosa 2 ) - xsina 2 ] (1.10) 
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EXAMPLE: In any xz plane of observation the fringe frequency is 

f x - slno^/i (1.11) 

At A * 0.5 microns, the fringe frequency for a ■ 10 arc seconds is approximately 

1 fringe/cm, while at o • 10 degrees, f is 3,600 fringes/cm. 


The fringe planes are given by 

x • ztan ( 02 / 2 ) + const. 


( 1 . 12 ) 


Therefore, they are parallel to the bisector of the angle between the beams. (This is 
true regardless of the observation plane.) Note that when a ■ 180° (the beams are mixed 
from opposite directions), the fringe planes are Derpendicular to the directions of 
travel of the beams. The fringes are separated by \/2. 


(b) The Plane of Observation is Normal to the Bisector of the Angle Between the Two Beams 
(<*1 “ - 02 ) 


Kx) - uj + u 0 2 + 2u 0 1 u 02 
cos [2nk(xsina)] 


0.13) 


The fringe planes are given by x - constant and are, therefore, parallel to the z-y 
plane. The fringe frequency, just twice that given in the previous case, is (2nsina)/x. 


Case 2 - Same Phase at the Origin, Approximately Same Wavelengths, ay * -02 

I(x,z.t) - u 0 ? + Uq 2 + 2u 01 u 02 
cos [2nxksino - &ut] 


(1.14) 


Two important features are to be noted. First, the lack of z dependence places the 
fringe planes parallel to zy plane as in Case 1b. Secondly, the argument of the cosine 
term in Equation 1.14 will be constant when 

which can be interpreted as a motion of the fringe planes perpendicular to themselves (x 
direction) with velocity 


v - Ad)/2knsina (1.16) 

Therefore, if a tiny, stationary scatter center is being illuminated by two such waves, 
the scattered light intensity will oscillate with frequency f, where f is the frequency 
difference of the two waves. If the particle moves in the direction of fringe motion, f 
will decrease. Motion in the opposite direction increases f. This principle is the 
basis for determining the particle velocity and direction. 

Case 3 - The Two Waves Travel Along the Same Line 


I(x,z,t) - u 0 ? + u 0 l * 2o 0 1 u 02 
cos[zn(k 1 - k 2 ) - Awt + a*] 


(1.17) 


At any position, z, the intensity varies at a frequency equal to the difference 
frequency of the two waves. This describes the wave illustrated in Figure 1.2. 


If the two waves have the same wavelength, then Equation 1.17 reduces to 

I(x,z,t) - u 0 | + u 0 2 + 2u oi u 02 co8 l4f] (1.18) 

and the intensity depends upon the phase difference of the two waves. Most 
interferometers can be explained in terms of the above equations. 


1.2.2 Spherical Waves 


The time independent amplitude in a spherical wave source located at r « 0 can be 
expressed as 

u(r) * u Q exp(ikr)/r (1.19) 

The quadratic approximation of this expression is 

u(x,z) - u Q exp[lij(x - *p) 2 ] (1.20) 

where x is the location of the source and x is the observation point. If two such waves 
are mixed, the resulting Interference fringes are curved. The amplitude is 
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( 1 . 21 ) 


u(x.z) - UQjexpIl-j^tx-Xp, ) 2 ]+ Uo^xpli^tx-Xpj) 2 ] 

and the intensity is 

I(x,z) - 

Note that when a 
plane of observation is 


1.2.3 The Moire Analogy 

A useful and instructive moire analogy can be constructed with grid patterns. 
These have proven so useful in instruction that a set of grids has been included in the 
appendix to allow the reader to construct his own set of transparencies. Representing a 
plane wave with a set of equally spaced lines in which the direction of travel is normal 
to the lines (as in Figure A1 ), one can then produce a transparency wi^h this as a 
master. When the transparency is superimposed on the master, a set of moire fringes can 
be observed as shown in Figure 1.4. Note that the fringes bisect the angle between the 
directions of travel of the two waves as predicted. As the angle between the two 
directions is increased, the fringe density increases. If the relative phase of the two 
waves is shifted by shifting position of the master, the fringes shift as predicted. 
Many of the intricacies of Interference can be illustrated by using the grids and a small 
amount of ingenuity. 


- 0 ? 


u 02 + 2 u 01 u 02 cos fx K*-* Dl ) 


„ u-‘ p W ' (*-*p2)' 1 (K22) 

spherical wave is mixed with a plane wave, the intensity in the 


I(x) 


2 2 
u 01 + u 02 


2u 


01 “ 02 ^ 




(1.23) 



Figure 1.4 Mixing of two waves at an angle. This moire pattern 

represents fringes seen from the top view of Figure 1.3. 

Figure A2 is included to represent a wave of slightly different frequency. By 
superimposing a transparency of A2 onto A1 , the beat waves described above are seen as 
shown in Figure 1.5. 



Figure 1.5 The beating of two waves of slightly different frequency. 

One of the wavefronts continues to advance ahead of the 
other moving in and jut of phase each time it advances by 
one-half wave. 
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Figure* A3 and A4 are Included to represent spherical waves. Now superimpose a 
transparency of A4 onto AA. The resulting fringes are illustrated in Figure 1.6. Note 
how the fringes move, as a relative motion between two light sources is simulated by 
shifting the transparency around. The mixing of two spherical waves of different 
frequency causing a beat wave Is illustrated by superimposing Figures A3 and AA as shown 
In Figure 1.7. We will return to the analogy later to explain holography and laser 
Doppler anemometry. 



Figure 1.6 The interference of spherical waves of the same wavelength. 



(a) Two waves of different wavelength (b) Two waves of different frequency 

emanating from separate points. emanating from the same point. 

Figure 1.7 Interference of two spherical waves of different wavelength. 

1•3 Interaction of Light with Matter and Boundaries 

1.3.1 Diffraction 

When light Interacts with matter or with a boundary, the resulting field becomes 
much more complicated. This constitutes the placing of boundary conditions on the 
Maxwell equations. The solutions then become series expressions or Integrals Instead of 
simple sinusoids. It is useful to think of diffraction as an interference phenomenon. 
When light passes over a boundary such as an edge, the light which strikes the edge is 





scattered In every direction and thla light mixes with the unacattered light. So any 
given point in a plane of observation receives light of two origins, scattered and not 
scattered. Since the distance to the point Is different for the two sources, the two 
coaponents of light will have different phases and can Interfere constructively or 
destructively, depending on the phase difference. Therefore, one expects to observe 
fringes In the vicinity of an edge (Figure 1.8a). There la no such thing as a "shadow" 
in the usual sense for coherent light, only a diffraction pattern. 


Consider the appearence of light after having passed through an aperture of 
diameter, d (Figure 1.R>,c). Near the aperture, the light distribution Is fairly uniform 
across the aperture. As one moves further along the optical axis, a series of 
diffraction rings can be observed. This Is called Che Fresnel region, or the near field, 
and the pattern will continue to exhlblc this appearance aa long as the distance Is less 
than about d 2 /i, a distance which la appropriately termed one far field distance or the 
distance to the far field. When the observation distance reaches and passes the far 
field distance, the appearance of Che diffraction pattern changes drastically. The light 
which was scattered from the edge remains In phase with the unacattered light over a much 
larger region of the observation plane. The lntenelty Is given approximately by 


I(r) 


j2 2J.(ird/xz) 2 

[ - -ba /i i — l 


(1 . 24 ) 


A cencral part of the diffraction pattern (sometimes loosely referred to as Airy's 
pattern), containing about 84X of the energy, dominates and, for a fairly large distance, 
Is actually smaller than the aperture Itself. Its radius Increases with z and Is given 
by 


r a » 1.22 xx/d (1.25) 

so that the half divergence angle, 9 d , of this portion of the beam Is given by 

tan» d - e d - 1.22 x/d (1.26) 

A similar result can be observed when light passes over an edge or an obstructing 
object. Light scattered from an edge has a phase difference, as, with unscattered 
light. For the small angle case, this difference along the z-axls is given by 

4 *(*)‘rf- °- 27) 




Figure 1.3a Interaction of light with edges. A diffraction pattern can be observed on 
a screen located at z. Light scattered from the edge (top figure) adds In 
phase with unscattered light when It travels a distance nx further 
producing a bright fringe. Light scattered along the centerline axis 
(bottom figure) produces a bright center when It Is In phase, having 
traveled nx further than unacattered light. After the distance z, where 
the path difference of scattered and unscattered light Is much less than x, 
the center stays bright. This is called the far field. 
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Figure 1.8b Optical arrangement (x - 633 nm) 



a. z* 0.2 cm b. z* 13.5 cm c. z-40cm 



d. z * 80 cm e. z • 120 cm f. z -160 cm 



g. z- 200 cm h. z- 240 cm i. z- 284 cm 


Figure 1.8c One-centimeter aperture (d “ 1 cm) 

where x Is the distance along the normal to the edge in the viewing plane. The intensity 
maxima occur when 

A*(z) m 2* N (N an integer) (1.28) 

When coherent light is used to produce images, the behavior of the diffraction 
pattern can be used to properly focus on the image. Unlike an incoherent image which 
exhibits a blurred edge when not in focus, the coherent image makes a transition in which 
the diffraction pattern can be observed to collapse into the focused edges as a focus 
position is approached. When the coherent image is badly o* t of focus, little change at 
all can be observed in the diffraction pattern as one moves toward the focus position. 
This point will be discussed further in the chapter on holography. 

1.3.2 Reflection and Refraction 

When light strikes a surface, some of the energy may be reflected from the surface, 
some transmitted into it, and some may be absorbed by it. If the surface is a perfect 
dielectric material, no absorption takes place as long as the intensity of the field is 
not strong enough to tear away the bound electrons. Flgpre 1.9 summarizes the 
interaction of the light with a dielectric surface. The angle of incidence and the angle 
of reflection are equal. The angle of refraction, Or is given by Snell's law: 
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Figure 1.9 Reflection and refraction of light at a dielectric aurface. Reflection and 
transmission depends upon polarization direction 

The incident ray, cranaaltted ray, and reflected ray stay in a single plane, the "plane 
of Incidence." If the light is divided Into Its polarization directions, one component 
In the plane of incidence and one normal to the plane of Incidence, then the two 
components can be treated Independently. For example. If all of the light Is polarized 
with a polarization vector In the plane of Incidence, the polarization of the transmitted 
and reflected light remains In this plane. 


A few simple special cases occur so often that It is useful to mention them here. 
There la an angle of Incidence, called Brewster's angle, at which light polarized in the 
plane of Incidence la totally transmitted (Figure 1.10). This angle Is given by 

a. ■ tan^n'/n (1.30) 



Figure 1.10 Refraction at the Brewster angle. Light polarized In the plane of 
incidence Is transmitted completely. 

When light travels toward a material of lower refractive index, there Is an angle 
called che angle of total Internal reflection, at which all light Is reflected 
(Flgnre 1.11). This angle Is given by 

alne c “ n'/n (1.31) 

These two conditions are used extensively In experimental optics. 

1.3.3 Propagation 

Because of diffraction, a light wave changes in Intensity as It propagates through 
space. The case of che uniform circular light wave han already been created above, and 
the Intensity distribution at some propagation distance, z, was given by Equation 1.24. 
Another case of practical Importance la the case of a light wave which has a Gaussian 
Intensity distribution. This particular distribution has Che unique and useful property 
that It maintains its form irtille propagating through space. The Intensity has the form 
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Figure 1.11 Total Internal reflection. Beyond a critical angle light traveling from a 
higher to a lower refractive Index la totally reflected. 

I(r) - I 0 exp(-2r 2 /a 2 ) (1.32) 

where Iq Is the central Intensity and a is the radius to the 1/e2 intensity. The radius 
increases according to the relation 

a( 2 ) - a[! + (- 41 7 ) 2 ] 1/2 (1.33) 


The propagation direction of a light wave will change If the refractive Index of 
the medium In which the wave la traveling changes along the wavefront. This can be 
explained by recalling that the direction of travel la normal to the wavefront and that 
the wavelength becomes shorter In a medium of higher refractive Index (Figure 1.12). The 
light wave la deviated by an angle 



Figure 1.12 Deviation of a light ray by a refractive Index gradient. The ray bends In 
che direction of the refractive Index gradient. 

where x Is normal to the original direction of travel. Assuming very small deviation of 

the ray away from z. Its phase is 

y(x.z) - 2. ^ ^ £ n(x, z')dz 1 (1.35) 

The second integral is called the optical pathlength fron 0 to z. It is the number of 

wavelengths existing between 0 and z times the vacuum wavelength. For example, since the 

wavelength is shortened in a medium with refractive index greater than unity, the optical 
pathlength becomes larger. In vacuum, optical pathlength is, of course, just 
pathlength. In many materials, n can be related to density, p, by the Dale-Gladstone 
constant, G. 

n - 1 + G p (1.36) 

so that 
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and 


dg G do 

<Tz 1 + G <Ix 

♦ (X,*) - ^ [1 + Gp(x, x 1 ) Jdx 1 


0.37) 


(1.38) 


These changes alter the aapllcude and the phase of the light over an observation 
plane. The use of this effect to observe flow fields will be described in Chapters 2 
and 3. 


1.3.4 Speckle 


When a light diffusing surface is viewed or photographed in incoherent light, it 
exhibits a more or less smooth appearance. When the same surface is viewed or 
photographed in coherent light, it exhibits a grainy appearance. This graininess is a 
statistical phenomenon that is called speckle. There are two categories of speckle, one 
caused by the limited resolution of the viewing system and the other caused by the 
limited extent of the illuminating system (Figure 1.13). When a finite sized diffuse 
surface is used to illuminate another surface that is located at a distance, z, the 
intensity on the second surface will exhibit speckle having a mean diameter of 


d sp " 1-22 (xz/D) (1.39) 

The phenomenon can be explained in the following way: any point on the observation 
plane is Illuminated by a confluence of light that originated at the many different 
points on the diffusor. Each of these points is characterized by a phase that is 
different from other points in a random way. There is a probability that the sum of the 
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a) Objective Speckle 



Figure 1.13 Speckle phenomena. 


sources will add up to produce a bright point in the observation plane; and, likewise, 
there is a probability that a dark point will be produced. 

When a diffuse surface is imaged, every resolution element in the image receives 
light from a large number of randomly phased sources that make up the resolution points 
(such as the high and low points that constitute the surface texture). There is a 
probability that the different components of light making up the resolution points will 
either add up sufficiently in phase to produce a bright point or out of phase to produce 
a dark point. The average speckle size, therefore, can be expected to be the size of the 
resolution of the imaging system. 


d_ - 1.22 xs/d (1.40) 
sp l 

The resulting useful resolution of the imaging system is, therefore, actually three 
or four times less than that given by the normal equation for resolution since the image 
will consist of speckles of that size. The highest useful resolution is attainable when 
direct, not diffuse, illumination can be used so that speckle does not occur. 
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1.4 Control of Lixht 

In this section, we shall summarize and discuss some of the most basic elements for 
the control and use of laser light. The special characteristics of laser light provide a 
spatial and temporal controllability that cannot be attained with Incoherent light and 
have fostered many new applications. In general, coherent light can be delivered into 
smaller spatial regions, at higher intensities. In smaller time elements, over greater 
distances, with more precise modulation, and, because of single wavelength, can be 
separated more easily from background radiation. The following discussion covers some of 
the devices and physics that make these Important characteristics possible. 

1.4.1 Lenses and Mirrors 

Lenses and mirrors transform an input light wave from one form, such as planar, to 
another such as converging. The input space Is called the object space, while the output 
space is the image space. The object space can be emulated exactly in image space only 
for very special cases; however, the approximations are useful especially for engineering 
design. 

The most basic function of a lens is the conversion of a plane wave to a converging 
wave. When a plane wave strikes a positive lens, it converges to a focal point that is 
located at a distance from the lens, f, called the focal length of the lens. When the 
wave strikes a negative lens, it appears to diverge from a point on the input side of the 
lens, also called the focal point. 

Points in the object space are approximately related to points in the image space 
by the thin lens equation. 


s' - sf/(s - f) (1.41) 

Where s is the distance from the object point to the lens, s' is the distance from the 
lens to the image point, and f is the focal length of the lens. 

The image is never exactly the same as the object. Three things limit the accuracy 
with which a wavefront can be changed from one form to a second desired form: 1) the 
design of the lens (such as the specification of material, surface curvatures, and the 
combinations of these), 2) the accuracy of manufacture of the lens, and 3) diffraction. 
The effects introduced by the first two are called geometrical aberrations, and the 
limits set by the third are called diffraction limits. 

A high quality, well designed lens, in which geometrical aberrations are small 
compared to the limiting ef feces of diffraction, is loosely termed a diffraction limited 
lens. A lens cannot be produced to be diffraction limited for all transformations. The 
most common diffraction limited lens is a collimator; it transforms a plane wave into a 
converging wave or vice versa. A parabolic mirror is a diffraction limited collimator. 
Lenses and mirrors must be illuminated from the correct direction to produce diffraction 
limited images. Additionally, lenses must be illuminated from the correct side. 

A converging wave passes through a transition region of minimum beam diameter 
wherein the wave is plane before beginning to diverge. The intensity distribution at the 
f ocus depends upon the intensity distribution across the lens. Figure 1.14 a-d 
summarizes the results of illuminating a lens with two different plane wave intensity 
distributions, uniform and Gaussian. Note that the effective focus diameter and position 
depends upon the diameter of the input beam. The lens performs a focusing action on the 
beam while diffraction causes a diverging action. The net result is that the minimum 
occurs at a different point than would be predicted from geometrical optics alone. Since 
many lasers produce a wavefront of Gaussian intensity profile, it is extremely important 
to understand Gaussian optics and the difference between geometrical optics predictions 
and the more accurate diffraction optica results. Geometrical optics assumes that light 
can be focused to an infinitely small point. This, of course, does not happen. The 
image distance in diffraction theory actually depends on the input beam diameter and 
intensity distribution. The following are some significant characteristics of beams with 
Gaussian intensity profile. 

(1) Propagation through empty space - the intensity profile remains Gaussian, a unique 
property of Gaussian intensity profiles. 

(2) Divergence from a source located at the front focal plane - the effective image 
position (minimum diameter) is at the back focal plane. 

(3) Truncation - the beam will not remain Gaussian if it is truncated. The Gaussian 
should not be truncated below two 1/e 2 diameters to retain Gaussian properties. 

(4) Energy removal by truncation - a curious result is that the percentage of total 
energy removed by truncation is the same as the percentage of maximum intensity at 
which truncation is done. 










(c) Diffraction optica vs. geometrical optics 
Figure 1.14 Focusing laser beans. 









Example: A Gaussian beam passes through a circular stop. The Intensity of the beam at 
the edge of the stop is one percent that at the center. Therefore, one percent of the 
energy will be removed from the beam. 


At the focal point, 84 percent of the energy passes through the circle with 
diameter d A for uniform illumination, and for a Gaussian profile 87 percent passes 
through the circle with diameter d(e“2). The distance in the focal region over which the 
waves are nearly plane is referred to as the depth of field af e . For the Gaussian case, 

4f e " f e f * /,aZ 0.42) 
1.4.2 Scanners, Deflectors, Modulators, and Shutters 

Laser beams can be scanned, deflected, and modulated at high speeds with mechanical 
and electro-optic devices. Mechanical devices include rotating or vibrating prisms, 
gratings, and mirrors that can scan a beam over wide angles at kilohertz rates. 
Vibrating mirrors are attached to electromagnetically driven components that oscillate at 
rates up to tens of kilohertz. The highest rates are attained when torsion elements are 
driven at resonance, but it is possible to drive these at programmable angles at 
kilohertz rates. Rotating elements can attain even higher rates. Holographic gratings 
or other holographic optical elements also have provided an extremely useful range of 
capabilities. 

Electro-optical scanners can deflect beams at megahertz rates. The most common of 
these is the Bragg cell, an acousto-optical device. In this device, the beam is 
diffracted from ultrasonic waves that are generated it a material by plezo-electrically 
driven crystals. 

Other devices modulate the polarization properties of light at megahertz rates. 
When used in conjunction with polarization filters or polarization dependent reflectors, 
the beam can be deflected or shuttered at the same rates as the polarization modulation. 

Other types of modulators modulate phase only. A Pockels cell can be used in this 
mode. Any element that changes optical pathlength can be used as a phase modulator. A 
somewhat unique phase modulator is an optical fiber that is stretched electromechanically 
or piezo-electrically while light passes through it. By wrapping the fiber around the 
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piezo-electric element, a full wavelength shift in the optical pathlength of light tra¬ 
veling through the fiber la attainable with a snail amount of voltage applied to the 
crystal. 

Coherent light can be shuttered more easily than incoherent light, because nore 
parameters are controllable for the shuttering process. Shuttering elements employ both 
mechanical and electro-optical techniques covering a vast range of shuttering times down 
to nanoseconds or less. One simple form of shutter employs a scanner, as described 
above, to sweep a focused beam across a small aperture. Shuttering times less than one 
microsecond are possible with this method. Such times are attainable even with 
mechanical shutters if the beam is focused to small enough size. When the size of the 
beam cannot be made small, mechanical methods are limited to the millisecond regime. 

The fastest shutters are electro-optical devices. A Pockels cell, when used as 
described above, can shutter a laser beam for times as short as nanoseconds. 

Pulsed lasers themselves can provide the shuttering action if all other light is 
blocked from the sensor. This can be accomplished by using narrow band filters that pass 
only the laser wavelength and a narrow band around it. The effect can be further 
enhanced if the room light is limited to a band that is widely different from the 
laser. Common fluorescent lights emit little or no radiation at wavelengths above 
650 nm, with the result that a high pass cutoff filter such as Scott RG650 can be used to 
cover the film or sensor in a fully lit room when a laser above this wavelength is used. 

1.4.3 Holographic Optical Elements 

If any two wavefronts of light are known, then a hologram can be produced that 
converts one of the wavefronts into the other. Holographic optical elements (HOE) are 
holograms that are designed to convert a wavefront of one specific type into a different 
desired wavefront. A HOE is made by using one of the two wavefronts as a reference wave 
and the other as the object wave. Then, when the hologram is illuminated with the 
reference wave, the desired wavefront emerges. 

HOEs can be used to perform the function of virtually all optical elements such as 
gratings, lenses, and mirrors, or even much more complex devices. Primary advantages 
include low cost, small thickness and weight, and the fact that they can be produced on 
any shape of surface. The major disadvantage is that they are usually not achromatic. 

1.4.4 Spatial Filtering 

Laser beams almost always contain diffraction noise, flashlamp or plasma 
illumination, and transverse mode structure. Much of this can be removed by a procedure 
of spatial filtering. The beam is focused by a high quality lens through a tiny pinhole 
(Figure 1.15). The intensity distribution of the focused laser beam depends on the 
wavelength, lens diameter, focal length, and intensity distribution across the lens. 



Figure 1.15 Spatial filtering. Rays that are not collimated 
in the input beam are removed by the pinhole. 

The pinhole diameter can be estimated by examining Figure 1.14. The aperture 
having the diameter of the Airey disk will pass 84 percent of the collimated light 
striking the lens, while removing moat of the uncollimated light. Actually, most laser 
beams of interest are Gaussian in intensity profile rather than uniform. For continuous 
wave (CW) lasers, pinhole diameters range from a few micrometers in diameter upward, 
depending upon f. For high power pulsed lasers, the lens focal length must be Increased 
to reduce the intensity at focus, requiring that the pinhole size be Increased. 
Otherwise, the intensity can become high enough to ionize the air or damage the 
pinhole. Air breakdown does not necessarily render the process useless. The plasma 
absorbs the light and shortens the pulse length and may eventually destroy the pinhole, 
but the process works at least until the pinhole is destroyed. The use of diamond 
pinholes lessens the problem. 
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Figure 1.16 illustrates a spatially filtered beam with various size pinholes. To 
emphasize the filtering action, a 10 alcroaeter dlaaeter wire was stretched across the 
back side of the 20 power alcroscope objective. Ulth a 30 micron dlaaeter pinhole, 
alaost all of the diffraction noise is removed. Ulth no spatial filtering, not only is 
diffraction noise present in the beam, but also multiple reflection Interference rings 
are seen. 



for Illustrative Purposes 



Figure 1.16 Spatial filtered beans with various pinhole sizes. Note that 

higher spatial frequency information passes Into the final Image 
as the pinhole gets larger. 

1.4.5 Illuminating an Object Field 

Three different classes of coherent illumination, direct, diffuse,and structured 
can be used. Diffuse Illumination suffers the problem of speckle, but provides the 
easiest viewing by eye. Direct illumination provides speckle free light, but is more 
difficult to view by eye. The most common form of structured light is the so-called 
sheet of light or light screen produced with a cylindrical lens. This provides a 
convenient way to view activity in a plane, without suffering the optical noise of light 
scattered from the remainder of the volume that is not under examination. Other forms of 
structured light Include multiple sheets of light or rays of light. 

Other choices of object Illumination include front lighting and backlighting. 
Backlighting provides for high resolution and high contrast viewing of the perimeter of 
an opaque object, while front lighting provides for the viewing of surface detail of the 
object. 
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1.5 Lasers 4 * * 

1.5.1 How They Work 

Lasers are characterized by two major features, the lasing medium and the physical 
arrangement that forma a cavity, or Interferometer, for selecting photons of a particular 
frequency and feeding them back Into the medium, where they are amplified before leaving 
the cavity. The technical knowledge required to build a laser became available with 
Einstein's theory of stimulated emission in the early 1900s. Nevertheless, the first 
operational laser was built by Maimon in 1960, nearly a half century later. Within 
months of the announcement of his achievement of lasing action in ruby, investigators 
from laboratories the world over reported lasing in a large number of other materials. 
Today, hundreds of commercially available lasers employ many different lasing materials 
to provide useful characteristics. In the following paragraphs, the individual 
phenomenological elements of a laser are discussed and then combined to produce a laser. 

1.5.1.1 Stimulated Emission - An excited molecule can decay to its ground state by any 
one of several different mechanisms. The decay must be accompanied by an energy transfer 
of the residual energy, aE, either to another molecule (radiationless transition) or to 
an emitted photon. The emission of a photon can occur in either a spontaneous process or 
a stimulated process. The probability for stimulated emission at radiation frequency, 
f " AE/h, where h is Planck's constant, is increased sharply when the excited molecule is 
illuminated by radiation very close to f. These photons are said to be resonant with the 
excited state, and their presence greatly shortens the average lifetime of the excited 
state. Stimulated radiation leaves the excited molecule in phase with the stimulus field 
(Figure 1.17). Transitions are possible only between selected energy levels. When a 
molecule is struck by a resonant photon, the probability for absorption of the photon by 
a molecule in its ground state is the same as the probability for emission of the same 
photon by an excited molecule. 
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Figure 1.17 Molecular excitation and decay. 

1.5.1.2 Population Inversion - The lasing medium can be excited in many ways, but the 
resulting excited material is the same. When a collection of resonant photons passes 
through an excited material, some of them will be absorbed and some of then will 
stimulate emission of more of the same photons, resulting in an amplification of the 
illuminating radiation. The condition required for a net amplification is simply that 
the number of stimulated emissions exceed the number of absorptions and other losses. 
Since the probability of stimulated emission and absorption is the same, this condition 
occurs when the number of molecules that populate the excited resonant higher state 
exceeds the number of molecules that populate the lower state. This condition is 
appropriately described as population inversion between the two states, since the 
condition Is the inverse of the usual stable condition. 

1.5.1.3 Amplification - When resonant photons pass through a material where population 
inversion Exists, the net energy in the traversing beam will Increase. The material acts 
as a light amplifier for photons that have energy equal to the difference of energy 
between the population inverted states. Amplifiers are used to increase the output power 
of a laser as well as to aid in the amplification and detection of weak coherent light 
waves. In some instances a laser nay incorporate three or four amplifiers in tandem. 

1.5.1.4 Feedback-The Oscillator - When the lasing material is placed between two 
reflectors^ the amplified right wave is fed back into the material where it is further 
amplified. If a light wave gains in strength after making a round trip between the two 
mirrors, the lasing threshold is said to have been exceeded and lasing action takes 
place. The two mirrors may be planes, curves, or combinations of plane ana curve. 

The space between the two mirrors is called the cavity. The energy in the cavity 
could become so high that It destroys the mirrors or the lasing material, but by making 
one or both of the mirrors partially reflecting, some of the light leaks out of the 
cavity (see Figure 1.18). Waves travel in both directions in the cavity, forming a 
standing wave whose wavelength satisfies L * qj/2 where q is an Integer. 
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(c) The laser cavity 
Figure 1.18 Lasing Action 
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1.5.2 Coherence Improvement and Control 

1.5.2.1 Frequency Broadening - Temporal coherence Is tantamount to single frequency; 
spatial coherence Is tantamount to simple, plane or spherical wavefront. Improving the 
coherence of a laser amounts to making the cavity strongly supportive of a single 
frequency, simple wavefront and strongly unsupportive of anything else. If the energy 
associated with the inverted state were the same for all molecules, then a single 
frequency would result. Since the energy is not the same for all molecules, the 
frequency of light emitted from an excited molecule upon decay, AE/h, is not the same for 
all molecules. The energy broadening mechanisms Include: 


(a) Natural Broadening due to the uncertainty principle, which states that the product 
of the state lifetime and the energy uncertainty must exceed Planck's constant, 
giving rise to a frequency broadening of typically about 15 MHz. 

(b) Pressure or collision broadening is the alteration of the energy of a molecular 
state, through collisions or close proximity with other molecules. This effect is 
most important in solid lasing materials or gas lasers having higher pressures ( in 
excess of one torr). 

(c) Radiation Broadening related to the finite width of the stimulated emission 
resonance curves. 

(d) Zeeman Broadening caused by magnetic fields which alter the energy states. 

(e) Doppler Broadening caused by the different kinetic energies of the molecules that 
are transferred to the emitted photon, typically of the order of several GHz and 
greater for higher temperature lasers. 

1.5.2.2 Transverse Electromagnetic Modes (T.E.M.) - The solution of the electromagnetic 
cavity equations is a series of eigenfunctions with indices p,l,q. Each of these 
functions represents a mode of operation of the laser and is called a transverse 
electromagnetic mode (T.E.M.). A laser may operate in many modes at once, since the 
general solution of the equation is a combination of the eigenfunctions; however, the 
lower order modes usually predominate. A laser can be forced to operate predominately in 
a single mode, usually p “ 0, 1 * 0, q, or the TEMqq mode. The various p and 1 modes are 
characterized by different intensity distributions about the optical axis of the laser. 
The amplitude of the p and 1 mode for a laser with cylindrical geometry is given by 


u (r.t) - U c (-^|) V (2^) exp (=4) 
r a a 

U Q ~ amplitude at the optical axis 
L* p ~ Laguerre Polynomial 
a - 1/e 2 radius 
♦ ~ angular coordinate 

Consider the p - 0, l « 0 case. Since L°(x) - 1 

U (r,.) - U(r > - U o “P t^) 

a 


. sin 
‘cos 


1 l * 


(1.43) 


(1.44) 


Therefore, the lowest order mode has a Gaussian amplitude and intensity profile. 

The frequencies of the various t, p, and q modes differ, the amount depending upon 
the laser geometry. If, for example, the laser mirrors are plane, the frequency does not 
depend upon t and p. In general, the various modes are not coherent with each other. 
The lowest order mode is most widely used because the Gaussian intensity distribution is 
less divergent, more easily focused, and is most easily used in calculations; the higher 
the mode structure, the greater the divergence rate. In practice, a laser is forced to 
operate in the TEMqq mode by making the cavity optics almost perfect, and by using 
apertures to restrict the lasing to a small diameter. 

1.5.2.3 Longitudinal Modes - For each molecular energy transition, a lasing medium can 
provide one broadened frequency or color at which the laser can, in principle, operate. 
All but one color can easily be removed by using selective mirrors or prisms within the 
cavity. Within the broadened frequency, however, lies a number of individual frequencies 
that must be eliminated to approach a single frequency from the laser. These are denoted 
by the q index on the TEMplq. If the lasing material is placed between two mirrors that 
bound the cavity, standing waves are supported between the two mirrors and their nodes 
occur on the mirrors (Figure 1.19). The condition placed upon the standing wave is that 
the distance between the two mirrors be an Integral number, q, of half wavelengths. If 
the medium is uniform this can be expressed as 

L * q a/ 2 - q c/2f q (1.45) 

fq lies within the broadened energy transition then the laser will support that 
frequency. 
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Figure I.19 Cavity modes and mode selectors. The standing wave supported by the cavity 
is an Integral nuiSber of half wavelengths. 


EXAMPLE: For a typical gas laser, say \ • .6 micrometers, then f • 5 « 10'* sec-1, and 
if L - 1 meter, then q > 3 x 10*. Adjacent longitudinal modes are separated in frequency 
by c/2L or 150 MHz in this case. If the frequency broadening for the given transition 
were 1.5 GHz, then one could expect to get up to 10 longitudinal modes from the laser. 


Figure 1.20 illustrates this effect. The curve la the broadened frequency spectrum, and 
the curve is the product of Che spectrum and the cavity efficiency as « function of wave¬ 
length. The width a is a measure of the selectivity of the two cavity mirrors known as 
finesse, to be diacusaed later. Only those modes within the area above a certain thres¬ 
hold can lase. The various modes are in constant competition with each other for the 
energy of the lasing medium and consequently adjacent modes rarely lase at the same time. 



Figure 1.20 Representation of possible resonant frequencies centered about one color. 






If the cavity contains a varying refractive index, then the wavelength varies 
within the cavity. The nuflfcer of half waves between the two mirrors bounding the cavity 
must still be an integer. Therefore, if the refractive index varies in a direction 
normal to the cavity axis, or if the two mirrors are not precisely parallel, then the 
frequency of light likewise varies across the cavity. This affect is illustrated In 
Figure 1.21, where the cavity supports two different wavelengths on the two different 
halves. This is one reason why it is difficult to produce a single frequency from a 
large diameter cavity. 
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Figure 1.21 Spatial variation of frequency by a nonuniform lasing medium. 

1.5.2.4 Single Frequency Lasers - To improve the laser to its ultimate in coherence 
requires that all but one o? the longitudinal modes be eliminated. This can be 
approached in one of several ways. The laser can be shortened such that the mode spacing 
is larger than the entire broadened spectrum, b. Namely, 

wc/n o L > b (1.46) 

This approach limits the cavity length and, consequently, the output power. 
Additionally, special attention must be given to insure that a mode does Indeed lie near 
the peak of the envelope curve, and that it does not move around under the curve leading 
to instabilities or power fluctuations (to be discussed in the next section). 

Other techniques for frequency selection depend upon the use of tuned reflectors 
(etalons) acting within the cavity or as cavity mirrors. One of the simplest is the 
tilted Fabry-Perot etalon which is inserted into the cavity. The transmissivity of an 
etalon has peaks at multiple wavelengths that depend upon the angle of incidence of the 
wavefront. These maxima are matched with one of the lasing frequencies by tilting the 
etalon. See Figures 1.22 and 1.23. The resonant frequency of the etalon can also be 
adjusted by changing the separation of the reflecting surfaces. The most common way to 
do this is to control the temperature of the etalon. 

1.5.2.5 Stability and Tuning - The latter discussion has been somewhat oversimplified, 
and a few qualifying considerations must be mentioned. From the previous discussion, it 
is clear that a small variation in nL, especially for a single mode laser, leads Co a 
fluctuation in the power output of the laser. Maintaining the optical pathlength, nL, 
between the two end mirrors, exactly constant is not trivial because of thermally induced 
variations in L and because of fluctuations in n in the cavity. 

The variations in n arise out of high frequency fluctuations in the plasma 
properties (for a gas laser), and little can be done to correct that. The long term 
stability of a laser is mostly related to the thermal and mechanical stabllllty in L. 
This is the primary reason why a laser requires a warm-up time to reach a constant 
average output power. 

Tuning of an etalon changes the gain at a given frequency as is shown in 
Figure 1.23. If only one frequency is lasing, it will continue to lase until the 
efficiency for an adjacent mode overtakes it. Then the laser will "mode hop" to the new 
frequency. Therefore, tuning the laser cavity is a discontinuous process. 

The resonant frequency of the cavity may be changed deliberately by changing the 
pathlength nL. The most common method is to mount one of the end mirrors on a piezo¬ 
electric crystal. The insertion of a material of time varying index of refraction into 
the cavity is a second method. 

When a fixed single frequency of light is required, a feedback from the laser 
output can be used to tune the cavity to hold the output frequency constant. Several 
error signals are available for the feedback. One might well conclude that the laser 
output is at a maximum when the etalon is tuned resonant to the center of the broadening 
curve. Then the output of the laser could be monitored and used to drive the cavity 
mirror to maintain this peak. This reasoning is partially correct. The broadening curve 
does not accurately represent the frequency distribution when lasing takes place. 
Saturation effects cause a "hole burning effect" in the frequency distribution, producing 
a gain curve (Figure 1.24) that is characterized by a dip in the curve (termed "Lamb 
dip"). When the laser is operated at this point, either an increase or decrease in nL 
will Increase the output power of the laser. This provides an information source for 
feedback to drive the mirror position. The frequency can be located and fixed by this 
procedure to at accuracy of about one part in 108. 
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Figure 1.23 Tilted etalon characteristic curve. 



Figure 1.24 Representation of lamb dip. 

More accurate methods for location and fixing the frequency of a laser include the 
use of the Zeeman effect to produce a similar ”dip" in the output versus frequency 
curve. The use of absorbing vapors and dyes in the cavity is yet another method to 
narrow wavelength. 

1.3.2.6 (^Switching - In the normal mode of operation, lasing action will begin as soon 
as a certain thresnold of excitation is reached. Lasing will continue as long as the 
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lasing medium Is "pumped". If lasing is inhibited, then the population inversion 
increases with more molecules remaining in the excited state. If the inhibitor is then 
suddenly removed, the amplification is enhanced for a short period of time, and the laser 
will emit a burst of light that lasts as long as is required to sweep the surplus of 
photons from the cavity, typically 10 to 100 nanoseconds. This switching of the cavity 
efficiency or "Q" from one value to another is commonly used in solid state lasers to 
generate high power, short duration pulses of light. 

Q-switches fall into two categories, light valves and beam deflectors. Beam 
deflectors inhibit lasing by misaligning the laser cavity or by deflecting the beam from 
the cavity. Light valves inhibit lasing by introducing an absorber or blocking device 
into the cavity. 

Either of the cavity mirrors can be employed as a Q-swltch by tilting it out of 
alignment to inhibit lasing. Rotating or piezo-electrically driven mirrors are often 
used for this purpose. Acousto-optleal beam deflectors can be used as Q-switches within 
the cavity. 

Electro-optical Q-switches are fast acting light valves. Two common varieties are 
Pockels cells and Kerr cells, which are polarization rotation devices, and, when used in 
conjunction with polarization analyzers can become selectively transparent or opaque to 
polarized light. The polarization rotation is achieved by the addition of a high voltage 
to the crystal, making it possible to turn the Q-switch on and off in fractions of a 
microsecond. Electro-optical Q-switches possess the distinct advantage that they can be 
turned on and off with extreme precision. They can be fired easily with a variability 
less than 50 nanoseconds. On the negative side, they are more expensive and less 
efficient than other types of Q-switches. Frustrated total Internal reflection is yet 
another mechanism used for high-speed light shuttering. 

Some chemicals possess absorption characteristics that make them ideal for use as 
Q-switches. The chemical, which must be highly absorbent at the laser wavelength, is 
placed in a transparent holder within the cavity. This inhibits lasing until the light 
absorbing solution becomes saturated (or bleached) by the Impinging radiation, depleting 
the molecules that are available for absorption. The solution then becomes momentarily 
transparent to the radiation. For ruby lasers, the most common chemicals are crypto¬ 
cyanine dye in alcohol or acetone solution and chlorophyll-D in mineral oil solution. 
The Schott Company RG-690 red cutoff dye filter can also be used as a Q-switch for ruby. 

Chemical Q-switches are the least expensive of all, and, because of their narrow 
resonance absorption characteristics, enhance lasing at extremely narrow wavelengths. 
Unfortunately, they cannot be precisely controlled, have a dependence on environment and 
history, and are not as repeatable as other types. Their variability in time of firing 
is typically about 50 microseconds. Used in combination with other types of Q-switches, 
they can, however, enhance the coherence of a laser. 

1.5.3 Commercially Available Lasers** 7 

1.5.3.1 Summary - Hundreds of different types of lasers are now available to provide 
light at wavelengths ranging from the ultraviolet to the far Infrared. Pulsed lasers can 
provide pulses as small as picoseconds and output powers in the megawatt range up to 
several joules per pulse. Commonly available, continuous lasers can provide powers up to 
tens of watts. The following discussion summarizes the available lasers, with an 
emphasis on the most usable lasers in aerodynamics and on relatively new lasers with 
unusual potential in this field. Table 1.1 summarizes lasers and typical characteristics 
of those of primary interest here. 

1.5.3.2 Solid State Lasers - Ruby and YAG (yttrium, aluminum, garnet) are the most 
common of the solid state lasers in this field. (Diode lasers are treated separately.) 
The dominant wavelengths are 694 run and 1060 nm respectively, and it is possible to 
frequency double either to produce 347 nm and 533 nm respectively. It is further 
possible to triple YAG to produce 353 nm. 

Because of the higher density of radiating molecules in a solid, much more energy 
can be produced per unit volume than in gas lasers. Solid lasers must be pumped 
optically by lamps or other lasers, whereas gas lasers can be pumped by passing a current 
through the partially ionized gas. 

Ruby is termed a three level medium for lasing, because the lasing energy arises 
from a transition from an excited state to the ground state. Initially, the molecules 
are excited to a large band of energy states that decay down to a dominant state, 
creating a population inversion with respect to the ground state. Clearly, this requires 
that at least one half of the molecules in the material be excited before lasing 
threshold can be reached. So ruby is a material that must be pumped "hard" for lasing. 

Population inversion Is much easier to attain between two initially unpopulated 
states. All that is required Is that, during pumping, more excitations go to the higher 
state. Such materials are termed four level lasing materials. YAG is 8 four level 
material, and it has a much lower threshold than ruby. 

Figure 1.25 Illustrates two common ways of pumping a solid state laser. In one 
case, a xenon flashlamp is wrapped around the lasing material; In the second case, a 
linear flashlamp Is placed at one focus of an elliptical reflector while the laser rod is 







Table 1.1 

SIM4ARY 07 AVAILABLE LASERS 


Type 

Medium 

Pulsewldth 

(seconds) 

Energy (Pulsed) 
Power (CW) 
Joules or Watts 

* Principle 
Wavelength 

Wavelength nm 

Coherence: 
Good, Fair 
Poor 

Gas 

Hade 

cw 

IO' 4 - 10*] 

633*, 1150, 3400 

G 


HeCd 

cw 

IO* 4 - 10-’ 

325*. 442 

G 


Argon 

cu 

10-f - 10 

330, 488*, 514* 

G 


Krypton 

cu 

10‘ 2 - 1 

330, 647*. 799 

G 


co 2 

cu 

10* 3 - 100 

9000, 10700* 

G 

Gas 

Argon 

10" 6 - IO' 4 

io- 4 

330 - 514* 

F 


Cu Vapor 

2 x 10*f 

.006 

510, 578 

P 


Au Vapor 

2 *« 10 Q 

.003 

628 

P 


Nitrogen 

10-8 . 10 -9 

up to .03 

337 

P 


C0 2 TEA 

Hr 2 - 1 

up to .01 

9000 - 11000 

F 


Excimer 

10' 7 - 10" s 

up to 1 

282, 308, 351 

F 


(several 






varieties) 





Solid 

Ruby 

10-8 - 10-8 

10-8 . , 0 

347, 694* 

G 

State 

YAG 

10- 4 - 10-8 

10-8 . 10 

265, 530. 1060*. 

1318 F 


Glass 

10' 6 - 10-8 

IO' 3 . 10 

263. 355. 525, 532. F 





1060* 



Alexandrite 

10- 4 - 10-| 

10*8 - , 0 

730 - 790 

P 


GSGG 

IO' 4 - 10-8 

10"8 - 10 

1060 

F 

Semi- 

GaAs 

CW 7 z. 

up to 10“2 

830 - 910 

P 

Conductor 

Lead Salt 

10‘ 7 - 10" 4 

up to 10“° 

3000 - 30000 

F 

Diode 

Tunable 





Tunable 

Ar or Kr 

CU 

up to 6 

250 - 1000 

F 

Dye 

Pumped 

a _ i 

watts 




YAG Pumped 

10-f - 10 * 

up to 0.1 

200 - 4500 

F 


Excimer Pumped 

10'8 . ]0 -3 

up to 0.01 

200 - 1000 

F 


Flashlamp 

10' 7 - 10' 5 

up to 0.001 

200 - 1000 

F 


Pumped 

a c 





N 2 Pumped 

10 8 . IQ' 5 

up to 0.01 

200 - 1000 

F 


Baflectori. B 



Hof lietor 


Figure 1.25 Methods of exciting the lasing medium 


placed at the other focus. Without Q-switchlng, the laser will begin lasing as soon as 
the threshold is reached, and will continue until the flashlamp ceases to maintain the 
population inversion. The process is illustrated in Figure 1.26. For a material like 
ruby, the lasing in this mode is somewhat erratic and the ou put consists of a series of 
spikes that are randomly spaced in time with a mean spacing of about five microseconds. 
For a four level material, the output is much smoother. 

A Q-switch inhibits lasing and the excitation of the material becomes more 
complete. The fluorescence lifetimes in ruby are relatively large (milliseconds), 
allowing the excitation to continue as the flashlamp operates. Q-swltching should be 
accomplished in less than the fluorescence lifetime, since further pumping does not 
effectively excite the material. Figure 1.26 shows that the Q-switch essentially drops 
the laser threshold and allows the stored energy to be swept from the cavity in a short 
time. 
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Figure 1.26 Laser output comparison with flashlamp current 

A laser can be repetitively Q-switched as long as 1) the lasing material is being 
pumped and/or has a sufficient population inversion, 2) sufficient pumping time is given 
between each switching, and 3) the proper temperature and threshold conditions are 
satisfied. Figure 1-27 illustrates the process. The output power will depend on the 
level of excitation of the material. As the pulses are pushed closer than a few hundred 
microseconds, in typical systems, the power in the second pulse will begin to diminish, 
depending on how much energy had been removed from the rod during the first pulse. 
Commercial lasers can push these pulses to about one microsecond separation for a few 
pulses. Repetitive Q~switch4.ng can be done with any of the devices described above. 
Including the chemical Q-switches, but chemical Q-switches are more erratic and harder to 
control than the other types. 

To produce light pulses with closer spacing than one microsecond has, in practice, 
been accomplished by combining the pulses of more than one laser. This can also be done 
by dividing the cavity into more than one independent cavity as shown in Figure 1.28. 
Each part of the cavity can then be independently Q-switched with an arbitrarily small 
spacing. 

High coherence and the TEMqq mode in solid state lasers is attained by aperturing 
the laser cavity to a beam diameter of less than 2 mm, operating the laser at near its 
threshold, and by using wavelength selective dyes or etalons in the cavity. All of these 
devices tend to reduce the total power output of the laser. To obtain both high energy 
and high coherence requires that the relatively weak pulse be further amplified. 
Figure 1.29 illustrates the most common methods for achieving such an output. One or 
more amplifier is added to the basic oscillator. A ten centimeter long ruby rod can 
amplify the intensity by about ten times. A ruby oscillator and two amplifiers can 
produce about two joules of highly coherent light in 20 nanoseconds. 

Among the newest commercially available lasing materials of potential Interest here 
is Alexandrite. This material can be tuned over a small range around 700 nm. Exactly 
how to use this to advantage is still open to the ingenuity of the readers. Another 
material that is showing promise is GSGG (gadolineum, scandium, gallium, garnet). This 
material has many of the properties of YAG and a lower threshold. The lasing wavelength 
is 1060 nm. The major prospect of the material lies in its potential for producing a 
lightweight laser, for example, for airborne applications. Finally, titanium doped 
sapphire is yet another new and promising material. 

1 *5,3.3 Gas Lasers - Much of the discussion of solid state lasers applies also to gas 
lasers. Gas lasers are typically of lower power output than solid lasers. They are 
largely pumped electronically by flowing an electric current through a plasma in a sealed 
tube. The conversion efficiency for converting electrical power into optical power is 














27 


FLASHLAKP 

CURRENT 


EXCITATION 
LEVEL OF 

THE RUBY 


LASER 

OUTPUT 

POWER 




PEAK POWER DEPENDS UPON SIZE AND EXCITATION 
LEVEL OF THE RUBY AND EFFICIENCY OF THE CAVITY 


Figure 1.27 Multiple Q-switching * ruby later 



Figure 1.28 Double-barrel, double-pulsing laser schematic 

commonly less than one percent. The end mirrors can be fixed directly to the plasma tube 
for low power lasers. For higher power lasers, the mirrors are separated from the plasma 
tube, and the ends of the tube are sealed with windows placed at the Brewster angle with 
respect to the optical axis to reduce losses in the cavity. This supports lasing of 
polarized light favored by the Brewster windows, leading to a highly polarized output. 
Q-switching of a gas laser does not significantly enhance power output, but the types of 
light valves described above for Q-switching are used in the cavities to modulate the 
power output. 

Typical power outputs range from milliwatts to watts. Tuned etalons are used in 
the cavities to enhance single frequency output when required. The most commonly used 
lasers are HeNe (633 nm), Argon (408 and 514 nm), and C0 2 (1060 nm). Less commonly used 
are HeCd (325), Kr (330 and 647), and N (337) (see Table 7.1,). 

1.5.3.4 Diode Lasers - Although they are not yet widely used in this field, diode lasers 
hold great promise Bor many applications because of their high efficiency and small 
size. Diode lasers are orders of magnitude more efficient than gas lasers. The typical 
output of a commercial diode laser ranges from microwatts to tens of milliwatts for a CW 
laser weighing only a few ounces. In the pulsed mode, such a laser can generate hundreds 
of watts of peak power in pulses that are 10 to a few hundred nanoseconds in duration. 

The cavity of a diode laser encloses the junction region of the semiconductor 
between two reflective faces. The cross sectional dimension of the cavity is extremely 
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Figure 1.29 Oscillator-amplifier lasers 

small and usually rectangular, commonly about one by ten micrometers. This leads to a 
highly divergent beam because of diffraction, the divergence being its greatest In the 
short dimension direction. So the output beam must usually be conditioned with lenses to 
make it usable. 

In applications where a small light source is required, a diode laser is often an 
ideal choice. One current limitation with diode lasers is their long wavelengths, 
starting at about 780 nm. This requires special detectors and prevents their use with 
photographic emulsions. Recently, diode lasers in the green have been developed, but 
they are not yet commercially available. This development could be extremely important 
for laser diagnostics. In other new developments, phased diode arrays are being produced 
with net output power capabilities of watts CW. 

1.5.3.5 Metal Vapor Lasers - In recent years, copper and gold vapor lasers have become 
commercially available and they fill an Important gap in flow diagnostics. They emit 
light in nominal 20 nanosecond pulses at pulse rates up to a few thousand per second at 
510 nm for copper and 628 nm for gold. This makes them excellent candidates for use in 
high speed photography, schlleren and shadowgraph. Unfortunately, they have extremely 
poor coherence. 

1.5.3.6 Excimer Lasers - The word "excimer" is extracted from the term excited dimer. A 
dimer is an unstable molecule formed at high temperature from two or more molecules such 
as a nobel gas and * halide (for example, xenon and fluorine). Flashlamp pumped excimer 
lasers operate in the ultraviolet at quite high pulsed output powers. With pulsewidths 
of tens of nanoseconds and energies comparable to ruby and YAG, they fill an important 
gap in flow diagnostic applications. They are com only used as a pumping source for dye 
lasers. 

1.5.3.7 Dye Lasers - Dyes are extremely complex molecules, characterized by many 
possible energy transitions. Consequently, when a dye solution is used as a lasing 
medium, many different wavelengths are possible. With a suitable coid>ination of dyes, it 
is possible to produce a laser that is tunable almost continuously across the entire 

spectrum. A dye laser comprises a dye cell through which the dye solution flows, a 

pumping source which may be either a flashlamp or another laser, a wavelength selector 

* which may be a prism or grating, and frequency selecting etalons to further enhance a 

f narrow band frequency output. When another laser is used to pump a dye laser, the 

wavelength of the dye laser can be no shorter than that of the pumping laser. Therefore, 
| excimer lasers are often used to pump dye lasers, since the resulting dye laser can 

\ produce light over a wavelength range from ultraviolet tc infrared. 

} 1.5.3.8 TEA Lasers - TEA is an acronym formed from transversely excited atmospheric, 

v Commonly operated at atmosphere with CO 2 as a lasing medium and in the pulsed mode, these 

■ lasers can fill an important gap where high power, short pulses of light at 1060 nm 

wavelength are required. Pulsewidths range from nanoseconds to microseconds, and powers 
range as high as hundreds of kilowatts. 






CHAPTER 2 

USING LASERS WITH CONVENTIONAL OPTICAL TECHNIQUES 


2.0 BACKGROUND 

Lasers are often better light sources than more conventional ones, but they do 
suffer some disadvantages that Bust be traded off against advantages. An arc source, 
flashlamp or even an Incandescent source can often outperform a laser depending on what 
is required. Therefore, it is important to choose a light source because of its features 
and not simply because it is a laser. 

In this chapter we describe applications of lasers in photography, shadowgraphy, 
schlieren, interferoaetry, and deflectoaetry, emphasizing those applications where lasers 
extend the methods beyond their usual range of applicability. These will Illustrate how 
the unique properties of laser light can Improve space and time resolution and improve 
the ability to separate signal from noise. Coherence, intensity, and controllability 
associated with lasers can be used individually or in combination to achieve such 
improvement. 

Coherence - Temporal coherence implies single wavelength. An event can be 
illuminated with laser light and observed or recorded through an extremely narrow 
bandwidth optical filter that highly attenuates all but the laser wavelength. Therefore, 
experimentation can be extended into cases involving much greater ambient light 
intensity. Additionally, high coherence simplifies both the theoretical description and 
the practical use of interferometry and has resulted in a widespread increase in its 
applicability. On the other hand, spatial coherence aakea possible the production of 
extremely small, high-intensity concentrations of light not previously available in 
optics, significantly Improving the spatial resolution limits of optical probes. 

Intensity - Extremely high-intensity Illumination can be produced from a laser. 
This is, of course, beneficial when high-intensity is needed, but this should not be 
confused with total energy. With few exceptions, lasers are not efficient energy 
sources. Flashlamps, arc sources, induction sources, and other forms of energy delivery 
should be considered equally as candidates where the task is the delivery of energy or 
the production of large quantities of light. 

Controllability - Spatial coherence, temporal coherence, and polarization, all make 
laser light much easier to control in space and time than Incoherent light. For example, 
almost all precision optical elements are easier to design for a single wavelength 
including reflectors, lenses, shutters, filters and modulators. Also, highly collimated 
light producible with spatially coherent light can be transported over large distances 
without losses. Further enhancing controllability are electro-optical modulators and 
shutters that can create pulses of light with great precision in time, placement and 
duration. 

Disadvantages - The problems of speckle, interference, and diffraction noise often 
render recordings made with highly coherent light extremely difficult to use. Unwanted 

diffraction patterns, "parasitic'^ interference fringes, and patches of light can mix 

coherently with data making the result an even more serious degradation than simple 
addition of intensities. In general, therefore, laser optics must be chosen more 
carefully. Scratches and flaws are more deleterious, and cleanliness is much more 

important with coherent light. It has been said that if natural light had been coherent, 
then the discovery of incoherent light would have represented a major breakthrough. 

2.1 Laser Photography and Videography 

Photography and videography in aerodynamics can be divided into the following 

tasks: 

1. Illuminating the object field, 

a. with light of a desired structure, 

b. with the required intensity, 

c. from the desired direction, 

d. for the required time, 

e. at the desired time. 

2. Recording the light that has passed through or reflected from the object field, 

a. at the required time, possibly synchronized with other events, 

b. with the required temporal resolution, 

c. with the required spatial resolution, 

d. while rejecting other light sources that are of no i/iterest. 

Section 1.4 described the various methods available to control coherent light, and 
Section 1.5.3 described the characteristics of available lasers for producing the 
light. This section will show by example how the methods can be combined to achieve the 
above with some rather spectacular results. The most advantageous uses occur when one of 
the above nine tasks is at an extreme. Higher intensities, more complicated structures, 
delivery at larger distances, with greater ambient light rejection, and more precisely 
timed delivery (than can be derived from nonlaser sources) are extremely useful features 
of lasers when photographic limits are approached. 
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2.1.1 Photography of Hypervelocity Objects 

Studies in aeroballiatlc ranges sometimes require accurate in-flight measurements 
of the position, contours, and surface conditions of high velocity models. Figure 2.1 
illustrates a system in use at the AEDC Range G*» 2 exemplary of the techniques of laser 
high-speed photography. In this case, the laser light is first diffused by a ground 
glass to reduce the effects of diffraction. Then It is concentrated by a lens into a 
region of Interest where it finally strikes a background screen. Six cameras, each 
having four lenses that are focused at a different depth in the field, cover the sample 
volume. Narrow-band filters positioned in the image plane remove the ambient light. 


Range Tank Vail 
(3 ■ Diameter); 


Filter 


4- x 5-In. 

Camera 


Ground ( 



Figure 2.1 Front-light / back-light laser photographic system. 

The specifications of the components in the system vary according to test 
conditions. For example, narrow-band filters range from one to 1000 angstrom bandwidth 
depending on model luminosity. Red sensitive films include two major types, Kodak S0243, 
having a resolution of 500 line pairs per mm and an ASA rating of 1.6 and Kodak Linagraph 
Shellburst, having a resolution of 100 line pairs per mm and an ASA rating of 125. In 
the system described here, the higher speed film must be used when F/numbers greater than 
45 are required. 

The lighting procedure used here is intended primarily to produce a front lighted 
photograph of the model in flight. This requires an advance signal to fire the ruby 
laser at the proper time. The first signal, provided to the laser when the model is 
still far upstream, fires the laser flashlamp at about one msec before the model 
arrives. A CW laser beam interruption then identifies the arrival of the model near the 
sample volume whereupon a firing signal Q-switches the ruby laser. Even if the model 
misses the front-lighted volume, it may still traverse a region where it is backlit. 
Clearly, the problem would be made easier if the exact position of the model were known 
apriori. Figure 2.2 includes a series of laser photographs to illustrate the types of 
data that can be produced in such a system as described here. Figure 2.2a shows the 
appearance of a model in flight at various downrange stations showing the change in 
surface condition and contour. 


Example: The 5000 meter/sec model shown will move 100 micrometers during the 20 nano¬ 
second exposure, representing a blurring of the resulting photograph. 


Figure 2.2b shows the laser photograph of a model in which a coolant is ejected from the 
nosetip in flight. Figure 2.2c illustrates the photograph of a model in which the 
background is striped to augment the visualization of flow. 

Other variations of the laser station described above provide for specific diag¬ 
nostic needs. For example, four lasers were cortlned to provide high-speed photographic 
sequences such as shown in Figure 2.3. These are valuable when a rapid time development 
takes place. Another example shown in Figure 2.4 is the combined front-lit photograph of 
the model surface and the schlieren photograph of its flow field. 
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Figure 2.2a» Erosion model 
Bottom figure 


In the 300 meter range. Top 
Is down range. Note eroding 


figure is early In flight, 
nosetip, and Impact craters. 


A 


* 
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Figure 2.2b. Transpiration cooled nosetlp In Che A£DC track G. Note coolant material 
flowing from the noaetip of the model, which Is traveling at 5000 meters/ 
second. The striped background aids in viewing the flowfleid. 



Figure 2.2c. Front-light photography of a hypervelocity 
projectile. 

Figure 2.2 Front light laser photograph of hypervelocity models 
produced in AEDC hypervelocity ranges and tracks. 

2.1.2 Videography of High-Speed Events 

The standard television scanning rate is 30 new pictures per second, a speed that 
has been chosen to satisfy viewing requirements. Such a slow speed has always severely 
restricted the use of TV in high-speed recording. One way to partially circumvent the 
slow-scanning rate is to use a camera that retains a latent image momentarily, produce 
the image with a short burst of light from a strobe or a laser and, thus, stop the high¬ 
speed motion at a still limiting framing rate of 30 per second. New development- in 
videography in the past decade, however, have extended the ability to frame at higher 
rates and to control the entire image recording process in many useful ways. 
Consequently, videography has begun to play an ever-increasing role in optical 
diagnostics. 
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Four laser camera capturing the Impact of a hypervelocity model 
with spherical particles. The model surface is five centimeters 
across and the spheres are three millimeters in diameter. Photos 
were made in A EDC hypervelocity track G. 


i; 

4 


Figure 2.3 







Figure 2.4 Combined front lit md laser schiieren photograph of a 5000 ro/sec model in 

250 torr air. Note the ability to see boundary layer transition on the model. 






The most significant factor is the development and commercial availability of 
charge coupled devices (CCDs). These sensing elements are configured in linear and two- 
dimensional arrays. The typical array comprises 15 micrometer square sensors spaced by 
45 micrometers center-to-center. Linear arrays containing up to 2048 detectors and two- 
dimensional arrays up to about 300 by 400 are available commercially. The arrays are 
exposed in one step, leaving the intensity information locked within each detector. Then 
the array is scanned electronically in a serial fashion producing a string of numbers or 
a set of strings of numbers that describes the intensity distribution. This process 
provides a very natural way to sample and digitize the information in an intensity 
distribution, casting it into a form that is optimum for controlling and processing by 
computer. 

The remarkable development of these devices owes much to the simultaneous develop¬ 
ment of small, inexpensive, but powerful computers. What is a fairly simple control 
process for a small computer today would have been prohibitive in cost and complexity for 
most laboratories only a few years ago. In this sense it is hardly meaningful to 
separate the sensor and the computer, since the computer ultimately takes the string of 
numbers, performs operations on it and then prepares it for presentation on a monitor or 
printer, possibly as an image or in some even more useful form. 

An exemplary application is in the study of the effect of aerodynamic flow fields 
on imagery 3 (aero-optics). This might be encountered, for example, in astronomy or in 
photography from an aircraft. Figure 2.5 describes an experimental setup. A diffraction 
limited point image is focused on a two-dimensional array CCD after having traversed a 
flow field. The camera system in this case was manufactured by the Spin Physics Division 
of Eastman Kodak Company. The entire array can be scanned at rates up to 2,000 frames'* 
per second. The data is stored on video tape and can be played back at adjustable rates 
for viewing. This provides the capability to examine instant replays of relatively high¬ 
speed events. 



Figure 2.5 Conceptual experiment for measuring aero-optic effects. The incoming 
wave from is distorted by the flow field, limiting the attainable 
resolution of an imaging system in the vehicle. 

Figure 2.6 illustrates the resulting data showing how the point source changes its 
shape under the influence of an aerodynamic flow field. As can be seen here, the image 
has been processed by the computer to find the size and centroid of the image auto¬ 
matically. The image increases in size (blur) as a result of the phase modulating 
effect of small scale turbulence and the centroid moves about (jitter) as a result of the 
larger scale turbulence. Models for the aberration of such images have been developed 
and have been partially confirmed in these types of studies. 5 

CCDs offer extremely good versatility in Image evaluation. For example, higher 
scan rates can be attained by scanning a smaller part of the array, trading off data per 
picture for a higher framing rate. The Spin Physics camera can frame at rates up to 
12,000 frames per second in this way. Furthermore, at lower scan rates It is even 
possible to process and display processed pictures In near real time. 
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Figure 2.6 Aero-optical manifestation in an image. Vibration moves the image around; 

large scale turbulence and aberrations distort the image; small scale 
turbulence tends to diffuse and blur it more severely. The time over which 
an image is observed affects the image differently for each case. 


Still higher recording rates and faster shuttering times can be provided by taking 
advantage of short laser pulse lengths and scanning optics. A CCD array can be divided 
into sectors that are exposed sequentially with a pulsed laser that is electro-optically 
scanned between the sectors. Various examples of the lasers described in Chapter One 
could provide exposures of less than 20 nanoseconds and repetition rates up to one 
megahertz. 

The field of CCD/computer imagery is still young and offers much promise for the 
diagnostician in aerodynamics. Also, this is likely to be the source of information 
handling capability that will be needed for the ever increasing amount of data generated 
by optical diagnostic systems. 


2.2 Structured Light Methods 

Of the three classes of lighting chat have been discussed, direct, diffuse, and 
structured, the most elegant is structured light. When structured light is used, the 
sorting or processing of information actually begins when the field is illuminated. This 
is accomplished by virtue of the function that is structured into the light. The 
following are some of the structures that have been employed and will be discussed in the 
following sections. 

t. Thin sheet-of-light, also called laser planing and screening. 

2. Two or more lines of light. 

3. Two or more sheets of light. 

4. Sinusoidal intensity distributions of light. 

2.2.1 Sheet-of-Light Flow Visualization 


The visualization of complex three-dimensional flows is especially difficult with 
ordinary means, because the planes that are deep within the flow field are masked by the 
ones that are closer to tne viewer. This difficulty is partially overcome with light 
sheet flow visualization by illuminating only the plane in which the flow is to be 
viewed. The typical laser sheet-of-light is generated directly from the laser beam by 
passing it through a cylindrical lens. The resulting sheet is typically one mm in 
thickness and diverges in thickness at the same rate as the original laser beam 
(commonly, about one milllradian) and diverges in width at a rate determined by the focal 
length of the cylindrical lens. Lasers have replaced arc sources in most facilities 
using this type of diagnostic tool because it is possible to generate thinner sheets of 
higher intensity light. 

Figure 2.7, showing a light sheet system employed in an ONERA wind tunnel 6 , is a 
typical system geometry. The sheet generating optics is usually mounted on a traversing 
system that allows the recording of adjacent planes in the flow. The flow itself must be 
seeded with light scattering material such as smoke. Light scattered from the plane can 
be recorded in the form of movies or stills depending on the phenomena under study. The 
laser Is usually a CW laser having a power of 100 or more milliwatts. Figure 2.8 is a 
sample set of data from the ONERA system showing the flow field in the wake of an 
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Model at Angle Laser Input 
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Figure 2.7 leaser sheet flow visualization configuration of ONKRA, 

aircraft model. This method Is extremely graphic in its ability to show vortex 
structure, because the scattering particles are actually spun out of the vortex leaving a 
dark region where the vortex is centered. Other flow features also shown by sheet-of- 
light viewing include the transition from laminar to turbulent flow, shock waves, 
boundary layers, and separating flow. 

2.2.2 Particle Image Displacement Veloclmetry (PID) 

Velocimetry, accomplished with the use of tracers in the flow, is one of the oldest 
velocity measuring methods. The technique consists simply of tracking the position of a 
tracer at two or more times and then drawing the associated velocity vector. The use of 
a sheet of light makes this a much more tractable method in almost anv flow. The scatter 
center (SC) must remain within the sheet for the measurement. The sheet is pulsed twice 
or more to record the time separated positions of the SCs or alternatively a camera can 
be shuttered in some other way to produce the two recordings. Suprisingly, the newest 
refinements of this method represent a state of the art in velocity measuring techniques. 
The new method, sometimes called speckle velocimetry is actually a multiply exposed image 
that can be free of speckle, and particle image displacement velocimetry is more correct 
terminology. 

Figure 2.9a from the work of Meynart* shows a multiple exposure image of a laser 
sheet produced in a convective flow with an illuminating apparatus similar to that of 
Figure 2.7. The SCs in this flow are 5 micrometer particles. Ten exposures of 20 msec 
duration separated by 3 second intervals were made. If any part of this photograph is 
illuminated by a small laser beam then a set of Young's fringes will appear in the space 
beyond the photograph. The fringes will be oriented normal to the flow direction and 
will be spaced Inversely proportional to the speed in the flow. It is possible to 
determine the local velocity component in the plane in this manner. In fact, the same 
methods described in Chapter 3 for analyzing interference fringes can be used to automate 
the process. In Figure 2.9b-c Meynart has also shown an alternate technique for 
processing the data. By Fourier transforming the image and spatial filtering, the 
isospeed contours result. 







(a) The light sheet is parallel to the model surface showing a lengthwise cross section of 
the vortex created by the edge of the model. Note the vortex beginning to collapse 
at the extreme right. 





(b) The sheet of light is perpendicular to the model surface. Four recordings are 

superimposed showing cross section of the vortices at various stages of development. 

Figure 2.8 Sheet of light flow visualization of a vortex shedding from a delta wing. 

This technique has been restricted so far to relatively low velocity flows with 
fairly large SCs. It has been used in gas and liquid flows up to a few meters per 
second. Figure 2.10 from the work of Lourenco et al,*° was produced with a 7-watt argon 
laser operating at the 514 nm wavelength and shuttered with a Bragg cell. The tracer 
particles were titanium dioxide, having a diameter range between 10 and 30 micrometers. 
The film used here was Kodak Royal-x having an ASA rating of 1250. Figure 2.10a shows 
the overall view of a seeded vortex. Figure 2.10b is the result of processing a part of 
the image to produce the velocity vectors. This was produced by passing a small probe 
beam through the negative, generating the Young's fringes, processing the fringes using 
techniques described in Chapter 3. and superimposing the resulting vectors on the 
original picture. 
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(a) Multiple exposure photograph showing cross section of the toroidal flow. 



(b) Filtered linage with fringes displayed with reverse contrast, verticle 
velocity component. Contour differences are 15.5 micrometera/sec. 




(c) Filtered image horizontal component, contour differences are 18.8 micrometers/sec. 

Figure 2.9 Laser sheet velocimetry of a convection torus in glycerin. 

A variety of techniques have been reported to improve data quality. For example, 
even though the original figure is produced on a fast, large-grained material, improve¬ 
ment in the quality of the fringes results if a secondary contact print is produced on a 
finer grained material. This tends to reduce the film grain noise and is an illustration 
of a case where the film grain causes a more serious problem with coherent light than 
with incoherent light. Young's fringe contrast also Is improved with the number of 
exposures up to about 10 exposures. 


Figure 2.10a. LlghC Sheet Velocimetry - Observe the four exposures of scatter centers 
in the vortical flow field. 

This technology is quite active and appears to have tremendous potential for new 
application and refinement. For example, the method may be applicable to a large number 
of gaseous flow problems if the proper pulsed laser is applied. The automation of data 
handling should also make this a much more useful method. 

2.2.3 Deflectometry and Hartmann Methods 

As described in Chapter One, when light passes through a flow field containing 
refractive index gradients, a refraction or bending of the light rays occurs in the 
direction of the gradient. This can equivalently be described as a phase shifting of the 
wavefront; the wavefronts are being compressed more in the region of higher refractive 
index. To determine the effect of the flow field on the wavefront, one can measure 
either the amount of deflection of the rays or the amount of phase shift. The former is 
called deflectometry while the latter is called interferometry. 

Deflectometry produces a direct measure of the actual deflection of individual 
bundles or rays of light passing through a flow field. The light structuring to produce 
the individually separable rays can be done in one of several ways. Figure 2.11 
illustrates one such way in which an aperture (Hartmann) plate is placed Immediately 
following the object field. The plate comprises an array of holes ranging in diameter 
from about one mm to one cm. In some cases the apertures are actually lenses in the 
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Figure 2.10b Enlargement of a portion of figure showing velocity vectors in the flow. 

plate that focus the rays of light onto position sensing detectors on the receiving 
side. These constantly monitor the centroid of the light bundle. The emergence of solid 
state CCD arrays and position sensors have made this technique tractable in manv flow 
diagnostic applications. 

The sensitivity of a Hartmann system depends on how accurately the sensor can 
locate the ray centroid and on the distance from the Hartmann plate to the sensor. 
Typically such a system is capable of resolving about one microradian of deflection of 
the light ray. This Is quite remarkable sensitivity. It corresponds to about one- 
fiftieth wavelength increase in pathlength per cm across the flow field. 

If a wavefront is known to have structure that is larger than the aperture spacing, 
the Hartmann system can be used to define the wavefront. However, caution must be used 
in interpreting results when the spatial frequency content of the wavefront is not so 
bounded. This Is a sampling system and is governed by the nuances of sampling theory. 
If the wavefront is undersampled, the resulting data can be plagued with misleading 
results. 

Figure 2.11b describes moire deflectometry. In this esse the light waves are 
structured by grids, strategically placed at one or more locations forming parallel 
sheets of light. If„a grid is placed on either side of the flow field, then without flow 
a linear set of moire fringes will result, 11 depending on grating frequencies and 
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Figure 2.11 HertsAnn and moire sethoda of seaauring the angle of refraction. 

orlencatlon relative to each other. A flow field has the effect of optically distorting 
the grid preceding the flqw field. When the distorted grid Is observed against the 
undlstorted grid, the solre pattern Is appropriately distorted. This arrangement does 
not require the light to be collimated; It requires only that both grids remain In focus 
In the viewing system. 

A different type of deflectometry results If the two ^grlda are placed after the 
flow field. Here, collimated light must be used. The moire pattern forms because the 
light, having been deflected, passes through the two different grids at different places 
and will be blocked by the second grid differently than light that has not been 
deflected. 

In moire deflectometry laser light is usually not advantageous because temporal 
coherence creates significant diffraction noise. It would appear that some of the lasers 
having poor temporal coherence but good spatial coherence such as copper and gold vapor 
lasers might prove useful In this type of application. 

These methods have not seen widespread use as of yet so It Is not clear how 
Important they will become in flow diagnostics. They offer the advantage of adjustable 
sensitivity (by rotating the grids), an extremely broad dynamic range, and utility in 
very strong flow fields. Perhaps with the combination of new vldeography methods, new 
lasers, and new analytical techniques, they could become powerful diagnostic tools. 

2.3 Shadowgraphs and Schlleren Methods 

Shadowgraphy and schlleren methods are the oldest forms of flow visualization of 
phase objects and, consequently, the methods are highly refined.It Is extremely 
difficult to produce with laser light the same cosmetic quality that Is routinely 
produced with conventional light sources. Even so, because shadowgraphy and schlleren 
methods can be performed so easily with lasers and In conjunction with holography and 
Interferometry, they are used extensively In flow diagnostics. 

Figure 2.12 Illustrates the two methods. In shadowgraphy, the light striking a 
refractive Index gradient, is refracted leaving a darkened zone. The recording must be 
made In an out-of-focus plane to achieve significant sensitivity. This Is one of the 
disadvantages of shadowgraphy. When laser light Is used, diffraction effects can be 
especially annoying. Schlleren methods remove the problem of out-of-focus Images by 
enhancing the Image of the refractive Index gradient In a totally different way. The 
collimated wove that has passed through the field of Interest Is focused by a lens or 
mirror onto a filter plane. For the simple knife-edge filter case, the process can be 
modeled simply as the blocking of the refracted ray in the filter plane. Consequently, 
the part of the Image containing refractive Index gradients appears dark In the final 
Image if the direction of the gradient Is normal to and In the direction of the knife 
edge. 

There is a significant difference In the operation of the filter plane depending on 
whether coherent or incoherent light Is used. When coherent light is used, this plane 
represents the two-dimensional Fourier transform of the amplitude distribution In the 
field of Interest. This makes for a wide range of possibilities for performing 
mathematical operations on the input image before it is finally viewed. Some of these 
were analyzed In detail In AGARDograph !#&.* Since they have ostensibly not reached any 
significant use In flow diagnostics (to the author's surprise), they will not be repeated 
here. Such operations include differentiation of the image, Hilbert transformation, and 
phase contrast Interferometry. Regardless of this potential capability, both shadowgraph 
and schlleren methods are almost exclusively used for qualitative work. 
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Figure 2.12 Shadowgraph and schlieren methods of visualizing 
a phase object such as a density wedge. 

2.3.1 Space Shuttle Schlieren System 

The NASA Spacelab Fluids Experiment System (FES) for the study of fluid dynamics 
and crystal growth in zero gravity is equipped with both holographic recording and a 
laser schlieren system. A 20 Mw HeNe CW laser is used and the resulting images are 
monitored and recorded on standard TV. The knife edge is adjustable along the optical 
axis and can be rotated about the axis. This system has been used in orbital experiments 
and will be used as an example here for typical schlieren system interpretation. 

In Spacelab 3 experiments, a triglycine sulfate crystal was grown in zero gravity 
to determine the quality and feasibility of growing crystals in a totally diffusion 
controlled process. (In a gravitational field, the process is heavily influenced by 
natural and forced convection.) By inserting a seed crystal into a saturated solution, 
the process could be regulated by controlling the temperature difference between the 
fluid and the crystal. When the crystal temperature is significantly higher than that of 
the fluid, the crystal dissolves. As the temperature of the crystal is reduced, a 
critical temperature is reached, depending on solution concentration, where the crystal 
begins growing. 

In this case, refractive index increases with concentration of solute. When the 
crystal is growing, the solution immediately next to the crystal becomes depleted and the 
refractive index gradient is positive in the direction normal to and away from the 
crystal face. When the crystal is dissolving, the reverse is true. The schlieren system 
served two major functions: 

1. To determine when the crystal was growing and when it was dissolving and, 

2. To determine if there were any residual convection currents in the cell. 

Figure 2.13 represents the system geometry and the matrix of possible experimental 
conditions. The diagnostics can be optimized by adjusting the knife edge according to 
the experimental condition. The knife edge should be positioned with its edge normal to 
the refractive index gradient so that the refracting ray (which bends in the direction of 
the refractive index gradient) strikes the knife edge. This requires rotating the knife 
edge by 180 degrees when the process changes from dissolution to growing. The growth of 
other faces of the crystal, furthermore, requires a rotation of the knife edge to 
maximize sensitivity for that particular refractive index gradient. 

2.4 Interferometry 

The basic relationships for most interferometers were derived in Chapter One. In 
this section the most common types of interferometry are reviewed and some of the 
applications in aerodynamics are described. In the last section, the effects of a flow 
Cield on a transmitted wavefront were determined by measuring the amount of the 
deflection of the rays of light that define the wavefront. Interferometry determines 
this effect by measuring the variations in optical pathlength through the flow field or 
equivalently bv determining the amount of the phase shift of the wavefront Impressed upon 
the wavefront by the flow field. In Interferometry it is assumed that the bending of the 
rays in the wavefront are negligible, that they pass straight through the field with a 
phase shift only. 
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Figure 2.13 Appearance of schlleren inage for various knife edge 
settings and density gradients. 

Where refractive Index gradients are so great that ray bending la not negligible, 
the data Interpretation near such a region la extremely complex. Some attempts at 
analyzing such data have been made*s but no general solution of the problem has emerged. 

Interferometers find uses In aerodynamics In many more applications than flow 
diagnostics. An extensive treatment of such applications Is outside the scope of this 
work; however, some of the basic configurations are mentioned briefly. 

2.4.1 Mlchelson Interferometer 

The Mlchelson Interferometer Is a widely used instrument which In Its simplest form 
employs a single beam splitter and two mirrors (Figure 2.14). Variations are used In 
phase and distance measurement, measurement of refractive Index, and measurement of 
optical frequency change. 

Vibration Measurement 


Commonly used to examine vibrational stability of stable tables this apparatus can 
be assembled In minutes even with poor optics. If the k vectors are accurately aligned, 
Equation 1.18 Is applicable. Motions of either Mj or M 2 causing a aL ■ a(L 2 - Lj) 
results In a change of phase difference as at the origin, 0, which Is given by 

4S - (2.1) 

Thus, a change In aL by X/2 causes the observed output to pass through an entire 
cycle of intensity. Likewise, a change In the refractive Index, n, in either path causes 
a shift in as allowing the detection of air currents. 

Rotation of either mirror misaligns the k vectors and Equations 1.8 or 1.10 are 
applicable. For example, as previously stated, a 10 second rotation of Mi produces a set 
of linear fringes spaced by one cm (x - .5). Linear fringes can be obtained only If all 
optical elements are flat. 

Testing Optical Surfaces 

If the optical components are flat, one can test mirrors, lenses, and windows. A 
window Is tested for surface flatness by Inserting the window Into either leg of the 
Interferometer. Fringes will be Initially straight. If both surfaces are flat, they 
will remain straight. If the surfaces are parallel, the fringe frequency will remain 
unchanged with or without the window present. The method is useful In testing flatness 
to about A/10 and parallelism to about one second. The Individual surfaces can be tested 
by replacing Mf or M 2 with the surface to be tested. 
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Figure 2.14 Michelson Interferometer. 

Measurement of the Acceleration of Gravity 

A corner cube reflector reflects light antiparallel to itself. The device shown in 
Figure 2.IS has been used to accurately observe the velocity of such a prism in free 
fall. Equation 1.18 applies. Assuming uqj - uq 2 

I(x.z.t) - 2u 0 2 [1 + cos |£s (L 2 (0)-L, - gt 2 /2))] (2.2) 


Example: The position of the prism can thus be determined at any time with an accuracy 
of better than, say 1/20. At some time t' the signal at the photosensor Is oscillating 
with a frequency, f 1 - gt'/i. At t' - t sec, g - 1000 cm/sec 2 , 1 - .5 microns then 
f' - 20 megahertz. It should be mentioned that the Instrument can be considered a 
Doppler measuring device, wherein the light beam In the L 2 leg Is Doppler shifted (In the 
example given, by 20 MHz). As such, Equation 1.7 applies giving an equivalent result. 
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Figure 2.15 Measurement of the acceleration of gravity. 

Measurement of Coherence Length 

Formation of fringes at the observation plane requires coherence of the input beam. 
Fringe contrast is defined for ideal interference as 

c “ 2u oi u 02^( u 01 2+ u 02 2 ) <2.3) 


Associated with every light wave is a coherence length. If the quantity aL Is Increased 
from zero, one will find that the fringe contrast is reduced below that expressed in 
Equation 2.3. Coherence length is commonly defined as f&L)*/2 where (aL)* is that 
difference which reduces C to zero. Coherence lengths vary from a few cm for ordinary 
lasers to many meters for specially designed lasers. 

2.4.2 Parallel Plane Interferometers 

A large variety of interferometers employ two or more nearly parallel reflecting 
planes to add beams that have reflected from the individual planes. Some of these are 
extremely simple and yet useful. Many of these were not practical devices before the 
advent of the laser because they employ unequal pathlengths or spatial mismatching of the 
split beams. 
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A staple MuurtMnc of the coherent length of e leeer can be accomplished with one 
beaa splitter Mounted parallel to the front reflector of the laeer (Figure 2.16). The 
front laeer reflector providea the aecond leg of the lnterferoaeter. When the dlatance D 
la one-half the laaer coherence length the fringe contract will be very poor. The 
lnterferoaeter can alao be uaed for obeervatlon of aaall varlatlona In the product nD. 



Figure 2.16 Parallel plate lnterferoaeter. 


Plane Parallel Plate 


A plane parallel plate (etalon) la commonly employed In a laser cavity to spoil the 
cavity for all frequencies except a discrete set. The net reflectivity of Intensity of 
such a plate having reflectivity R at each face Is given by (assuming a perfect plate, no 
losses) (Figure 2.17) 



Figure 2.17 Geometry for Equation 2.51. 

R . 4R sin 2 [2,nhcos (e ’)!x 1/(1-R) 2 
p 1 + 4R sin 2 [2,nhcoa (8 , )/x)]/(1-R) i 


A plane wave will not lose any energy to reflection If m> - 2nh cos s', where m is an 
Integer. 


Example: e - 0, h * 0.5 cm, \ “ .5 u- Transmission peaks occur at approximately 
every .25 angstroms, for all wavelengths satisfying x - 2nh/m. 


Such a plate can be inserted Into a laser cavity and tilted until Its resonant 
wavelength matches one of the longitudinal modes of the laser. The only losses at that 
wavelength are then cauaed by scattering (not reflection) from the plate, and by the 
drift of the laser wavelength away from resonance. The remaining longitudinal modes may 
not all be eliminated entirely but they are reduced considerably In Intensity below the 
resonant frequency. The coherence length of a laser can be Increased In this manner with 
a very Inexpensive piece of glass; however, commercially available etalons for this 
purpose are generally of very high quality and are expensive. 

The transmissivity of such a plate Is given by 

2 

T - 1/1 + ^ ,t S (2, nh cos »'/*) a (2.5) 

P (1-R) Z 

A typical plot of the transmissivity versus 2»nh coss'/i is sketched In Figure 2.1®. As 
can be seen, as the reflectivity of the surfaces approaches one, the plate becomes more 
and more selective In transmitting specific wavelengths. 
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Figure 2.18 Transmissivity of a parallel plate, 

Fabry-Perot Interferometer 

With the discussion above, the remaining is simply a discussion of hardware and 
terminology, A Fabry-Perot interferometer consists typically of two plane parallel 
reflecting surfaces with an air spacing (n - 1, Equation 2.5) and a variable separation, 
h. The variation in h is usually accomplished (1) thermal expansion of the spacers, 

(2) pressure control in the air space, (3) separators constructed of piezo-electric 
crystals that can be expanded electrically. The familiar terminology for such an 
interferometer includes. 

(1) Free Spectral Range - The interferometer has a transmission peak in every free 
spectral range of frequency. It has a peak transmissivity for many different 
frequencies. From Figure 2.18 we see that this is given by 

Af - c(velocity of light)/h cos e* (2.6) 

This is typically a few gigahertz or a wavelength spread of a few hundredths angstrom. 

(2) Finesse - Usually defined as the free spectral range divided by the width of the 
transmissivity curve at half maximum, or 


Finesse - i —(2.7) 

1 / o-Rr 

Finesse varies from zero to infinity and is a measure of the instrument resolution. The 
free spectral range is essentially divided into F resolution elements. The maximum 
achievable finesse in practice is about 200. 

(3) Peak Transmission - Practical instruments are not loss free. Peak transmission is 
defined as maximum transmissivity, Tp» for a real interferometer. 

(4) Contrast Factor - Peak Transmissivity/Minimum Transmissivity 

To this point we have discussed only parallel Illumination. If a diverging (point 
source) illumination is employed, one observes a set of circular fringes for each 
wavelength. The transmissivity repeats itself each time nh cos increases by one-half 
wavelength. 

Other variations of the interferometer employ spherical reflectors or one spherical 
and one flat reflector. Laser cavity mirrors in this respect constitute a Fabry-Perot 
interferometer. 

2.4.3 Mach-Zehnder Interferometer 

Figure 2.19 represents the geometry of the widely used Mach-Zehnder interferometer. 
Two laser beams are generated by beam splitting. One beam is directed over a path which 
is uniform in refractive index (Path 1) and is subsequently used as a reference for the 
interferometer. The second beam passes through the interest field before it is mixed 
with the reference. The interest field is imaged to an observation plane. The mixing of 
the two waves characterizes the refractive index changes through changes in the inter¬ 
ference fringes. 

Infinite Fringe Interferometry - With no disturbance in the interest field the linear 
fringes observed at the output can be adjusted out by precisely aligning the k vectors 
(assuming high quality optics, windows, mirrors, and beam splitters). This results in 
essentially a single fringe which ia infinitely wide. The common practice in explaining 
a resulting fringe pattern when a disturbance is present is to apply an approximation of 
Equation 1.7 which reduces it to a form similar to Equation 1.18. Assume that beam two 
passes through the interest field. Refractive index gradients change its direction 
resulting in a nonzero value of ag- Equation 1.7 predicts the proper set of fringes only 
if the wave remains planar, since the equation is valid only for plane wave Interference. 
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Figure 2.19 Mach-Zehnder interferometer. 

Wave number two is, however, in general, modulated by the Interest field into an 
extremely complicated wave. Wherever the angle «2 is extremely small, it can be 
neglected if one still accounts for the phase shift. We can apply Equation 1.7 if the 
term kx sin <*2 is replaced by a*(x), since, in general, 02 will be a function of x. 

2 2 

I(x) - u Q1 + Uq 2 + 2u oi u 02 C08 [*♦<*> + A *o) (2.8) 

The term &f(x) represents the increased amount of phase shift suffered by the ray 

passing through the disturbed field. When the interest field is undisturbed 44 (x) is 

zero and a + q is the undisturbed phase difference of the two mixed waves. A+(x) can be 

calculated as 


7 1 

A*(x) - -j 3 / (n(x.z') - n Q )dz‘ 

z 0 

where zq and zj bound the interest field. 


(2.9) 


The integral in Equation 2.9 defines the optical pathlength difference between 
disturbed and undisturbed interest fields. A particular fringe (constant intensity 
lines), therefore, locates contours of equal pathlength through the field. The simplest 
case occurs when n is not a function of z. Equations 2.8 and 2.9 reduce to 


I(x) 


U 01 


+ u. 


02 2 + u 01 u 02 co8 (f i l n (*>' n ol 


z r z ol + A *ol 


( 2 . 10 ) 


In advancing from one fringe to the next, the optical path length increases (or 
decreases) by [n(x) - n^Jaz " X. Since nn and az are known constants, calculation of 
n(x) and application or Equation 1.36 leaas to density p(x). The function G depends 
mildly on temperature and upon the gas mixture. For gases, it is not significantly 
dependent upon wavelength. 


For air G * 1.07 x 10~ 2 3cm 3 . For az * 20 cm, x * .5 microns, a one fringe increase 
represents a density change of 2 x 1017 molecules/cm 3 . This gives the measurement 
resolution in the order 10 ”* to 10“ 2 atmospheres at room temperature. 


The next simplest case is that of the axlsymmetric body, namely n (x,z) * n(x 2 +z 2 ). 
Equation 2.9 reduces to (Figure 2.20) 



Figure 2.20 Geometry for Equation 2.56. 
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4 *(x) 


R [n(r> ■ 


n Q ]r dr 

. 2 , 1/2 


( 2 . 11 ) 


Solution of this integral ia known as radial inversion. The procedure amounts to 
division of the circle into circular shells in which n is assumed a constant. The 
Integral is then solved first for the outer shell, then for the first two outer shells 
then three and so on, (like peeling an onion). A variety of computer programs and data 
reduction procedures are available. 


Figure 2.21 is an interferogram of this type. The high pressure (high refractive 
index) region behind the bow shock of a projectile in flight manifests itself in a fringe 
system that can be analyzed through application of Equation 2.11. 

Finite Frlnae Interferometry - The method juat described is not always useful when the 
net phase shifts in the interest field do not exceed «, it is useful to misalign the 
original k vectors to provide a set of reference fringes. Again, we make an 
approximation that the change in aj is small and apply Equation 2.7 giving 


2 2 

I(x,z) - u Q1 + Uq2 + 2 uqjUq 2 cos k [z (1 - cos a 2 ) 
- x sin a 2 + »♦(*) - 4*o) 


( 2 . 12 ) 


A linear system of fringes exists on the interferogram which, according to 
Equation 2.12 are equally spaced by i/slna 2 when no disturbance (a+(x)) is present. The 
presence of a disturbance changes this spacing by a fraction of a fringe, p 

p(x) - a*(x)/2, (2.13) 



Figure 2.21 Three fringe interferogram of a Cal. 30 cone cylinder at Mach 2.5 in air at 
one atmosphere pressure. 








Figure 2.22 Two cm. diameter sphere M - 10, 76 torr air, 
finite fringe lnterferometery. 

The quantity p(x) Is measured and Equation 2.12 Is applied to Equation 2.9 (or 2.11 for 
axlsymmetrlc fields) to determine n(x) and ultimately p (x). Figure 2.22 Illustrates a 
finite fringe Interferogram applying these principles. 

2.4.4 High Sensitivity Interferometers 

When the Interest field Is an extremely weak phase modulator, steps must be taken 
to either Increase 4s(x) or to Increase the accuracy of measuring It. For example, to 
accomplish the former, the wavefront can be passed back and forth through the field many 
times. A Fabry-Perot Interferometer becomes a powerful detector of phase change, when 
the phase object Is placed In the space between the reflectors. In fact, a laser Itself 
Is perhaps the most sensitive detector of phase shift of all Interferometers when the 
disturbance is placed in the laser cavity. The output power of the laser Is extremely 
sensitive to phase changes within the cavity. This method has not seen significant 
application, however. 

A _more straightforward method Increases the accuracy of the measurement of 
ss(x).i 4 I s One such method measures directly the intensity of the mixed reference and 
test waves (Figure 2.23). We return to Equation 2.8 which has the Intensity modulation 
term 

■ 2u oi u 02 cos [is(x) + 4S 0 ] (2.14) 


If we limit the equation to small a»(x) and set “ -/2 Equation 2.14 can be expanded 
to give 


!„<*) • 2u 01 u 02 l 1 ' k 4*(x>] (2.15) 


and the measured intensity Is proportional to the phase shift Af. Physically this means 
that the output signal Is nulled before the interest field is Introduced. Such an inter¬ 
ferometer Is orders-of-magnltude more sensitive than those described In the previous 
sections. The Interferometer must either be traversed across the field or other beam 
pairs introduced to define the entire function &s(x). 
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Phase Shift 
Element 



Figure 2.23a Subfringe interferometer 


Simple measurement of the intensity will not ordinarily be satisfactory in 
practice. This would require the laser to be noise free and single frequency. The 
problem of laser noise and beating can be solved by detecting the signal with what is 
known as a balanced detector. 

Smeets and George 17 have developed this instrument to a highly refined state and in 
the above reference they describe 14 different optical arrangements for dealing with 
practical problems. Figure 2.23b illustrates one such arrangement. A laser beam is 
passed through a quarter wave plate which shifts, say, the vertical component of polari¬ 
zation by w/2. The beam is split by a Wollaston prism into two normally polarized 
components which are collimated by a lens. They pass through a test cell, one passing 
through a phase disturbance, +. The two components are again combined by a Wollaston 
prism after which they can be expressed in the following way: 

vertical exp(i(wt ♦ * + ir/2)] (2.16) 

horizontal exp[i( w t)] (2.17) 


Quarter Wave 
Plate 


Flow 

Region 


Interference 

Filter 

Diodes 



Prisms 


Figure 2.23b Smeets and George interferometer. 

They are then passed into a third Wollaston prism which is tilted such that its vertical 
direction is 45 degrees with respect to the first two. This separates the single beam of 
light into two components which have a combination of the above two components. The 
upper component can be written as the sum 


exp[i(<jt + t ♦ n/2)J + exp i wt (2.18) 

The intensity at the photosensor is, therefore, 

2 [1 + sin ♦] (2.19) 

The lower component can be written as the sum 

exp[i ( u t + * - */2) ] + exp i wt (2.20) 

with the resulting intensity 


2[1 - sin ♦] (2.21) 

The configuration shown subtracts the two intensity signals providing an electrical 
slgnal. 


i 
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E - 2 sin + 


( 2 . 22 ) 


This method of signal detection is known as balanced mixing. It removes the bias 
term from the intensity, tends to cancel out constant noise terms and, thus, improves the 
accuracy and sensitivity of the measurement. It is the combination of the x/4 plate with 
the last Wollaston prism which produces the opposite signs before the two sine terms 
above and provides for the bias cancellation. 

One variation of this interferometer has been used to detect extremely weak phase 
changes in boundary layers with sensitivities for pathlength changes as small as .01 
angstroms 10 (Figure 2.23c). This allows the observation of subsonic flow disturbances, 
sound waves, and transition from laminar to turbulent flow. In the method shown, the 
signal is derived from light reflected from a surface, an ideal condition for observing 
flow over airfoils. 

2.4.5 Phase Shift Interferometry 

Phase shift interferometry is a relative new technique that solves many of the 
problems in standard interferometry, especially the ones caused by coherent noise. 
Futherraore, the method is readily adaptable to microcomputers and new electro-optical 
components, therefore, simplifying data reduction. Writing Equation 1.8, which describes 
the Mach-Zehnder interferometer, in the following form 

I t (K) * I fl + I b cos (a*(x) + * 0 ) (2.23) 

Let us consider that a controllable phase shift is added to the Interferometer so that +q 
can be accurately set at will. There are three unknowns, I a , 1^, and a$(x). These can 
be solved if three equations are produced. This can be done by setting three different 
values of $q (usually 0 and * 2tj/3 ) and measuring the respective values of Ij(x), 
l 2 (x), and ^(x). The solution of the equations for a$(x) is then 

4 ,(x) - Tan " 1 /3 (I 3 - 1 2 >/(2Ij_ - I 2 - I 3 ) (2.24) 

The three interferograms can be presented to a CCD array which immediately extracts and 
stores the values of Ij(x), ^(x), and Ij(x), whereupon the computer can solve the equa¬ 
tion for a*(x). This process is much faster and more forgiving of noise than locating 
fringe centers to determine a$(x). Also, it is no longer necessary to interpolate 
between fringes, a process that can involve significant error if the sampling is not 
properly done. Of equal importance is the automatic addition of sense to the data when 
the known values of $q are added. Normal fringe tracing methods in contrast to this 
method cannot determine if a fringe shift represents an increase or decrease in optical 
pathlength. 

One drawback to the method is that three interferograms are required for each 
condition. This drawback can, however, be thought of as an asset in providing a way to 
average out noise by taking advantage of the controllability of the interferograms. 



Figure 2.23c High sensitivity laser interferometry for boundary layer transition 

detection. The upper Wollaston prism produces two separate beams that are 
focused to the model surface. Backscattered light from the two points is 
differentially phase shifted by turbulence in the region of the two 
points. The lower Wollaston prism recombines the beams into two beams that 
are adjusted to cancel each other by the feedback system. 
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A second requirement Is chat Che intensities In Equation 2.24 be linearly related to the 
actual Intensities In the lnterferogram. Therefore, an intermediate storage of the 
lnterferograa must remain linear. The only other specific requirement Is that sampling 
be done at least once between each fringe so that the fringe count is not lost. 

The method of phase shift interferometry has not been available long enough for it 
to have become a widely used method in aerodynamics; however, the apparent potential of 
this procedure la almost certain to lead to widespread use as soon as experimentalists 
fully exploit the method. 

2.4.6 Heterodyne Interferometry 

Another relatively new method Is heterodyne Interferometry, which Is essentially 
the most general case for phase shift Interferometry. In this case the fringes are 
continuously shifted in time while monitoring the intensity in the lnterferogram. We 
return to Equation 1.7, which describes the Interference of two beams of different 
wavelength. In general, the Interference pattern Is a set of fringes that sweep over the 
field, such that at any point, the intensity varies at a frequency, s u , the frequency 
difference of the two Interfering waves. The desired quantity for measurement la &y(x), 
which Is manifested as a phase shift on the time varying Intensity at any point In the 
lnterferogram. 

In heterodyne Interferometry, the desired phase difference representing the unknown 
phase object is measured by monitoring the phase shift of the intensity signal sampled at 
each point In the lnterferogram by a small detector. An alternative to scanning a single 
detector is the use of a detector array. The only requirement 1 b that the array density 
be high enough that at least one detector exists between each fringe In the 
lnterferogram. 

There are a number of advantages In making the phase measurement In this manner. 
First, the accuracy with which phase can be determined Is considerably Improved over 
standard fringe tracing methods. Accuracy la limited by the ratio of the diameter of the 
detector to the fringe spacing and values better than one degree or wavelength/300 are 
commonly reported. Secondly, since the signal Is a high-frequency electronic signal. It 
Is more tractable to electronic filtering, noise rejection, and averaging than other 
methods. Finally, Interpolation between fringes la not required since the signal exists 
everywhere In the lnterferogram. 

Even though this method offers tremendous potential. It has not achieved widespread 
use In flow diagnostics at this writing. It seems evident that either this or the method 
of phase shift Interferometry will eventually be Important In flow diagnostics as the 
technology Is developed and passed on to the users. 

2.4.7 Local Reference Wave Interferometry 

A primary difficulty In the use of Interferometry Is the stability requirement 
imposed on the system. In wind tunnel applications this is further complicated by the 
usual requirement of transmitting the reference wave over or under the wind tunnel 
maintaining not only the stability but also eliminating the density gradients In the 
entire path of the reference wave. Local reference wave Interferometry attempts to 
overcome this problem by extracting the reference wave from the data wave that has passed 
through the flow field. 

To be Interpretable quantitatively, an lnterferogram must have a reference wave 
that has a known phase distribution, preferably a uniform phase distribution such as 
planar. This can be done by taking a part of the data wave and spatial filtering all of 
the phase information from It. Figure 2.24 Illustrates how this can be done. The method 
clearly works best when a large part of the flow field Is undisturbed. When this is the 
case, the aperture should be so arranged to pass Into the reference wave only the part of 
the data wave that Is undisturbed. 

The error In this type of interferometry can be estimated by examining the part of 
the lnterferogram from which the reference wave was extracted. If the fringes in that 
region are straight, then the error is small. If they are curving, Chen the error is of 
the order of the fringe variation. 
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Figure 2.24 Local reference wave Interferometer. The reference wave must be spatially 
J filtered by stops or other filters to remove all phase information. 
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If no part of the data wave is undisturbed then spatial filtering must be 
applied. One method of accomplishing this is with the commercially available, Smartt, 
point diffraction Interferometer. In this system the beam is focused to a point after 
having passed through the field of interest. At the focus, a mask scatters (diffracts) 
light from the centermost part of the beam to produce a reference wave while passing the 
remainder of the light as a data beam. The other means of performing the spatial 
filtering is to actually split the wave into two waves, one of which is then passed 
through a pinhole filter to serve as a reference wave. 

Either method is extremely difficult to accomplish in practice without adding 
significant errors into the interferogram, errors that will be completely overlooked if 
the technique is not fully understood. (The interferogram will look perfectly normal). 
If the pinhole or diffraction mask is sufficiently small to produce an acceptable 
reference wave, the alignment is extremely critical and vibration sensitive, and the 
reference wave is extremely dim, requiring similar intensity reduction of the object 
wave. Therefore, the temptation is to increase the pinhole size, producing what looks 
like a better interferogram. The error in the local reference wave interferogram 
resulting from spatial filtering can be bounded by knowing the diameter of the spatial 
filter. A ray of light that misses the center of focus by 1.22 fx/D (the diffraction 
limited focus radius, see Section 1.4.4) represents one wavelength of error in the 
original wavefront. Clearly, unless the reference wave is produced by filtering with a 
pinhole that is nearly the diffraction limited focus diameter, then ambiguities in excess 
of one wavelength can be anticipated in the resulting interferogram. 

An exemplar wind tunnel application of local reference wave interferometry is the 
demonstration by Bachalo*? of such a system in producing high-speed interferometry movies 
of the flow over an airfoil in a transonic flow at the NASA Ames Research Center. The 
same principles have been applied in local reference wave holographic interferometry, 
described in Chapter 3. 

2.4.8 Diffraction as a Type of Interferometry 

Many important features of diffraction phenomena can be analyzed by applying 
principles of interferometry, and one may consider a diffracting edge to be, perhaps, the 
simplest of all interferometers. Consider the aperture of Figure 1.8. Light scattered 
from the edge has a phase difference with unscattered light. For the small angle case 
this difference along a normal-to-the-edge is given by Equation 1.27. The intensity 
maxima or fringes are defined by Equation 1.28. 

Similarly, a circular aperture is an interferometer (see Figure 1.8). When the 
distance, z, along the optical axis is less than r 2 /x (where r is the aperture radius) 
intensity peaks occur along the axis each time 

Z n - r 2 /ox (2.25) 

t-he reparation between the maxima can be determined bv 

AZ m - r 2 mx - r 2 /(a + l)x (2.26) 

At one far field (m • 1) this separation is 

JZ, - r 2 /2x (2.27) 

An interferometer of this type can be used as a vibration monitor or as a linear 

measuring device simply by monitoring the intensity scattered along the z-axis by a 
surface. 

2.4.9 Reconstruction of the General Phase Object from Interferograms 

An interferogram defines the phase of a wavefront after it has traversed the flow 
field along a given projected direction. In general, every different projection produces 
a different interferogram. The task is to reconstruct the phase object from the projec¬ 
tions. This reduces to the definition of the so-called Radon transform. 20 For radial 
symmetry conditions the Radon transform reduces to the Abel equation. 

During recent years, related image reconstruction techniques have been successfully 
applied to several fields such as medicine, astronomy, electron microscopy, nuclear 
magnetic resonance, geophysics and optical interferometry. An extensive review can be 
found in Reference 21. 

There are several difficulties associated with the reconstruction of phase objects 
from their projections: 

(a) No sharp boundaries can be defined for phase objects, (with the exception of 
shock waves). The definition of the reconstruction region is based on 
a priori information about the flow field. In general, the contribution of 
the phase object outside the reconstruction zone must be assumed to be 
negligible. 

(b) Relatively large errors are present in the interferometry data. The optical 
path differences are recorded as interference fringes. The accuracy of the 
fringe data ordinarily may not be better than 1/2 of the fringe spacing. In 
addition, for an unsymmetrlc density field, the number of fringes for each 





projection varies for different view angles resulting in a nonuniform data 
set. The ideal projection data is a set of nonoverlapping, equally spaced, 
parallel rays covering the whole reconstruction region. 


(c) Only the relative phase shift of the projection can be measured by the 
interferogr an data. There is, therefore, a possibility of incorrect 
identification of fringe numbers. This results in a set of "slightly" 

inconsistent data. 

The mathematical reconstruction method employed in such cases should not be overly 
sensitive to the "noisy" or inconsistent data. 

Although the Radon transformation 


t(p,e) ■ // n(x,y)fi[p -xcos(s) -ysin(e)Jdxdy (2.28) 


provides a rigorous solution to the problem of reconstruction from projections, the 
solution is uniquely determined only by an Infinite set of perfect projections. Here, i 
is the optical path length, n is the refractive index, 6 is the Dirac delta function, e 
i8 the projection angle, p is the coordinate along the projection plane, and x and y are 
the coordinates describing the reconstruction region. In practical problems, the 
discrete nature of projection data and the unavoidable measurement errors may result in 
the failure of reconstruction. Most tomography codes are based upon some approximation 
of the application of the Radon transform with practical "tricks" for reducing the 
effects of such errors. 

In this section, the iterative refinement method of least squares solution for 
tomographic interferometry developed by Tan 21 , is presented. The effects of projection 
number, limited viewing angle and the measurement error on the reconstructed image have 
been determined for this solution through numerical experimentation. As an example 
application, the results of the reconstructed density field around the tip region of a 
revolving helicopter rotor blade that was reconstructed from 40 interferograme are 
presented. 

Iterative Refinement Method of Least Square Solution 

In principle, the ART (algebraic reconstruction technique) algorithms are the 
schemes for solving a large system of linear equations. The reconstruction region is 
divided into a square grid (M - m x m cells) and the refractive index within each cell or 
pixel, ny, is assumed constant. The Radon transformation then reduces to a set of 
discrete rinear equations 

l y X w ij<P' e) n ij " l( p. 0) (2 * 29) 

where Wjj(p.e) are weight factors determined from geometric relations. It should be 
noted that only those factors associated with pixels through which the projection ray 
passes are nonzero. For k different projection angles, each having n nonoverlapping and 
equally spaced parallel rays, the number of equations N, is given by N - k x n. 

Equation (2.29) can be written in the form 

CX+L-0 (2.30) 

where, L and X are N and M dimensional vectors respectively, and C is the coefficient 
matrix with N X M elements. Equation (2.30) can be transformed Into a symmetric, 
positive definite, normal equation given by 

C T L + C T C X - 0 (2.31) 

In practice, due to measurement errors and other inconsistencies, Equation (2.30) becomes 

CX+L-R (2.32) 

where R is the residual vector. An approximate solution is sought for which the 
Euclidean norm of the residual vector iRi is minimum. The Euclidean norm is defined as 

iRi - /y R • K\ (2.33) 


The solution of the normal equation gives the standard least square solution of 
Equation 2.30. 

Procedures for obtaining solutions to a large system of linear equations can be 
found in a number of different scientific areas, 22,23 Iterative techniques have been 
widely used. They begin with an initial estimate and repeatedly modify the estimate 
until some threshold condition is satisfied. There are different ways to modify the 
estimate. The effectiveness of a method depends on the matrix character. Considering 
the features of matrix C (large and sparse) and relatively large measurement errors, the 
so-called method of conjugate gradiants is found to be applicable for image 
reconstruction. 22 
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The computer procedures of the conjugate gradient method have been detailed.2 1 

Some of the results of applying the algorithm are presented here. This procedure runs on 

an IBM AT microcomputer, and is limited to 20 x 20 mesh size. 

It can be proven? 1 * that: The norm of the residual vector R decreases with 

increasing iteration number. Thus, in principle, this code is a technique minimizing the 

residual R to the level determined by the inherent errors associated with the computer 
truncation and measurement data. Therefore, the final reconstruction accuracy is 
determined by data noise. Here the errors associated with the measurement of the optical 
path length (fringe lines) are the major attributor to the data noise. 

Tan has performed extensive numerical analyses to evaluate the convergence of the 
procedure to the correct values. This is accomplished by starting with a predetermined, 
exactly defined distribution producing synthesized interferograms• The study has 
evaluated the effect of errors in the measurement from the interferogram, the number of 
projections, and the total view angle. In this way, precisely known errors could be 
introduced numerically. From these studies the following can be concluded. 

• The convergence rate is increased as the projection ray number, N, is increased. 
The corresponding inversion error for a sufficient number of iterations is, in 
general, independent of N. 

• The normal of the residual vector, R, in all cases, unconditionally converges to a 
limit which is a strong function of the input data error. The resulting error for 
the first few Iterations is independent of the input data error. However, for more 
accurate input data, additional iterations result in more accurate results. This 
underlines the Importance of the accuracy of the data set. It also proposes a 
guideline for determining the maximum number of Iterations based on the rate of 
change of the residual. 

• Typically, a set of data can be Inverted to convergence to the correct value to 
within a few percent with five to ten iterations. 

• Limited-look angles, (>30°) do not have a significant effect on the convergence 
rate or the accuracy of the recontructed data. Acceptable results were obtained 
for a total view angle of 30°, Care must be taken to ensure that sufficient infor¬ 
mation is recorded on the interferograms when the view angle limitations are 
imposed. 

An application of this code is presented in Chapter 3.2.2.4. 
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CHAPTER 3 

AERODYNAMIC HOLOGRAPHY 


3.0 BACKGROUND 

A variety of good holography books are available*and, therefore, the author 
makes no atteapt to improve the already existing coverage of basic holography. However, 
for sake of completeness, a short treatment of the holographic process and its 
terminology is presented here. Also a number of important developments that have not 
reached today's textbooks are incorporated here. This chapter emphasizes the application 
of holography in aerodynamics, especially the application to particle and flow 
diagnostics and the analysis and Interpretation of holographic data. 

Holography has satisfied a need that is so basic to the study of particle and flow 
fields chat its use has continued to grow steadily since the first application. Even 
though the data handling problem has not been entirely solved, the number of applications 
requiring the features of holography has continued to fuel development and use. With 
present technology, holography is the only method for accurately recording a three- 
dimensional image of a dynamic event. This is particularly useful, for example, in the 
microscopic examination of combustion and explosive events, droplet breakup and formation 
mechanisms, and velocity measurement of droplet fields. Holography has likewise made 
some types of flow visualization possible that have not been achievable by other methods. 
Specifically, the applications of interferometry to directly compare two flow fields and 
applications of interferometry in facilities with poor optical windows are outstanding. 

Holography was invented in 1947 by Dennis Gabor 8 who produced holograms of a 
variety of microscopic samples using partially coherent light since lasers were not 
available for another fifteen years. His application exploited the magnification 
properties of holograms and not the three-dimensional imaging properties. Another twenty 
years passed before holography saw application in a field environment In the studies of 
fog droplets by Thompson and coworkers . 9 ' 10 

The first application to flow diagnostics was by Brooks et al . 1 1 , who used 
holographic interferometry to analyze the flow field around projectiles in flight. Since 
then, virtually every test facility in the world with a need for flow diagnostics has 
incorporated holography. One of the first major test facilities to incorporate 
holography into an operational wind tunnel was the USAF Arnold Engineering Development 
Center. *2»13 By the early 1970s, systems were in operation in most of the other 
aerodynamic test facilities of the NATO countries . 14 

The development of data handling techniques lagged far behind what was required to 
compete with conventional photographic techniques. In recent years, great strides have 
been made in removing this obstruction. Available computer power, new codes, and 
detector arrays have played a key role. 15 

Flow visualization holography has not replaced more conventional methods, but 
rather has become a complement. When a distinct set of conditions or requirements exist, 
flow visualization holography is often the best or least expensive or, quite often, the 
only way to obtain a desired set of data. It is otherwise usually, but not always, the 
more difficult method since it involves less refined, more complicated hardware. 

3.1 Introduction to Holography 18 

3.1.1 How Holograms Work 

Ideally, a hologram of a scene can be thought of as a window into the past, through 
which a precise image of a scene can be viewed even after the actual scene is physically 
removed. Having the hologram is, from a viewing standpoint, equivalent to having the 
actual scene behind the hologram frozen in time. The radiation field on the viewers side 
of the hologram has the same information content in either case. This radiation field, 
which is unique to a given object field, can be extremely complex. It is synthesized by 
illuminating the hologram with a simple wavefront, which, through the process of 
diffraction, is converted into the more complex one. We are referring to the coherent 
radiation field which originally was reflected from or passed through the object field. 
Radiation generated by the object field must be excluded, since it is not coherent and it 
simply adds a bias to the coherent radiation; it is not involved holographically. 

Holograms can be made by many methods on many types of material. They can be 
generated by computer, drawn by hand, or produced by photographic methods. The 
discussion here is limited to optical holograms. A hologram can be made by recording the 
intensity of the sum of two or more waves, one of which must *>e a reproducible wave (the 
reference wave), and the other of which is the modulated wave carrying phase and 
amplitude information (the object wave). 

Holography will be explained both as a diffraction process in general and as a more 
simple reflection process for some cases. A hologram is a complex diffraction grating 
that diffracts the light of a reconstructing wave to form an image of the object field. 
Holograms can be understood most easily in terras of point objects. Once the hologram of 
a point source is understood, then more complex objects can be explained in terms of a 
superposition of many point holograms in a linear process. 
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Figure 3.1 Illustrates the foruatlon and reconstruction of a hologram of a point 
object formed with a collimated reference wave. Note that as the angle between the 
object and reference wave becoues larger, the spatial frequency of the Interference 
fringes likewise Increases (see Chapter 2). When the developed recording Is relllumln- 
ated, the wave la diffracted wore by the higher frequency diffraction grating. The 
diffracted rays of light converge on a point to the right of the recording, and they also 
diverge froa a point on the left aide of the recording. The latter Is called the virtual 
Image, while the former Is called a real Image and la said to be the conjugate Image of 
the latter. Because diffraction gratlnga exhibit this dual directional character, 
holograms always possess two Images, one of which Is the conjugate of the other. 

Figure 3.1b Illustrates this process in terms of reflection. Recalling that fringe 
planes form as a bisector of two Interfering waves, with a spatial frequency that 
increases with angle between the two waves, an emulsion will be exposed to light sheets 
as shown. A processed silver halide emulsion will contain deposited silver layers In 
these regions that will act as mirrors to a reconstructing wave, diverging the llRht from 
a point where the original object point was located. This does not explain the simul¬ 
taneous reconstruction of the conjugate point; however, a real Image can be observed to 
occur from this process if the hologram la illuminated from its opposite aide. 



Figure 3.1a Diffraction Model. Note that the fringe frequency Increases as the angle 
between the two Interfering beams Increases In forming the hologram. In 
reconstruction, the higher frequency grating diffracts the light through 
larger angles. 




Figure 3,1b Geometrical Optics Model. The fringe planes create planar deposits of 

silver in the emulsion that act like tiny mirrors to reflect the light at 
different angles during reconstruction. Note that the fringe planes bisect 
the angle between the object and reference waves. 

Figure 3.1 Holography of a Point 
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Figure 3.2 defines a general geometry for forming and reconstructing a hologram of 
a point object. A first order analysis for thin holograms, in which the angles between 
object and reference wave are small, results in the following set of equations for 
locating the reconstructed image. 
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When the object and reference waves are collimated, the holography equations become 
simply the relationships between angles. 
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Figure 3.2 Geometry for the formation and reconstruction of a hologram. 

3.1.2 Classifying Holograms 

Holograms are classified according to the recording material placing relative to 
the object and reference wave, the direction of the lighting of the object field, the 
type of lighting, and the use of lenses In the recording. Figure 3.3 summarizes this 
taxonomy. Two basic types of holography are ln-llne (made with the object wave almost 
parallel to the reference wave) and off-axis (made with a relatively large angle between 
the object and reference wave), as illustrated In Figure 3.4. 

An in-line hologram can be made by passing a simple wavefront (plane or spherical) 
through the field, which must be about 80 percent transparent, and recording the inter¬ 
ference pattern of the scattered (object) and unscattered (reference) waves. Off-axis 
holograms are made by mixing the scattered light with a mutually coherent wave that has 
taken a path around the field. Off-axis holography is much more versatile with respect 
to the type of object illumination. 

HOLOGRAPHY LIGHTING LIGHTING RECORDING 

TYPE DIRECTION TYPE TYPE 


> 



Figure 3.3 Hologram Classification. Shown here are twenty-seven 
different classifications of a hologram. 
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Off-Axis Holograa 




(a) In-Line and Off-Axis Hologram of a Point Object (b) Forming the Image with the Hologram 
Figure 3.4 Holography of a point object. 

An object field can be front lit or back lit, with direct, diffuse, or structured 
light. The most common type of structured light is the so-called sheet of light produced 
with a cylindrical lens and a small diameter collimated bean of light. Direct lighting 
is derived by expanding or diverging light directly from the laser. 

Imaging with lenses or mirrors before recording is commonly done to impart a 
convenient size, location, or magnification to the field of interest, or to relax the 
required recording capabilities of the hologram. This is discussed in more detail later. 

Other classifications describe the process of diffraction, the character of the 
recording material, and how the information is coded before recording. The diffraction 
process can be caused by amplitude transmission variations (amplitude hologram), phase 
shifting (phase hologram), or by pathlength difference on reflection (surface relief 
hologram). The process can occur in approximately a single plane (thin hologram) or In a 
volume (thick or volume hologram), and it can occur on transmission (transmission holo¬ 
gram) or on reflection (reflection hologram). More than one hologram can be stored on a 
single plate (multiplexing), or they can be stored in narrow strips by coding with 
cylindrical optics (integral holograms) or coded by filtering one of the dimensional 
dependences from the object wave (rainbow holography). 

Clearly, a large number of possibilities exist. All of these have found use in 
some form in aerodynamic applications and there are many variations yet to be tried. 

3.1.3 When to Use Holography 

When the evaluation of diagnostic methods for a particular problem includes 
holographic techniques as candidates, one should ask if holography offers more than other 
recording methods. Holography is essentially an information buffer memory between the 
experiment and the data analysis. Holographic data storage often makes a standard 
technique applicable that would not otherwise be; when conventional Imagery cannot store 
sufficient optical information, holograms sometimes can. 

How many photographs are stored in one hologram? A resolution element, R, can be 
stored in a photograph if the element lies within a field depth R 2 /2>i near a single 
focused plane (see Figure 3.5). A hologram with effective diameter D can resolve R if 
the element lies within a distance of RD/ X of the hologram. The ratio of these two is 
the number of photographs in the hologram, N p ,clearly a very large number. 

Np - 2D/R (3.3) 


EXAMPLE: A 10 cm diameter hologram is produced to cover a field of particles having 

5 micrometer minimum diameter. This hologram could contain up to 40,000 different, 
useful photographs. 


Holography is especially helpful in dynamic cases where it is reaulred to record 
high resolution optical information ever a large volume in a short time. But, in 
stationary cases, holograms may not necessarily be superior in storing the 3-D data, 
depending on how it is to be processed later. When a wavefront is to be analyzed for 
phase information, holographic storage is convenient since the hologram retains all 
optical information including phase. (But holography is not the only way to store phase 
Information.) When the wavefront is extremely complex, such as that from a diffuser or 
from very poor optics, holography is likely to be the best if not the only candidate. 
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Figure 3.5 The sample volume in a hologram. The number of different photographs in a 
hologram can be estimated by dividing the photographic depth of field at a 
given resolution into the depth over which the hologram can achieve that 
resolution. 

When a variety of wavefront diagnostics is desired, holography provides a method to 
store the wavefront until the diagnostics stage takes place. This is especially true 
when the experiment time is either short or extremely expensive, since data reduction can 
be done at a separate time and place at the scientist's convenience. 

Finally, when ambient light levels are too great for photography, holography can be 
used to extend the recording range. The hologram acts like a coherent light filter. 
While ambient light does fog the hologram, it plays no direct part in the reconstructed 
wave. The reference wave intensity can be raised to greatly exceed ambient light level 
at the hologram. When mixed with the object wave coherently, it provides a tvpe of 
amplification. 

3.1.4 Hologram Quality 

The quality of a hologram can be defined by three parameters: 1) diffraction 
efficiency, 2) resolution or spatial bandwidth, and 3) signal-to-noise power ratio or 
smallest retrievable object intensity. Diffraction efficiency depends upon the 
holographic fringe visibility, which is maximum when the object and reference waves are 
of equal intensity. (This is not where S/N is largest.) In some materials, diffraction 
efficiencies are above 90 percent, resulting in extremely bright holograms. Brightness 
is an important parameter in holographic gratings and in display, but it is of little 
importance in diagnostics. 

Resolution is determined by experimental geometry, hologram size, and quality of 
optical components. The latter requirement can be replaced by a capability to accurately 
reproduce the reference wave used to make the hologram, as well as to reproduce the 
optical geometry in which the recording was made (to subtract out aberrations during 
reconstruction). 

Signal-to-noise ratic (S/N) is determined by the type of emulsion and by the qual¬ 
ity of technique such as cleanliness, linear recording, using liquid gates, and refined 
processing methods. Table 3.1 summarizes the techniques for improving resolution and 
signal-to-noise. 

3.1.5 Recording Materials 

Three types of recording materials are important in aerodynamic holography, silver 
halide emulsion, thermoplastics, and photochromies. The highest quality holograms are 
still made on silver halide emulsions, but the processing requirements of these materials 
greatly restrict their use. Thermoplastics and photochromies have been used only in very 
recent years, but they offer tremendous advantages to be described. Other types of 
recording materials that have seen little use in this field will not be discussed. A 
detailed discussion of recording materials is included in Smith's book. Holographic 
Recording Materials. 16 
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Table 3.1 



TECHNIQUES FOR IMPROVING 

HOLOGRAM QUALITY 

WHERE 

SIGNAL-TO-NOISE 


APPLIED 

RATIO 

RESOLUTION 

DESIGN 

1. Use as few optical elements as 

1. Hologram size must be large 


possible. Low scatter elements 

enough 


2. Baffel stray light 

2. Emulsion resolution high 


3. Spatial filter beams 

3. Low object/reference beam angle 


4. Choose an emulsion with low 
noise, high MTF, high Gamma 

5. Have reference wave at least 

4. Use flat plates 

5. Place object close to hologram 


10 times object wave 

6 . Use direct (not diffuse) light 


6 . Object reference angle is 

7. High quality optics (flatness. 


bisected by plate normal 

parallelism, good lenses) 

8 . Image and magnify the object 
before recording 

9. High coherence 

DURING 

1. Keep everything clean in 

1. Keep everything stationary to 

RECORDING 

the optical train 

A/I 0. 


2. Maintain uniform intensity 

2. Keep windows thin and go through 


wavefronts 

3. Put as little as possible in 
the object field 

4. Expose optimally 

them straight 

DURING 

1. Develop to transmissivity of 

Same as for high S/N 

PROCESSING 

.16 - .20 

2. Shorten development time 

3. Maintain cleanliness 

4. Final distilled water wash 


DURING 

1. Keep optical train free of dust 

1. Orient hologram position 

PLAYBACK 

2. Uniform beams 

precisely 

2. Use real image corrected by pro- 


3. Use a liquid gate 

4. Illuminate only the part 

jecting backwards through any 
optics originally used 


of hologram required 

3. Use correct wavelength 


5. Spatial filter 

4. Duplicate original reference 
wave 

5. Spatial filter 


3.1.5.1 Silver Halide Emulsions - Figure 3.6 summarizes the properties of a typical 
photographic emulsion. The optical noise scattered from a beam of liRht incident upon 
the emulsion is greatest in the transmission direction, and the ratio of transmitted to 
scattered light begins to level off at about 30 degrees at a value of about 10” . Note 
that the diffraction efficiency peaks at an intensity transmission coefficient of about 
0.16, a number which has been found to apply to most photographic materials. The amount 
of noise scattered from the illuminating beam decreases with photographic density since 
it is absorbed in the emulsion. The diffraction efficiency peaks where the slope of the 
transmission versus exposure curve is maximum. 

Figure 3.7 compares the characteristics of a variety of photographic materials. 
Note that the peak diffraction efficiency is about 3 percent. The fastest emulsion is 
Kokak SO-253, which is also the lowest in resolution. The most commonly used materials 
for ruby laser holography are Agfa 10E75 and 8E75, while Kodak 120-02 is most widely used 
for HeNe laser holography. Kodak 649-F, first used for holography, produces excellent 
quality holograms but is much too slow to be of use generally. 


4k 
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Figure 3.6 Characteristics of silver halide emulsions. 

DIFFRACTION 
EFFICIENCY l 



Figure 3.7 Characteristics of common holography materials. 

In producing the hologram, we are interested in four things related to the film: 

(1) maximum diffraction efficiency for a bright image, 

(2) high signal-to-noise ratio, 

(3) image resolution for high definition, and 

(4) minimum recording time. 

The diffraction efficiency is simply the percentage of the reconstruction reference 
which is diffracted into the reconstructed object wave. The sources of noise include 
scattered light from the emulsion grain, surfaces of the emulsion backing and other 
surfaces, and contamination in the system. Parameters that affect desired hologram 
characteristics include the exposure of the hologram, its development, and the choice of 
a reference to object beam intensity ratio. 

(1) Exposure - For most films, the maximum diffraction efficiency occurs when the 
amplitude transmittance at the average exposure level lies between .4 and .5 for all 
object to reference wave beam intensity ratios. This corresponds to an intensity 
transmittance of .16 and a photographic density of approximately .8. The diffraction 
efficiency can easily vary by a factor of five over the hologram if the exposure of the 
hologram varies by a factor of five. 

(2) Development Time - A vast amount of control is possible during development. 
Normally, a shorter development time Improves the signal-to-noise ratio, since the 










64 


diftract ion efficiency increases slower Chan noise beyond normal development time. This 
means that it is better to overexpose the hologram in the first place and decrease 
development time to improve the signal-to-noiae ratio. Also, this leads one to the 
conclusion that choosing a higher resolution film does not necessarily improve the 
signal-to-noise ratio if too long a development time is required. For example, one might 
otherwise choose Agfa 8E75 over Agfa 10E75, since the resolution of the former is 
higher. The conclusion here is that such a choice produces an improvement only if the 
development time remains short. 

(3) Reference to Object Beam Intensity Ratio - The simplified analysis of the 
holography process usually makes the assumption that a linear relationship exists between 
the exposure energy and the amplitude transmission of the developed hologram. When such 
a relationship does not exist, the recording process is said to be nonlinear. Non- 
linearly recorded holograms produce not only the normal first order holographic image, 
but also higher order Images. These images, in themselves, do not normally constitute a 
problem. However, further analysis shows that even first order holographic images in 
nonlinearly recorded holograms have added noise terms. Therefore, this results in a 
lower signal-to-noise ratio. 

A wide linear range can occur in holography processes when the reference to object 
ratio is 10 to 1 or greater. When the reference to object ratio is less than this, care¬ 
ful development time is required to keep the recording linear. The brightest hologram is 
not necessarily the linear recording. In fact, the highest diffraction efficiency occurs 
when the reference to object ratio is equal to unity, so brightness is not the correct 
thing to optimize. Therefore, the hologram which has been optimized for resolution is 
normally not as bright as the hologram maximized for brightness. 

Development time should be shortened to achieve the photographic density on the 
hologram of .8 by increasing the exposure during recording, and by carefully monitoring 
the hologram during processing. The normal way of processing a hologram is to develop 
and inspect the hologram occasionally during the development, using a suitable safety 
light. This practice can be made more accurate by comparing with a properly developed 
hologram placed side by side with the hologram under processing. A more scientific and 
less subjective way of making this decision employs an infrared light emitting diode and 
receiver, inserted into the development tank. As the hologram develops, constant 
monitoring of the transmissivity can be made without effecting the hologram emulsion, 
which is not sensitive to the infrared radiation. 

3.1.5.2 . Thermoplastics - With the exception of lasers, the thermoplastic recording 
device (TPD)is probably the most important development for the field of applied holog¬ 
raphy. TPDs provide a method to record and develop holograms in place electronically. 
Figure 3.8 depicts the principle of operation. A transparent conductive film on a 
transparent substrate is coated with a photoconductive film and a thermoplastic layer on 
top. 


CHARGING 


SECOND CHARGING 


'////////////////, 


3 — 


THERMOPLASTIC 

PHOTOCONDUCTOR 


A UNIFORM CHARGE AS 
DEPOSITED ON THE THERMOPLASTIC 
PHOTOCONDUCTOR. 


EXPOSURE 
LIGHT LIGHT 


V///////////////,\ 


l^ r . — I 


'////////////////A 




RECHARGE INCREASES THE ELECTRIC 
FIELD ACROSS THE EXPOSED AREA OF 
THE THERMOPLASTIC. 


DEVELOPMENT 




HEATING OF THERMOPLASTIC 
CAUSES PERMANENT DEFORMATION. 

THE IRREGULARITIES UILL DEFRACT 
LIGHT TO RECREATE THE ORIGINAL IMAGE 


EXPOSURE TO LASER FRINGE 
PATTERN REDISTRIBUTES THE 
CHARGES THROUGH THE 
PHOTOCONOUCTORS. 


ERASURE 



ERASURE IS DONE BY 
CONTROLLED HEATING. 


RECORDING-ERASURE CYCLE OF 
THE THERMOPLASTIC HOLOGRAM. 


Figure 3.8 Thermoplastic Recording Principle of Operation 
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To make a recording, a positive charge layer la added to the thermoplastic. 
Exposure to light redistributes the charge ao that the capacitive forces compressing the 
therooplastlc are greatest where the exposure occurred. A second charging enhances this 
effect even more. Heating the therooplastlc allows It to flow under the electrical 
forces. Finally, cooling the thermoplastic causes it to set permanently In Its deformed 
state, resulting In a phase hologram. The entire process takes a few seconds. 


Figure 3.9 shows the characteristics of typical TPDa. The sensitivity of this 
material is similar to that of Agfa 8E7S emulsion; however, the resolution Is, In 
general, worse than that of silver halide emulsions. The diffraction efficiency is 
considerably higher since a phase process Is Involved. Since the material Is Insensitive 
before charging and after development, It Is easier to use. 



mm 



EXPOSURE EXPOSURE 


Figure 3.9 Characteristics of typical thermoplastics 

TPDa suffer some disadvantages, however. A high voltage process Is Involved In the 
charging of the plate, the materials are sensitive to moisture, dust, and temperature, 
and the Initial cost of a system Is relatively expensive. The available formats are 
limited to 35 mm and 70 mo. Probably the most serious limitation is the relatively low 
slgnal-to-nolse ratio compared to most photographic materials. 

Two system types are commercially available. One of these is a glass backed 
erasable material that allows the recording, erasing, and reusing of the plate up to 
about 1,000 times. This can be extremely convenient and cost effective, If the required 
data can be extracted from the hologram shortly after recording. The second type has a 
flexible backing and Is not erasable. This Is convenient and cost effective if each 
hologram la to be saved. 

3.1.6 Holocameras and Reconstruction Systems 

3.1.6.1 General Requirements - A holocamera is a Bystem of lenses, mirrors, lasers, and 
other optics for making holograms. In aerodynamic holography, holocameras are specially 
designed to meet the requirement. Choosing from the array of possibilities described in 
Section 3.1.2, the system must be engineered to accomplish the following tasks: 

(1) Properly Illuminate the field of Interest. 

(2) Place the interest field or Its Image In a suitable place relative to the hologram. 

(3) Achieve an acceptable F-number for the desired resolution. 

(A) Deal with environmental factors. 

(5) Meet coherence requirements. This is done by splitting the laser beam into two 
beams, using one as a reference wave and the other as an object wave, and bringing 
the same rays back together as accurately as possible after having traveled the 
same distance. 

(6) House the film or plate transport, laser, and optics. 

(7) Produce a hologram of acceptable quality. 

The reconstruction system la the system of optics that provides for the replay of 
the hologram Into a data reduction system. It may be part of the holocamera or It may be 
a completely separate system. Ideally, the hologram should be illuminated with the same 
wave as the original reference wave, but practical considerations often prevent this. 
So, when a hologram is made with a pulsed laser, a CW laser Is commonly used In the 
reconstruction system, and It Is not always possible to match wavelengths.’ 
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Changing the laser wavelength, angle of illuminatton, and radius of curvature for 
hologram reconstruction introduces aberrations Into the final image that must be taken 
into account. Such aberrations are similar to the Seidel and chromatic aberrations 
present In lenses. A hologram that is illuminated at the wrong angle will exhibit 
astigmatism and coma, while one that is Illuminated with the wrong radius of curvature 
exhibits spherical aberration. An aberrated reference wave translates with a linear 
transfer function directly into the reconstructed image. 

A reconstruction system should Incorporate the following elements: 

(t) Optics to produce a reference wave like the one with which the hologram was made. 

(2) Mechanical devices to position the hologram in the reference wave as it was 
originally recorded. These should be adjustable over few degrees. 

(3) A scanning system to scan through the three-dimensional image. 

(4) Controls that allow the viewer to scan the hologram while examining the image. 

(5) A recorder, such as a closed circuit TV and VCR or a photographic camera (or both). 
Some reconstruction systems require special features such as the following: 

(6) Adjustments that allow the positioning and location of the hologram to wavelength 
accuracy. 

(7) Fixtures to hold and position two or more holograms relative to each other. 

(8) Allowance for more than one reference wave. 

(9) Independent phase and intensity control on one or both of the reference waves. 

(10) Computer controlled traversing and other mechanical adjustments. 

(11) Features that allow the reconstruction of the conjugate wave backwards through the 
original recording optics for the purpose of aberration removal. 

The reconstruction system must be compatible with the holocamera and must be 
designed to produce an image that is sufficiently bright, noise free, aberration free, 
the right size, and in the right location for analysis. An extremely important function 
of this system is to cast the data into a useful form. The images may be recorded as 
photographs, or analyzed directly by an image analyzer. So the quality of the optics 
must be generally as good as the required image quality. 

3.1.6.2 Typical Holocamera - Figure 3.10 illustrates a common holocamera design, using a 

pulsed ruby laser as a recording laser and a collinear HeNe CW laser as an alignment 
laser. The CW laser can also be used for reconstruction within the holocamera. This 
particular system uses a double-barreled laser to provide two separately Q-switched beams 
of light. (Most holocameras use only single beam lasers.) Note that the object and 

reference pathlengths are matched. The object field is imaged with a one to one imaging 

telescope to a position near the hologram. In a system like this where the same rays do 
not come back together at the hologram, it is important that the lasers have good spatial 
coherence. This type of holocamera is capable of producing three-dimensional images of 

the object field with a resolution of about 10 micrometers. Multiple pulsing of the 

laser and use of both beams make possible the recording of multiple 3-D images, with time 
spacings from almost zero up to about 800 microseconds. 

A ruby laser equipped system can make holograms typically with a 20 nanosecond 
exposut', which freezes the motion of most fields of Interest. Other laser candidates 
for this type of system are YAG lasers and dye lasers. Ruby lasers were historically the 
first used, but YAG lasers are gradually becoming more competitive. Compared with the 
maximum one multiple pulse per second of commercial ruby lasers, YAG lasers can pulse at 
rates up to about 30 per second, making them ideal for making movies. 

3.1.6.3 Viewing Reconstructed Images - A key element in reconstruction and analysis of 
holography data Is the development of skill and technique in focusing coherent images. 
This can and should be developed upon images which are first generated directly without 
holograms, so that holographic effects can be understood separately. The holocamera in 
Figure 3.10 has a built-in alignment laser that projects a continuous wave helium-neon 
beam of light through the system in the same space that the ruby laser will ultimately 
transcend. This alignment laser can be used for more than alignment, providing informa¬ 
tion about the reconstructed image, even before the hologram is recorded. This is done 
by placing a stationary object, in the field of view that will ultimately be occupied by 
the experimental object, and examining the image that is to be presented for holographic 
recording. This image can, by removing the frosted glass from the hologram holder, be 
projected into a working space which is directly accessible to the researcher. There¬ 
fore, many of the properties of the image can be studied in detail directly, without 
being contused by holography effects. 

It is useful, especially for the novice, to practice focusing with this directly 
produced image using wires or sharp edges placed at the position of the object. This 
image can be examined in a variety of ways. It can be positioned directly on the face of 
a TV camera and examined on a CRT monitor. This is, perhaps, the most ideal method for 
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viewing, since no safety considerations need be taken. The Image can be viewed on the 
frosted glass of a macro-camera simply by focusing the camera on the actual spatial image 
which appears behind the hologram position. Finally, for the experienced observer, the 
simplest method is to view the spatial image directly with a microscope or an eyepiece 
and no frosted glass. In this manner, the viewer can observe the actual image free of 
speckle, which is introduced when a diffuse surface is used. The image can also be 
examined in this mode with a macro-camera, if a clear focusing screen is used in the 
camera as opposed to the usual frosted screen. 
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EXAMPLE: When viewing these images directly, safety is a primary consideration. The 
maximum acceptable emission limit established by the American National Standards 
Institute shows that, under almost all circumstances, this direct Image may be safely 
viewed with the proper CW laser. Acceptable Emissions, taken from Table 4 of the above 
guide, state that for the helium-neon wavelength for emission durations between 10 
seconds and 10,000 seconds, 3.9 x 10” joules of energy are permissible. A 4 milliwatt 
alignment laser, used in the holocamera described above, allows approximately 1,000 
seconds viewing the direct laser beam at the hologram position before exceeding the 
maximum exposure limit. In the holocamera system, the beam reaching the hologram has 
been expanded from 2 mm to approximately 100 ran, reducing the laser intensity by a factor 
of approximately 2500. ____ 

To develop skill in focusing an image, one must understand a few characteristics of 
coherent light. Consider Figure 3.11. A sharp edge is imaged by lens into a space where 
it is examined in fine detail bv a viewer. It is always best to search for a sharp edge 
or a small object to determine focusing effects most critically. The object is imaged as 
a shadow and in the exact image plane the shadow has very sharp edges with extremely fine 
diffraction effects. 
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Figure 3.11 Focusing on a coherent image of an edge. 

In the actual plane of the object, light is scattered by the edge both into the 
shadow region and into the region where light passes unscattered by the object. The 
scattered light and how it varies as the image goes in and out of focus can be observed 
with a magnifier in the image to determine with high precision the actual focus position 
of the image. In the illuminated region, this scattered light mixes with unscattered 
light, forming a diffraction pattern which at first is extremely fine. As one proceeds 
further and further out of focus, the diffraction lines become broad and more 
distinguishable. 

One difficulty in focusing occurs when one is extremely out of focus. When the 
observed plane is badly out of focus, a motion toward the focus plane will not produce a 
drastic change in the diffraction pattern. Only when one nears the focus plane do the 
collapsing diffraction patterns and the illumination scattered into the shadow become 
clearly observable with rapid change such that an accurate focusing can be determined. 
Therefore, one must first find approximately the focus plane and then examine the micro¬ 
scopic details of the image edges. 
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A common mistake made by the novice la to expect that since parallel light la 
passing over m object, there la no focus plane. This would be true only if there was no 
diffraction. Diffraction allows one to locate, very precisely, the focus plane of the 
image. Similar focusing characteristics can be observed for wires, straight edges, and 
particles, and it la useful for an observer to practice focusing on such targets. 

Focusing is best done on an object that is strongly diffracting, such as an edge. 
Phase objects also diffract light, but except for step changes in refractive index, they 
are more difficult to focus because they are not strongly diffracting, and the effects 
described In the previous paragraph are not as easily observed. Therefore, when the 
field of interest contains no opaque objects, it is often difficult to locate exactly 
where the focused image is. This Is another reAson why it is Important to locate the 
approximate position of the Image, relative to the hologram, before the hologram is 
actually made. 

3.1.6.4 Optimizing for Data Reduction - A principal advantage for holography is that the 
entire wavefront passing through tKe interest field can be recorded and reproduced 
exactly. In practice, however, real effects make it important to choose a configuration 
that has been optimized for the highest quality and easiest image analysis. Holograms 
can be made without lenses; however, it is also highly advantageous to use the best 
properties of holograms and lenses In combination to get the best results. Even so, 
every additional optical element added to a system increases optical noise. For this 
reason, one should limit to a minimum the number of lenses, windows, and mirrors required 
to produce optimal results. 

The accuracy with which reconstruction angles are to be matched depends ultimately 
upon the desired resolution. Image resolution in the range from 5 to 10 microns does not 
require a highly critical alignment of the hologram. Image resolution below 5 microns 
requires that this angle be aligned very critically. This is commonly done by making the 
hologram of a resolution chart and adjusting the angular position of the hologram until 
the resolution is optimized in the reconstructed image. 

It is usually helpful to position the image at a place that is not located in the 
overlapping reference wave, making it simpler to view without interference by the 
reference wave. It is also desirable to place this image at a location where it is 
easily accessible to a recording camera. At this stage, it may even be desirable to 
produce an image size that is the same size as the format of the camera. Again, this 
must be traded off with resolution. It is easiest to obtain high resolution by magni¬ 
fying before holographic recording, thus making the overall reconstructed image size 
quite large. For this reason, large format cameras are desirable for reconstruction. 

Further magnification of the reconstructed image is often useful, and is limited 
only by the quality of the hologram itself. The simplest high resolution image that can 
be obtained is the in-focus image in the plane of the hologram (known as image plane 
holography). This image, while actually reconstructed as a three-dimensional image, is 
(in the exact plane of the hologram) also a photograph and, therefore, is not as 
seriously affected by the holographic parameters as are images out of the film plane. 
This image can, be observed with incoherent light, and suffers least the type of aberra¬ 
tions such as chromatic, spherical, and astigmatism, commonly associated with holographic 
reconstructions. As a general rule, the further the distance from the hologram to the 
image plane, the more critical the hologram parameters (such as alignment) become. 

3.1.6.5 Lens Configurations - In an imaging holcamera, such as the system described 
above, the possibilities For positioning and magnifying the image before and after 
reconstruction are many, and it is important to understand how these affect the ultimate 
appearance of the data. It is extremely important to know in advance the probable 
location of the image. This can save a vast amount of confusion and delay during the 
reconstruction. It also allows one to optimize the reconstructed data for whatever 
recording or viewing system is used. Three classes of imaging systems must be 
considered. 

(1) A collimator pair which images one-to-one in three dimensions. 

(2) A collimator pair which magnifies such that the magnification in the image field is 

uniform everywhere, but the magnification is not the same in the depth dimension as 

it is in the planar dimension. The lateral magnification is the ratio of the focal 

lengths of the two collimators used. The longitudinal magnification is the square 
of this ratio. 

(3) Use of the single lens to image and magnify before reconstruction. 

In the following discussion, each of these lens systems and its imaging properties 
is described in detail. For the first part of this discussion, the hologram itself is not 
mentioned. The reader must only remember that the hologram itself only recreates what 
exists optically, so, in any of these discussions, the hologram position could be at any 
place in the entire diagram. Ideally, with or without the hologram, the image space 
stays the same. One can consider the image space as being produced by the alignment 
laser and the holocamera lenses when the hologram is not present, or being produced by 
the hologram when it is set up in a reconstruction system. From an optics point of view, 
when the optical wave from the hologram is examined, one cannot discern whether it has 
come from the hologram or from the imaging system in the holocamera, since the two should 
be identical. 
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Figure 3.12 describes the imaging properties of a collimator pair which images one- 
to-one in three dimensions. This system allows the largest field of view and, because 
the image space is identical to the object space, is the simplest system to analyse. 

There are three separate cases that must be considered for use with a holocamera 
when the two matched lens collimator pair is used. Case 1 , which is illustrated in 
Figure 3.12a, has'a sample space defined by a volume that is two focal lengths long and 
one lens diameter wide. The hologram can be positioned anywhere on the right-hand side 
of the second lens. In some cases where an Image plane hologram is desirable, the image 
of the sample volume exists on both sides of the hologram, as shown in Figure 3.12b. As 
previously mentioned, it is sometimes convenient to position the hologram such chat all 
of the sample space appears on only one side of the hologram. When the image exists on 
both sides of the hologram part of it is a real image and part of it is a virtual 
image. The real image can actually be presented directly to a frosted screen for viewing 
or directly to a lensless camera. The virtual part of the image must be relmaged into a 
real image before viewing on a frosted screen or recording. 
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Object Space is represented by A, B, C, D, E, F, G. The corresponding image 
(congugate) points are primed (A’, B*. etc.) respectively. Note that the relative 
positions of the object and image space are fixed regardless of the location of 
the lens pair! In lens Position 3 note that the object and image spaces overlap. 

Figure 3.12a One-to-one imaging In three dimensions with a matched collimator pair. 

The part of the image that is virtual can be made real in holography simply by 
reversing the direction of the reference wave and projecting backwards through the 
hologram as shown in the third figure of Figure 3.12b. In this case the reference wave 
comes from the right and the viewer sees the reconstructed image on the left-hand side of 
the hologram. The part of the image that was virtual now becomes real. This process is 
often useful when presenting an image directly to the fac<» of the vidicon, without using 
other lenses, where part of the field of interest happens to exist behind the hologram in 
a virtual space. This procedure allows one, by turning the hologram over, to convert 
this virtual space into a real space. 

The two lens image methods just described applies both lenses before the recording 
of the hologram. In an equivalent technique, one can apply one of these lenses before a 
holographic recording and the other during holographic reconstruction. This is illus¬ 
trated in Figure 3.12c. The recording is made with the hologram in a position just 
forward of where the second collimator would have been placed. The diverging wave is 
then recorded in the hologram. During reconstruction, the reconstructed wave will be 
that diverging wave. The collimator that was used during recording can then be used as 
an exact matched collimator to collect this diverging wave and recollimate it, producing 




a one-to-one image of the original sample volume. The points A', B', C' in this image 
correspond to the original A, B, and C in the sample volume. 




RECONSTRUCTING THE ORIGINAL WAVE 
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RECONSTRUCTING A CONJUGATE OF THE ORIGINAL WAVE 

Figure 3.12b Holography with a one-to-one imaging system. 

To remove any imaging defects (see Figure 3.12c), the conjugate wave traveling 
antiparallel to the original wave is projected backwards through the lens used during 
recording, and the image appears once again in the original sample space. 

Figure 3.13 illustrates the process of imaging with magnification. Two separate 
cases need to be considered--one in which the image space is a constant magnification of 
the object space, and the other in which the magnification of the object space depends 
upon the location in that space. The top diagram illustrates the imaging for the former 
case. Here, two collimators with different focal lengths are used, taking a collimated 
wave from the object space and reproducing a collimated wave in the image space. The 
lateral magnification (x,y) is the ratio of the two focal lengths. Magnification in the 
z dimension can be seen from this diagram to be the square of the magnification of the x 
and y dimensions. For this lens, therefore, magnification in the x and y dimensions is 
3. while magnification in the z dimension is 9. For example, the distance BC is Imaged 
B , C', nine times the original object dimension. 

This type of lens configuration is extremely useful because the Image space is 
magnified by a constant value independent of the position in the original object space. 
The volume lying from one focal length in front of the short lens up to its surface is 
Imaged to the right-hand side of the large lens in a space is convenient for producing 
holograms. As with the previous lens, the hologram can be positioned at any place in the 
image space. Precisely the same method for handling the reconstructed image for the 
matched collimator pair can be used here. If any part of the image is a virtual image in 
the reconstruction, it can be converted to a real image simply by reversing the direction 
of the reconstruction reference wave, a practice which is achieved normally by turning 
the hologram over. The stretching out of the image space in this reconstruction does 
make the analysis of the reconstructed image require more space, if the object Itself has 
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considerable z dimension. When such is the case, it is often better to reconstruct the 
wave back through the lens pair to create an image space identical to the original object 
space. 



CASE 2: RECORDING WITH ONE COLLIMATOR 



RECONSTRUCTING ONE TO ONE WITH THE SECOND COLLIMATOR BEING SUPPLIED DURING 
RECONSTRUCTION. 



Figure 3.12c Using one lens during recording and the same lens during reconstruction 
to achieve one-to-one imaging. 

When an object field is magnified with a single lens, the amount of magnification 
depends upon the position of the object. The closer the object is to the focal point of 
the lens, the greater the magnification. When a very thin object field is used, little 
problem is created, since the amount of magnification is easy to determine simply by 
knowing where the object or image field is. The equations are provided in Figure 3.13. 
The magnification is simply the ratio of the image to object distance; or, when either 
the image or the object distance alone is known, one of the other two relationships for 
magnification can be used. For example, if only the Image distance is known, as might be 
the case in observing a reconstructed image, the magnification is provided by the 
difference between the image distance and the focal length divided by the focal length. 

When a larger sample volume is desired with magnification, it is sometimes 
convenient to produce a constant magnification image space. This can be done either 
during recording, as Illustrated in the top figure of Figure 3.13, or during recon¬ 
struction as Illustrated in Figure 3.14. 

By studying Figures 3.12 through 3.14, it is easy to see how, by using different 
lens configurations, the image space of the hologram can be placed almost anywhere. 
Without some predetermination of this image space location, it is likely to wind up in an 
inconvenient space and require reimaging. Locating the image space can be done by using 
the relationships shown in these figures, or as described previously, by inserting the 
calibration object in the object space and observing the location of its image relative 
to the hologram. 
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Lens 2 




IMAGING WITH A SINGLE LENS PRODUCING MAGNIFICATION WHICH VARIES WITH AXIAL POSITION. 


Figure 3.13 Imaging with magnifications. 

3.1.7 Holographic Interferometry 

Interferometry is the mixing of two mutually coherent waves to produce interference 
fringes that can be used to determine the phase distribution of one of the waves relative 
to the other. In holographic interferometry (HI), one or both of the waves Is derived 
from a hologram. 

3.1.7.1 Techniques - Basically, there are about seven different types of HI (see 
Figures 3.15 and 3.16) . These include: (3.15A) interferometry with one nonreconstructed 
wavefront serving as the interferometry reference wave, (3.15B) real time holographic 
interferometry, (3.15D) '’common path" interferometry (e.g., wave shearing), (3.16C) 
interferometry with two time separated wavefronts from a single plate (e.g., double 
exposure or double pulsed), (3.16A) Interferometry of wavefronts from two holograms 
(e.g., double plate or sandwich), (3.16B) interferometry o c wavefronts reconstructed with 
two reference waves. Heterodyne HI applies two reconstructing reference waves of 
different wavelength, producing moving fringes. Phase shift interferometry phase shifts 
the reconstructed waves in steps by independently phase shifting the reference waves. 

3.1.7.2 Equations - In the following, we analyze the most general case, where the 
interfering waves are derived from two different holograms made at different times, 
separated by az , moved laterally and tilted with respect to each other between 
exposures, and made with two different reference waves. The resulting general equation 
can then be specialized to describe all of the cases illustrated. 
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CASE 3: RECORDING WITH ONE MAGNIFIER 



Image of 
U Flglrtai 
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RECONSTRUCTION BACK THROUGH THE ORIGINAL PATH 

Figure 3.14 Holography with Premagnification. 

The general equation will describe double plate HI. By setting az equal to zero, 
double exposure, double reference wave HI results. By setting the two reference waves 
equal, double exposure HI results. By allowing one of the reconstructed waves to emulate 
a real time wave, real time HI results. 

In this computation, we are Interested in the effect of the relative position of 
the holograms on the interferogram. We assume that a plane reference wave is used with 
an originally plane object wave that has been modulated by a phase factor eU c ( x ), and 
recorded first on a Plate 1 , then on a Plate 2 which is located behind Plate 1 
(Figure 3.17). Then the two holograms are illuminated by a plane reference wave 
simultaneously reconstructing object waves from both plates. Let us fix the z position 
of Hologram 2 and assume that the position of Hologram 1 can move after the recording, 
and compute a one-dimensional (x) hologram. 

We ultimately require the general equation of interference of the two waves 
UqJ and UqJ as a function of az . The following are the object and reference wave 
amplitudes*in the front plane of the respective holograms. 



U 02 - u 02 el ^ (X ’ ln t02) + 1>02 + IS(X ' “ Can * 02) + l *° 2(X> 
U r2 - u r2 e l I* (* « ln *r2) + 1 *r2 + l8 (* * tan *r2) + l *r/ X 


(3.4) 


(3.5) 
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D. TWO REFERENCE WAVE SHEARING 

Figure 3.15 Holographic interferometry with singly exposed holograms. 


U 0 , - u 0 , e 1 f* (* sln *01 - az C08 *01) + *02 (3.6) 

U r , - u r , e 1 7* ( x 8ln *rl - 4Z 008 *rl) + *r2 (3.7) 


where ♦ . and a . are phase modulations of each wave. contains the flow fi 
information plus any aberrations due to poor optics, while > ^ describes the affect 
poor optics on the reference wave, and 8(x) is the aberration introduced bv Hologram 1 
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make the hologram, object and reference waves are added, giving 
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(3.8) 

(3.9) 


Now we assume 
exposure. 


hologram which has an amplitude transmissivity which is linear in 


hi 


(3.10) 


During reconstruction, the reference wave is returned to the hologram and is modulated by 
Thf 
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The second term contains the reconstructed wave that is of Interest. Namely 
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(3.11) 
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C. DOUBLE PUL8ED HOLOGRAPHIC INTERFEROMETRY 


Figure 3-16 Holographic interferometry with multiple recordings. 

Now we examine the mixing of the two waves at the second hologram surface. We 
assume that the wave from the first hologram propagates with negligible diffraction over 
the distance az . The amplitude of the reconstructed wave from the first hologram at the 
second hologram surface is (Figure 3.17) (assuming the same wavelength is used) 

U 012 ’ “of e 1 I* cos *01 + e ( x ‘ **' tan *01) (3.13) 

We now mix UqjJ and U^' and square to get the interferogram Intensity. 
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Figure 3.17 Double Plate Holographic Interferometry. 
Substituting 3,12 and 3.13 into 3.14 leads to 


I, - I, + 2Re Ur U_, U„„ U. 


ri U r1 U 01 u r2 U 02 


tz' cos * 01 


Xe -l8 (x - A*' tan * 01 ) e t(*' 1 - i i2 ) 


(3.15) 


Get U Q £ and U r ^ from Equation 3.4 - 3.7. 

To include the effect of moving Hologram 1 In the x direction during reconstruc¬ 
tion, note that the hologram transmissivity remains the same except for an x translation 
of ax. Namely 

T hi - U r ; (x + AX) U 01 (x + ax) (3.16) 

and this must be used in the above equation. And 

Ij - I g + 2 u rf u r , (x + ax) u 01 (x + ax) u r J u r2 u Q2 (3.17) 

times 

2 

cos y 1 (x sin - az " cos * r1 - (x + ax) sin * r1 ) 

+ AZ COS * r l 

+ (x + ax) sin * 01 - AZ cos * 01 
- - sin * r2 ) 

- *(sln * 02 + az- cos * 01 ) + » 02 - 8 (x - az tan * Q2 ) 

+ *01 (x + 1x5 ■ *02 (x) + *r1 (x) " * r 2 (x) + *r1 * *r2 
+ 8 (x - az' tan * 01 ) 
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Ij - Ig + I g cos [ Y 1 (* (sin * r f - sin , r , + sin * f2 - sin * r2 (3.18) 

+ sin * Q) - sin * 02 ) + 4X (sin - sin * r ,) 

+ (4z - 4Z-) (cos * rI - COS S 0 ,)) 

+ 1*0,(x + 4x) - * 02 (x)) + i* rI (x) - S',(x)] + [* r2 - s' 2 ] 

+ S (x - 4Z' tan So)) - B (x - 4z tan . Q2 ) ] 

Note that we have still retained two different reconstruction reference waves so that the 
equation is general. 

Now consider some special cases. 


Case 1 

(a) Single plate (az ■ ax - 0) 

(b) Same recording reference waves (# r1 ■ an< * 

(c) Single reconstruction wave and *^1 " *' 2 ^ 

(d) Same original object wave (*qj * 

The resulting interferogram is described by 

I x - I B + I s cos [AA 0 (x)] (3.19) 

Note that all of the aberrating effects have cancelled out of the equation. Also note 
that reference fringes (finite fringe interferogram) could have been added by changing 
conditions between the two exposures, so that * j * ♦ 2 or *01 * ^02 or 

Equation 3.19 describes double pulsed holographic inxerferfcroetry, and shows how it 
depicts the phase change a$q(x) between the two exposures. 

Case 2 - Double Reference Wave Holographic Interferometry (DRWHI) 

Assume: 

(a) Single plate (az - ax - 0) 

(b) Same original object wave (*Q 1 - * 02 ) 

The interferogram is described by 

l I ~ h + l s 008 fr 1 ( & 8ln *rl + 4sln *r2 + sln *01 - sl "*02) (3>20) 

+ 4* 0 (x) + . r1 (x) - .',(x) + 4 r2 - t' 2 ] 

Equation 3.20 describes an extremely powerful technique. An interferogram 

containing the data, which is a$q(x) , is produced, and it contains a completely 
controllable set of reference fringes. The controls include: 

(1) Spacing adjustment of the reference fringes during reconstruction by adjusting the 
relative angles of the two reference waves. 

(2) Orientation of the reference fringes. 

(3) Phase shift of the reference fringes by adjusting relative phase of the reference 
waves. 

(A) Fringe contrast adjustment during reconstruction by relative intensity adjustment 
of the reference waves. 

Note that the phase errors in the object and reference waves carry through directly 
into the final interferogram. However, if the same reference waves are used for 
recording and reconstruct ion, the error cancels. 

Case 3 - Double Plate Holographic Interferometry (DPHI) 

Assume: 

(a) The special case where identical reference waves are used for recording and 

reconstruction 


♦rl - M ■ *t2 ‘ *r2 " *r and *r1 " *r2 " M " *r2 





(b) The holograa norual bisects the angle between the object and reference wave 


(♦r “ '*01 " '*02) 

The resulting Interferograa la described by 

Ij " Ig + J s eoa ( 2ax ^® in *01^ + 4 * 0 ( x > + 48 ) (3.21) 

where at Is Che difference In the phase error Introduced on the object bean by passing 
through the first hologran fron that Introduced onto the reconstruction reference wave 
for Che second. If the hologran does not change Its thickness between recording and 
reconstruction, this can be zero. Also note that the reference fringe phase can be 
shifted by a relative lateral adjustnent between the two holograms during reconstruction. 

Case 4 - Double Holographic Interferonetry with Tilt In One Plate 

(a) Consider a slight tilting of one of the holograms during reconstruction, say 
hologran Number 2. 

(b) Assume az ■ az' and a* “ 0 

and for Identical reference waves during recording and reconstruction 

(♦rl " *r2 _ and *r1 " *r1 _ *r2 “ *il) 

Since the second plate la tilted during reconstruction 

♦r2 ' *r2 

Let 


*r2 • *r1 + 8 


(3.22) 


where s Is the tilt angle (See Figure 3.18). 



Figure 3.18 Reconstruction with a film tilt error. 
Kquatlon 3.18 becomes for small e 

2 

*1 “ *B * J S eos I T 1 I x ( 8 cos *r1 + * In *01 ‘ * ln *02 ) 
+ ( co ** r i ' C08 *oi ) 8 } ♦ 4 *o fx > ] 


Now consider a case where * 0 j - » Q2 and - -y 0) . Then 

2 

*1 ■ *B + *S 008 i IT x ( 8 co8 *r1 ) + 4 ♦o (x) 1 


(3.23) 


(3.24) 


Note that the phase shift across the lnterferogram Is minimized with tilt when the object 
and reference waves are symmetrical about the hologram normal during recording (l.e., 
when * r , - -eg, ). 
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Case 5 - Real Time Holographic Interferometry (RTHI) 

The reconstructed wave, made at time zero, Is mixed with a real time wave that has 
passed through the region of Interest, providing an interferogram that depicts the phase 
difference of the real time wave and the time zero wave. Equation 3.18 can be applied 
directly with some analogy. 

This Is normally done by viewing the illuminated real time scene through the 
hologram. Assuming that the hologram is not moved after Its original recording. 
Equation 3.18 reduces to 

Ij - r B + I 8 co * if* ( xs In +' - sin# r + - sln*Q 2 ) + A*(x)] (3.25) 

Note that If the same reference wave Is used, all aberrations cancel. Even the 
aberrations introduced on the waves by the hologram substrate are removed, since they are 
the same for both interfering waves. This important result allows high quality inter¬ 
ferometry to be done through relatively low quality photographic film and glass plates 
without added error. 

Case 6 - Single Exposure. Live Reference Wave HI 

A degenerative case of real time HI is that in which an interferometer reference 
wave derived from a laser is mixed with a reconstructed data wave. If the reference wave 
is passed through the hologram, then Equation 3.25 applies. This method, which has been 
used to cancel otherwise troublesome film aberrations, has been termed the FACT (Film 
Aberration Correction Technique) method. The advantage is that single exposure holograms 
can be used for holographic interferometry with considerable fringe control available 
during data reduction. 

Case 7 - Local Reference Wave HI (LRWHI) 

Transferring the reference wave around the field of interest can sometimes present 
impossible requirements, for example, in large wind tunnels, long range cases, and 
holography between two moving platforms (such as aircraft). For such cases, the 
reference wave can sometimes be transmitted through a part of the field of interest that 
does not add phase information to it, for example, through a laminar part of the flow 
field. Indeed, even if phase information is added to the reference wave, it can be fil¬ 
tered from the wave before it is mixed with the object wave, as described for local 
reference wave interferometry in Chapter 2. 

In LWRHI, the object wave is split into two or more waves after having traversed 
the field of interest, one serving as the information wave, the others serving as the 
reference waves. It is usually best, at this point, to remove from the latter that part 
of the wave that has traversed the strong part of the flow field by aperturing, if 
possible. After this operation the wave may be a satisfactory reference wave. This can 
be determined by looking at the resulting interferogram which characterizes the region 
from which the reference wave was extracted. 

If there is no weak part of the flow field, then it is necessary to spatial filter 
the phase information from the wave by passing it through a pinhole spatial filter to 
make it a useful reference wave. Just how much error this adds to the final 
determination of phase can be seen by examining Equation 3.18, which contains a term that 
represents the aberration of the reference wave. 

Case 8 - Glint Holographic Interferometry (GHI) 

A degenerative case of LRWHI occurs when some part of the object or object field 
contains a suitable point source (glint) that can be used as an object wave (a mirror can 
actually be attached to the object). Consider the holography of the surface of a fast 
moving object, such as a projectile that has a shiny part. Light from the region of the 
shiny part may serve as a reference wave. 

Deriving the reference wave from the object wave may provide several benefits. 
Producing a hologram of the surface of the projectile may, under normal conditions, 
require extremely short exposures to insure that the holographic fringes remain station¬ 
ary during the recording. However, if the reference wave is derived from the object 
itself, both the object and reference waves are shifted equally in phase by the motion of 
the object, greatly relaxing the recording conditions on stability. 

Of perhaps even greater Importance is the impression of the same phase modulation 
on the object and reference waves by the media between the object and the recording. 
This results in the cancellation of such information in the recording as though it were 
not there at all. This might have application in holography through a strong flow field, 
or in the holography of an object at long distance through a disturbing atmospheric 
condition. This principle has been used in the holography of burning droplets in rocket 
exhausts . 17 

3.2 Applications 

In the following, a selected number of applications are presented to illustrate the 
use and state of the art of holography in the aerospace field. Holography is now being 
used in hundreds of laboratories. It is not attempted, however, to provide details and 
credits for all of these here. Examples have been contributed by the laboratories 
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visited or otherwise contacted, and this does not necessarily mean that they represent 
the originators of the various techniques shown. 

3.2.1 Particle Field Holography 

The study of dynamic particle fields was the first real application of holography, 
and is still one of the most important scientific appl nations. Holograms can optically 
freeze a three-dimensional, high resolution image of a ^article field that can be studied 
microscopically long after the actual particle field ceases to exist. Multiple exposure 
holography further allows the study of the dynamics of the field. There is currently no 
other way to perform such a three-dimensional study. 

The field has been enriched by many new refinements and techniques since its first 
use in the early sixties. Included in the new technologies are the following: 

(1) Methods for aberrations corrections to provide diffraction limited imaging with low 
quality optics. 10 

(2) Methods for subtracting out optical noise normally present in holographic images. 19 

(3) Holographic movies. 20-22 

(A) Methods to automatically extract the data from holograms. 23-25 

(5) Multiple wavelength holography to encode exposures. 20 ' 27 

(6) Multiple reference wave holography to provide two Independent looks at a phase or 
particle field in the same hologram. 20 

(7) Thermoplastic recording for in-place electronic processing. 29 

Many types of particle fields have been studied using by holography including the 
following: 

(1) Ice, snow, and water droplets in the field of meteorology 30 

(2) Pollution particulates and fibers 31 

(3) Fiber and particle generation in, for example, the fiber insulation 
industry 32 

(A) Spray nozzle development and characterization 33 

(5) Particle formation and breakup 34 

(6) Explosions 35 

(7) Impact phenomena 30 

(8) Two phase flow 37 

(9) Combustion 30 

(10) Crystal formation and growth 39 

(11) Microbiology 40 

(12) Cloud study from aircraft 41 

(13) Flow diagnostics by particle tracking 42 ' 43 ' 44 

3.2.1.1 Spray Nozzle DiagnosticB - Figure 3.19 Illustrates data taken from a doubly 
exposed hologram of a spray nozzle. 23 Though the difficulty of presenting three- 
dimensional data in a two-dimensional format has not been overcome, this figure illus¬ 
trates one technique of presenting such data. An overall shot of the field is recorded 
at relatively low resolution to provide a global look at the subject field. Then 
microscopic higher resolution Inserts are added to show specific details. As previously 
computed, such a hologram contains many such photographs. The illustrated recording is 
from a 10 cm diameter hologram. So there are about AO,000 different 10 cm diameter 
photographs available (with 5 micrometer resolution) in this hologram and many more at 
lower resolution. 

3.2.1.2 Particle/Shock Wave Interaction - Figure 3.20 Illustrates the study of the 
rebounding of particles through a shock wave after impact with a hypervelocity 
projectile. 36 This is a classic case in which high resolution analysis of a particle is 
required at two different times in a dynamic situation and where the particle location in 
the volume is not known. A photographic solution would have required many cameras 
focused at different depths in the field of Interest. 

3.2.1.3 Holographic Subtraction - Figure 3.21 Illustrates the use of holographic 
subtraction to improve the signal-to-noise ratio in a reconstructed image. 45 In this 
example, the hologram is made with two reference waves, one for the first exposure and 
one for the second. The reconstructed images of each exposure are Illustrated in a 
and b. In c the two images are reconstructed simultaneously. The optical noise Is quite 
severe in all of these. If, however, during the reconstruction process, one of the 
reference waves is phase shifted by 180 degrees, then the backgrounds will subtract and 
the Images will reverse from dark to bright. The improvement In S/N is evident. 

3.2.1.A Data Reduction - Particle field holography is still severely limited because of 
the difficulty in extracting the data from the hologram. Some image analyzers exist that 
augment data extraction, but, as of yet, no commercial Instrument can handle completely 

















Addition: Reconstruction of t + 2 Subtraction: Reconatructlon of T - 2 

Figure 3.21 Using holographic subtraction In double exposure particle holography 
to improve slgnal-to-nolse ratio. 

the problem of distinguishing In and out of focus Images and the severe coherent noise 
commonly associated with such Images. 

The most basic data reduction procedure Is to scan the Image volume while counting 
and sizing Individual particles. Particle field images are usually observed on a closed 
circuit TV monitor or In photographs. In either case, It Is Important to understand that 
the depth of field Is a function of particle size, so that the sample volume must be 
corrected for particle size when particles are counted In a single plane. 

An alternative method Is to examine the scattered light field from the hologram. 
It can be shown that the Fourier transform of a scattered light field can yield the size 
distribution. Using this principle, Ewan 21 * (1984) and Hess et al.*3 have demonstrated 
that the size distributions of particle fields recorded in holograms can be extracted 
almost Instantly. 

The limitations up to now on holography data reduction, and the requirement to 
record at relatively slow rates, Is somewhat confining. Other particle sizing Instru¬ 
ments that can provide real time data at high rates excel for some requirements, while 
holography la considered to be the desired method to attain global looks and occasional 
full field examinations. 

3.2.2 Flow Visualization Holography 

When a light wave passes through a flow field, the optical pathlength is longer in 
regions of higher refractive Index, so that the parts of the wave that pass through such 
regions are phase shifted, and the emerging wavefront Is modulated with phase Information 
that characterizes the flow field. As shown In Chapter 1, the refractive Index Is 
directly related to density In gases so the gas density Is sometimes derivable from a 
measurement of phase of the vavefront. 

Shock waves are characterized by a discontinuity In refractive Index, and they 
strongly diffract light. Therefore, they can be focused In 3-D like any object. Smaller 
refractive Index gradients refract the light at smaller angles and cannot be observed In 
an ordinary Image of the field; their effect on the phase of a wavefront must be observed 
with other methods to be described. 







In fLov visualization, since the object is usually a phase object (for example, a 
distribution of refractive index), off-axis holography is almost always used. A 
wavefront that has been modulated by the flow field can be recorded in its entirety for 
later analysis. The later analysis may comprise any of the ..conventtonal methods, such as 
shadowgraph, schlieren, deflectometry, Hartmann, moire, or interferometry (see 
Chapter 2). Choice of the method can be postponed, and, indeed, all of the methods can 
be applied. Moreover, many variations for each method can be applied. The seven types 
of interferometry can be applied, the shadowgraph focus can be varied, the schlieren 
knife edge can be adjusted, and the deflection sensitivity can be varied. 

A unique advantage of holography is that a wavefront modulated by the optical 
system only can be subtracted from one that has both optical system and flow field 
modulation, removing the effects of the optics. This can be extremely valuable in large 
systems, where cost of large diffraction limited windows and mirrors can be prohibi¬ 
tive. This principle can be carried further to subtract effects of a steady part of an 
object field, leaving an unsteady part which might otherwise be overshadowed by the 
steady part. Figure 3.22a illustrates a case in which a weak unsteady flow is 
overshadowed by a strong steady flow and can be seen clearly only after subtracting off 
the steady flow effects. Note the vortex street in the model wake, the unsteady shock 
wave, and the free stream turbulence in Figure 3.22b. One can zoom in on specific parts 
of the flow fieLd with very high resolution as illustrated in Figure 3.23. while other 
parts of the flow field do not require such resolution for analysis. 



Figure 3.22a Potential flow field. 





Figure 3.22b Potential field holographically subtracted leaving unsteady field. 


3.2.2.1 State of the Art System - Figure 3.24 illustrates a conceptual YAG laser based, 
state of the art holographic flow visualization system. YAG lasers have become popular 
in this application because of their repetition rate, which is higher than ruby lasers, 
the first used in this application. Other pulsed laser candidates are dye lasers and 
Alexandrite, which all have pulse durations in the twenty nanosecond range. CW lasers 
can be used if the entire system Is isolated from vibration, or if the laser can be 
properly shuttered to sufficiently short exposures. In general, the system would include 
two individually modulatable reference waves, a thermoplastic recording device, and an 
on-line image processor to analyze the interferograms. The methods illustrated in Figure 
3.15 and 3.16 would then be applicable. 

A good example of this type of system has been developed by Craig et al. 46 and is 
in use in the NASA Ames Research Center Two Foot Transonic Wind Tunnel. The system 
employs a frequency doubled YAG laser that produces 20 millijoules of energy per pulse. 
Craig has recently demonstrated the on-line capability of the system, by employing a 
thermoplastic device that can reconstruct data directly into an automated fringe 
reduction system developed by Tan. 4 ® 

This basic type of system has been in general use in many wind tunnels since the 
early 1970s 13 * 49 ' 63 , but only in recent years has hardware and technology evolved to the 
point of on-line data capability. The applications described in the following sections 
are largely variations of the system shown in Figure 3.24. 

3.2.2.2 Shock Wave/Boundary Laver Interaction - The interaction of shock waves with 
boundary layers is an extremely important process in aerodynamics. This process is also 
too complex to completely describe analytically. Therefore, experimental aerodynamics 
roust provide quantitative results for the development and evaluation of theoretical 
models and computer codes. Holographic flow visualization has played an Important role 
in providing such results. 
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Figure 3.23 Wake region illustrating high spatial resolution. 
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Figure 3.24 Conceptual state-of-the-art holographic flow visualization system. 
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Figure 3.25, from the work of Havener 61 ”®3 in the study of shock wave/boundary 
layer interaction, illustrates various capabilities offered by holographic flow visuali¬ 
sation. These also provide a good comparison of three different flow visualization 
techniques; shadowgraph (a), infinite fringe interferometry (b), and finite fringe 
Interferometry (c). The interferograms were produced by the double plate methods, one 
hologram being of the empty wind tunnel and the other being of the Mach 3 flow field. 
Figure 3.25d provides a more detailed view of a shock wave/boundary layer Interaction 
with an interpretation by Havener. The axlsymmetrlc parts and the two-dimensional parts 
of the flow field can be analyzed quantitatively from these interferograms, but the 
three-dimensional parts cannot. 



Figure 3.25a Shadowgraph - Flow visualization holography of a Mach 3 flow. 

One hologram provides an array of flow visualization types. 



Figure 3.25b Infinite Fringe Interferometry - Flow visualization holography 
of a Mach 3 flow. One hologram provides an array of flow 
visualization types. 
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Figure 3.25c Finite Fringe Interferometry - Flow visualization holography 
of a Mach 3 flow. One hologram provides an array of flow 
visualization types. 



Figure 3.25d Separation of a turbulent boundary layer from a flat plate 
due to a shock wave impinging on the boundary layer. 
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work of Rolls-Royce 77 ' 78 , using double pulsed hologram interferometry of the flow between 
compressor blades in a turbomachine. The flow is made visible because the position 
changes between exposures, as shown In Figure 3.28. This type of holography is used to 
locate the shock wave position in 3-D. It cannot be used to quantify the density. 
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a. Conditions at the edge of the boundary layer 
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b. Measurement in the Inviscid flow 
Figure 3.27 Diagnostics of the transonic flow in a channel. 



Figure 3.28a Method by which the holographic system was applied to the first 
stage compressor fan during a test. 
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movement' of 
blades during 
2 msecs 


shod* movement 



Movement of a normal shock between two pulses with a 2msecs 
separation 

Figure 3.28b Movement of the compressor blade fan assembly between two 
holographic pulses. 



Figure 3.28c Photographic reconstruction showing the shock structure at the 
tip of the compressor fan. 

Oftentimes, the experiment must be simplified to fit the capability of the measure¬ 
ment to obtain meaningful quantitative data. The flow cascade has been constructed in a 
two-dimensional form for this purpose. Figure 3.29 illustrates a double exposure HI of 
such a case, included with a computation of the flow field for comparison 611 . 

Another such example is shown in Figure 3.30. 7 * This is the flow field Inside an 
aerodynamic seal. This illustrates the appearance of vortices in an interferograa. 

3.2.2.4 Three-Dimensional Flow - As described in Chapter 2, tomographic solutions now 
exist for solving 3-b flows I? enough interferograms are made at the correct viewing 
angles through the flow. Only a few such systems are currently in operation; however, 
the number is likely to increase in the near future. A good example of such a system is 
in operation in the U.S. Army Aeromechanics Rotor Laboratory at Amea Research 
Center. 80 * 












Figure 3.29 Holographic Interferometry of transonic flow In cascades. 

The system Is a 60 cm diameter system of the type Illustrated In Figure 3.24, using a 
pulsed ruby laser. The object beam passes over the tip of a rotor blade, recording the 
flow In che region around the tip, (Figure 3.31). 

By traversing the syateu laterally, the view angle through the flow is varied 
automatically. The system must be fired by a signal that Indicates that the blade Is In 
the proper position. Holographic interferograns are typically recorded every five to ten 
degrees. For the views that contain large density gradients, store recordings are 
needed. Figure 3.32 *• is a sample of the data taken in this facility. Note the 
appearance of the shock wave at different angles of view. In some of the angles of view, 
the shock wave edge Is not visible. 

This type of system produces so much data that some type of automation in data 
reduction is Imperative. Considerable efforts to this end have been successful in 
producing systems that can digitize and analyze data from an interferogram In a few 
seconds. These will be described in a later section. 

The reconstructed density distribution (p/p 0 ), produced by the tomographic 
reconstruction code applied to 40 projections Is shown In Figures 3.32c-d for the heights 
z • 0.5, 1, 1.5 and 2 inches (chord C - 3 Inches; the leading edge of the blade 
Y/C - -0.5; aspect ratio - 13.7). The optimum number of Iterations was 3 to 6 for the 
experimental data. The computation time on an IBM PC was approximately 3 minutes per 
Iteration. 

3.2.2.5 Other System Types - Some applications require special system design. Examples 
Include the requirement to Increase sensitivity, to reflect light from the model, or to 
look at very large flow fields. Figure 3.33 Is a system that reflects light from the 
model or a panel behind it. If the surface can be treated with a retroreflective 
coating, then very large areas can be recorded. Recently, holographic optical elements 
have been devised for the treatment of surfaces for flow visualization. 
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Figure 3.30 Flow field In a pneumatic bearing. 

Figure 3.34 la e ayetea that la dealgned to achieve improved eeneltlvlty. The 
object bean paaaea through the flow field nore than once, adding modulation to the wave 
each time It paaaea through the object. The object wave with the correct number of 
paaaea la aelected holographically by pathaatchlng. This unique coherence filter forma a 
hologram of only one of the many object wevee that etrlkee It. The othera aImply tend to 
fog the film with light that la not mutually coherent with the reference wave. To 
compenaate for the fogging, a critical part of the technique la to bleach the hologram 
before reconatructlng the data. Figure 3.35 la a act of data taken with the ayatem to 
llluatrate thla point. The object la e heated ball with natural convection around it. 
Aa the number of paaaea through the field doublet, then ao doee the fringe denelty In the 
lnterferograa for eaaentlally the same flow field. 
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Figure 3.32b Interferogram of flow over a rotating blade, e * 176,5°. 

3.2.2.6 Data Reduction - Figure 3.3* ’'epresents the flow of data in the holographic data 
reduction process. The first two steps, producing and reconstructing the holographic 
image, have been described above. The remaining steps are described here. 

For years, publications reporting uses of holographic interferometry focused on 
feasibility demonstrations, and the actual data output of such research was quite 
limited. More recently, attention has been directed to the complete extraction of data 
from the holograms. 82 •83 Furthermore, recent laser and recording technology makes 
possible the recording of thousands or even tens of thousands of holograms, placing an 
even greater demand on automated data reduction. 

The basic information contained in an interferogram is the distribution of the 
relative phase of two waves. This information can be extracted in essentially three 
different ways. 

(1) Analysis of a fixed interference fringe set. 

(a) Locate the coordinates of the center of interference fringes. 

(b) Identify a reference point from which to start. 

(c) Identify and number fringes integrally starting with zero at the reference 
point and adding (subtracting) where the phase is known to have Increased 
(decreased). Each succeeding fringe represents an increase or decrease in 
the phase difference of the two waves by 2 w. 

(d) The dependence of intensity in the fringe pattern is sinusoidal and it is 
most practical to locate maxima or minima, although, in principle, any equal 
intensity locus between any two maxima represents an isophase difference 
contour. 

(2) Analysis of a moving interference fringe set. (Heterodyne Interferometry) 

(a) The fringes are swept at a known rate (by methods already described), while 
monitoring the intensities at discrete points over the Interference region. 

(b) Locate a reference point and arbitrarily define the phase as zero. As the 
fringes sweep over this point, a detector will produce a time varying 
sinusoidal signal. 


















Figure 3.34 Multiple pees holographic flow visualization. 



(a) SINGLE PASS 


0>) DOUBLE PASS 


(c) FOUR PASS 


Figure 3.35 Multiple pass holographic Interferometry. 
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Figure 3.36 Data flow In holography. 

(c) Moving out from this point with a second detector, a sinusoid of the same 
frequency, but different phase, will be observed. The phase difference In 
the two time varying signals Is the desired phase difference between the two 
wavefronts. Each time this phase difference increases through 2 », the 
signal appears Identical to that observed at the reference point. 

(d) Numbering fringes becomes less <fi a problem, since the sign of phase shift 
la observable. 

(e) Interpolating between fringes Is no longer a problem, since they sweep over 

the detector. The resolution of the measurement is determined by the ratio 

of the detector size to fringe spacing. 

(f) The entire process can be achieved with one fixed and one moving detector 
or, alternately, a fixed detector array. 

(3) Analysis of the Fringe Set by Introducing known phase shifts while observing 

Intensity changes at discrete points. (Phase Shift Interferometry) 

(a) A detector is positioned at a reference point where the phase is arbitrarily 
set at zero. 

(b) Moving out from the reference point in steps of &y less than 2 • the 

Intensity la given by Equation 3.18. 

(c) The unknowns are l c , I g , and sy. By introducing two known phase shifts, 

three equations are oroauced. 

(d> The three equations can then be solved for ay in terms of Ij. 

(e) A moving detector or detector array can scan the interferogram. 

Unfortunately, the intensity in a real lnterferogram is much more complicated than 
that expressed in Equation 3.18, bebause of optical noise and parasitic Interference. A 
real lnterferogram Is characterized by: 

(t) Broken, discontinuous, and split fringes 

(2) Varying contrast 

(3) Extraneous fringes 

(4) Diffraction noise 

(5) Broad cloud-like fringes 

(6) Closely spaced fringes 
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(7) Lack of a known reference position 

(8) Unknown fringe sign 

(9) Uneven background 

(10) Regions blocked by opaque objects 

(11) Speckle noiae 

Therefore, considerable preprocessing is sometimes required before the task of 
locating fringe position or measuring phase difference is practical. 

Signal-to-Nolae Improvement - In addition to noise in the interferogram, electronic 
noise in the sensors also contributes to the final results. One of the most direct 
methods of improving aignal-to-noise ratio is signal averaging. Image analysers can be 
programmed to average many data frames at TV rates, thus improving picture quality. 
Next, spatial smoothing is useful to reduce high frequency noise due, for example, to 
speckle-like noise. 

Thresholding and Nonuniform Background - Many of the techniques commonly used in 
pattern ancl 1 ine recognition fail in interferometry because the intensity is sinusoidal 
and the background is nonuniforra. This presents a problem in thresholding. Procedures 
for dealing with nonuniform background Include the following: logarithmic detection, 

normalization and referencing to a mean intensity across the interferogram, which is 
first established before thresholding begins, normalizing to local Intensity averages, 
and, finally, subtracting off the nonunlform background. 

A floating threshold method, used by Becker 84 *, requires two passes over each line 
across the interferogram. The first pass searches for extrema. With the found extrema, 
a step-like threshold function is then defined, which is the mean gray value between 
adjacent extrema. The aecon d scan is referenced to this function to produce a binary 
fringe pattern. This produces good results on fringe fields, even with severe background 
intensity variation and low contrast. 

The problem can also be approached optically, in double reference wave holographic 
interferometry, by phase shifting one of the reconstructed waves by 180 degrees to 
subtract away the background entirely. 

Interactive Preprocessing - Until artificial intelligence is applied to this 
problem"] there is little hope of producing a fully automated system. Image processing 
systems fortunately allow a direct interaction from the operator to connect broken 
fringes, give the correct sign to fringes and to guess where a diacontinuoua fringe 
should proceed. What can be automated is partial recognition and flagging of the 
occurrence of such events to allow an operator to intervene. 

Coordinate Extraction - After the above preprocessing, the fringe pattern is a 
binary pattern. Then this pattern can be reduced to a line pattern where lines are the 
centers of the black or white fringes. Alternately, the edges or transition from black 
to white can designate the fringes. Each scan line, therefore, provides a y coordinate 
(for the line), an x coordinate for the fringes. With each x coordinate, a fringe number 
is required. Each adjacent scan must then be referenced to the previous scan, so that 
like numbered fringes can be identified. 

Because fringe density becomes very large in some regions, a magnification must be 
performed to maintain accuracy. This further complicates the numbering of fringes. A 
number of rules exist for fringe numbering. 

(a) Adjacent fringes differ in number by zero or one. 

(b) Fringes of different numbers cannot touch each other. 

(c) Fringes cannot end inside the field of view (except where they close on 

themselves). 

(d) When integrating over a closed line, the fringe number differences add up to zero. 

Phase shift interferometry and heterodyne interferometry both offer promise for 
making automated data reduction of holograms easier yet. The output signals appear to be 
more amenable to filtering and noise removal. This is one of the most promising and 
rapidly growing areas in holographic interferometry today. 
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CHAPTER 4 
LASER VELOCIMETRY 

4.0 BACKGROUND 

Laser velocimetry (or aneraometry) baa become one of the moat popular laser 
diagnostic tools in the field of aerodynamics. Various methods are routinely employed in 
virtually all aerodynamics research facilities, universities, and companies using fluid 
mechanics. Organisations like ONERA (France), DFVLR (Germany), NASA (USA), ISL 
(Germany/France), RAE (United Kingdom), and the military services supporting aerodynamics 
have invested heavily in elaborate, computer controlled LV systems for use in wind 
tunnels and flow research facilities. These have undergone continuous evolution now for 
over twenty years. 

The first principles of the methods and applications to simple flow fields are 
simple and straightforward; however, the general analysis and rigorous interpretation of 
the physics involved is extremely complicated, and a number of controversies still exist 
in the field. Among such controversies are the methodologies for proper seeding of the 
flow with scatter centers (SC) that accurately track the velocity fluctuations and the 
interpretation and removal of biasing errors. 

Most laser velocimeters can be described as time-of-flight instruments. They 
function by measuring the time taken for a SC to traverse a structured and fully 
characterized amplitude or intensity distribution in space. The simplest structuring of 
light is, perhaps, the one in which two focused spots, with an accurately known spacing, 
constitute the sample volume (Figure 4.0) first suggested by Tanner 1 *2 and reduced to 
practice by Schodls-6 and Smart. 7-1 ° The task is then to determine the time required for 
a scattering center to traverse the two points. These systems have been called two-spot 
systems, laser transit anemometers (LTA), and laser transit velocimeters. 



Figure 4.0 Illuminated volumes a) LDV, b) LTA. 

A somewhat more complicated configuration exploiting the Doppler effect to 
determine velocity of the SC. As will be shown, this can also be described as a tlme-of- 
flight method. The laser Doppler velocimeter (LDV) or anemometer (LDA) can take on many 
different configurations. Most of these can be modeled by the so-called fringe model, 
which will be discussed in detail. The sample volume is structured by crossing two or 
more focused coherent beams of light. The crossover volume contains a set of well- 
defined interference fringes that can be used to track an SC as it traverses the 

volume. The set of fringes will be "real" or "virtual", depending on the configuration 
used. The bgsls for these systems was set forth in the works of Yeh and Cummins * 1 and 

was reduced to practice by many workers such as Rudd 12 , Durst and Whitlaw 1 *, and 

Brayton 1 *' 1 ® and Mayo. 17 Stevenson, et al 18 have published a bibliography of earlv 

works. 

4.1 Laser Doppler Methods 

4.1.1 Theoretical Basis 

4.1.1.1 Doppler Theory - When a coherent beam of light of frequency, f 0 , strikes a 
moving SC^ an observer sees scattered light with frequency f Q + f D , where fn is the 
familiar Doppler frequency shift of the light. The shift occurs because the SC has a 
velocity component in Che Jc Q direction (Figure 4.1), causing it to receive wavefronts at 
a faster or slower rate tnan f Q ; and also because it has a velocity component with 
respect to the observer, causing either a compressing or spreading of the scattered 
wavefronts in the k s direction of the observer. 
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<•> 



Scatter Center Illumination 



Figure 4.1 Doppler shift of radiation by a moving scatterer. 


The observed phase of a wave, scattered by a stationary SC, will be changing at a 
rate u 0 (or 2af Q ) rad_ians per second. However, the phase at a moving SC is seen to 
change at a rate Wo - v . it and the observer sees still a different rate of change of 
phase given by Wo - v • k Q ° v • k . The Doppler frequency is, therefore. 


f - f - 


v 

“27 


(V l o) 


(4.1) 


Note that the Doppler frequency is proportional to the velocity component along the 
vector difference of the illumination and observing directions. 


EXAMPLE: An SC traveling 100 cm/sec is illuminated normal to its velocity (a * o) with 
an argon laser (x ■ .5 microns) , and is observed at an angle, b , of 60 degrees. The 
Doppler frequency is 100 cm/sec x cosine 60/x , which is about one MHz. 


An alternative method is to illuminate the sample volume with two beams separated 
by some angle. The two beams are Doppler shifted by different amounts, and the detected 
Doppler difference can be used to determine the velocity component along the vector 
difference between the two illuminating directions. Equivalently, the scattered Doppler 
shifted radiation can be detected at two different angles and, as will be shown, this 
Doppler difference is proportional to the velocity component along the vector difference 
of the two directions. We will show how to analyze each of these cases by both a 
Doppler/heterodyne analysis and a fringe model analysis. The systems that collect the 
scattered light, and measure the amount of the Doppler shift with reference to the 
unshifted light are referred to as reference-scatter (RS) systems. Those that use two 
components of scattered and Doppler shifted light are called Doppler difference or dual¬ 
scatter (DS) systems. 
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The remaining problem is to measure the Doppler frequency. Optical frequencies lie 
in the range of 10 14 to 10 15 Hz and cannot be measured directly with electronic devices. 
Essentially, two ways are used: 1) interferometry, and 2) heterodyning (or "beating down" 
the frequency, by mixing a second wave of slightly different frequency with it to produce 
a lower frequency signal that can be analyzed electronically). 

4.1.1.2 Fringe Models - In wave theory, the phenomena that involve frequency shifts can 
often be described not only with time domain equations, but equivalently with spatial 
domain, phase shift, or interferometry equations. The Doppler effect is such a phenom¬ 
enon, and the so-called fringe models for LDV are examples of interferometric analyses of 
the Doppler effect. Fringe models have played an extremely important role in the 
evolution of LDV because they have simplified the understanding, design, and use of the 
technique. Consider the configuration of Figure 4.2, in which the SC is illuminated by 
two focused, crossed beams that have equal intensities and Gaussian profiles. Referring 
again to Figure 1.3, it will be recognized that the crossover region will contain a set 
of stationary fringe planes. Figure 4.3 summarizes the character of the fringes and the 
sample volume. As an SC passes through the fringe planes, it scatters light whose 
intensity is modulated at a frequency 


f d “ v y/*)sin(e/2) 


(4.2) 



Figure 4.2 Dual scatter laser velocimeter, heterodyning at the scatter center. 


V Gaussian Radial 
^Intensity Distribution 



Figure 4.3 Enlarged cross-sectional views of the beam crossover region. 

The so-called geometrical sample volume, if determined by the transmitter optics, 
is represented by the volume wherein good fringes exist. Along the z-axis and in the xv 
plane, the fringe contrast is ut»>'tv. For laser beams with Gaussian profiles, the optical 
intensity varies in such a manner that the fringe contrast in the sample volume is not 
unity anywhere else. A point SC scatters light uniformly in all directions, and it 
scatters proportionately to the intensity to which it is exposed. The form of the 
detected signal is plotted in Figure 4.4, for traces of a point SC traversing the sample 
volume in the xz plane. 

The geometrical properties of the sample volume can easily be determined from 
Figure 4.3 when the transmitter is the defining system. The surface on which the 
Intensity falls to 1/e 2 of its maximum value is an ellipsoid. 
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EXAMPLE: 9 - 6 degrees, x - .5 microns, lens f/no f t /2b - 200. Then the sample volume 
diameter is .013 cm by 25 cm long, and the frequency to velocity conversion factor is 
2 kHz/cm/sec. 



Figure 4.4 Intensity distribution in plane of crossed beams. (Assumes beams of equal 
intensity and polarization vectors linear and parallel). 

The geometrical properties of the sampling volume can also be determined by the 
receiver. This is easiest to see when an off-axis detector is used. To determine the 
actual geometrical sampling volume, one can image the receiving aperture into the 
transmitter-formed sampling volume. The signal can come only from the overlapping 
region. 

So far, we have referred to the geometrical sampling volume as being determined by 
the transmitter or receiver. The actual volume from which data can be taken is much more 
complicated than this. Consider a case in which SC of a single size are used. The 
volume from which data is taken depends upon the detectability of the scattered light. 
This is determined not only from the geometrical properties of the system, but also from 
the electronic properties of the detecting system, for example, by where the detection 
threshold is set or by other restrictive signal processing. The actual sample volume can 
be modified electronically by deliberately placing conditions on the signal. For 
example, the signal can be required to contain a set number of fringes above some 
threshold. One might also choose to accept only those signals having contrast greater 
than some value. Clearly, the actual sample volume depends upon the SC size for anv 
geometrical configuration. 

Consider the system of Figure 4.5. Here the sample volume is illuminated by a 
single beam, but light is collected at two different angles. No real fringes could be 
observed in such a sample volume. Nevertheless, it is useful to introduce the concept of 
virtual fringes, fringes that would exist if the two received beams of light were back- 
projected to form a sample volume. We will show that the signal from such a svstem is 
accurately described with this model for a wide range of system conditions. 

An interesting reciprocity relationship exists between the laser and photodetector. 
The laser and photodetector position can be switched for any given system geometry 
resulting, in an alternative configuration with similar geometrical properties. 
Figure 4.2 and 4.5 are examples of reciprocal systems. 

The parameters of an RS system can be determined to a large extent by using the 
same fringe model. The method is to examine a virtual sample volume that is formed by 
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Figure 4.5 Dual scatter laser velocimeter, heterodyning at receiver. 

the virtual crossover of the input beam and the scattered beam. The scattered vector 
changes within the sample volume, however, so that the Doppler frequency varies over the 
sample volume leading to the expression "Doppler Ambiguities". This volume will contain 
a set of virtual fringes of the proper spacing and orientation to evaluate the system. 
However, it is more accurate to examine fringes at the place where mixing effectively 
takes place. In some DS systems, and in all RS systems, this occurs at the photo¬ 
detector. Therefore, the signal from such a system can be accurately modeled as an 
interference fringe that is formed by the mixing of a scattered, Doppler-shifted wave 
with a second reference wave producing a beat frequency. This model is useful in 
determining a variety of instrumental effects such as misalignment and poor quality 
optics. 

4.1,1.3 . The Moire Analogy - To augment the fringe model. Durst and Stevenson!* adapted 
the moire/interferometry analogy to LDV. By using the grids supplied in Appendices 1-4, 
many of the characteristics of LDV sample volumes and signals can be illustrated. Most 
of tjie effects that can be modeled with the fringe model can be demonstrated with the 
moire analogy, (see Chapter 2). Practical effects such as beam misalignment, vibration, 
frequency modulation of one of the beams, multiple SCs in the sample volume, wavefront 
distortion, and heterodyning of signals can all be lucidly illustrated. These and other 
effects are left as an exercise for the reader. 

4.1.2 Measuring the Doppler Frequency 

4.1.2.1 Direct Optical Frequency/Wavelength Measurement - As shown in Chapter 2, 
interferometry can be employed to determine optical wavelength and, thus, to measure 
Doppler shift. The first interferometer to be adapted to this measurement was a Fabrv- 
Perot interferometer. 20-21 Figure 4.6 illustrates the configuration. The interferometer 
is swept through its free spectral range by piezo-electric elements that vary the mirror 
spacing. The interferometer is illuminated by both Doppler-shifted and unshifted 
radiation, the latter serving as a reference. A velocity measurement is made each time 
the interferometer passes the Doppler-shifted light, typically up to kHz rates. The best 
available sensitivities for these interferometers are in the order of 100 MHz, limiting 
their use to higher velocity measurements. This limitation has apparently discouraged 
the use of Fabry-Perot systems to any large extent. 

To partially overcome the restrictions encountered with the Fabrv-Perot interfer¬ 
ometer, Smeets and George 22 “ 30 developed a more sensitive Michelson (Section 2.4.1) 
interferometer/spectrometer, that allows the direct measurement of Doppler frequencies to 
resolutions limited only by the stability of the laser Itself. The system is Illustrated 
in Figure 4.7. A single laser beam illuminates a sample volume which contains moving 
SCs. The radiation scattered by the SCs ig Doj>pler shifted by an amount that is 
proportional to the velocity component in the k - k direction. The scattered radiation 
is collected by a lens, and is then conducted 8 to the interferometer by a single mode 
optical fiber. 

The light bundle leaving the optical fiber is collimated with an objective lens. 
After having been linearly polarized, it passes through a Pockels cell and then enters a 
Michelson interferometer. The Michelson interferometer has a number of important 
features: 1) its two legs have a path difference, dl ; 2) it employs a polarizing beam 
splitter to direct one polarization into the short leg and the other into the long leg; 
and 3) it employs a second polarization beam splitter to produce two interferences that 
are complementary to each other. The relative phase of th * two complementary signals can 
be adjusted so that they cancel. A change in wavelength of the incoming light wave can 
be exactly compensated for by a phase shift produced by the Pockels cell. In this 
manner, the velocity is proportional to the voltage required on the Pockels cell to 
maintain the null condition. 

The velocity range of this technique has essentially no upper limit, but, as 
discussed in Chapter 1 (Section 1.5.2), lasers have inherent frequency Instabilities that 
will limit the lowest velocity measurable. With some degree of stabilization, frequency 
can be stabilized to one part in 108, representing a low end value of 6 MHz for the 
0.5 micron wavelength. This translates into a minimum velocity of 3 meters/second 
without additional laser stabilization. 
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(•) Schematic Diagram 



(b) Typical Oscilloscope Trace 


Figure 4.6 Laser velocimeter employing Fabry-Perot interferometer. 

This LDV represents a significant advance and should open the wav to a wide variety 
of new applications. Unique advantages include the following: 


• No upper limit on velocity range 

• Fast, direct, and continuous velocity readout 

• Works well with heavily seeded flows 

• Works on large particles and surfaces 

• Much simpler transmitter, less affected bv density fradlents and lower laser 
coherence than other methods 

• Easily adaptable to fiber optics 

4.1.2.2 Heterodyning Methods - Equation 1.7 describes the intensity of the sum of two 
waves that have different frequencies. Namely, when two such waves are combined, the 
resulting time averaged intensity, averaged over times that are large compared to the 
optical period, is modulated by a beat frequency that is the difference between the two 
original frequencies. If the frequencies are sufficiently close to each other, then the 
resulting beat frequency can lie in a range that is measurable by conventional 
electronics. This is called heterodyning and la routinely used in electronic slgnsl 
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Figure 4.7 Michelson Spectrometer LDV of Smeets and GeorRe. 

processing. Therefore, to measure floppier shift electronically, one can mix the laser 
beam coming from the laser, with the Doppler shifted, scattered radiation to produce a 
beat signal at the Doppler frequency. Commercial electronic equipment now measures up to 
about 200 MHz, which covers, essentially, the velocity range of interest to 
aerodynamicists. 

The different methods of heterodyning and the place at which heterodyning actually 
occurs create distinctively different systems and system requirements. The most 
important of the distinctions is that between systems which employ coherent detection or 
incoherent detection. Figures 4.2, 4-5, 4.8, and 4.9 show the differences. 



Figure 4.8 Reference scatter laser velocimeter, heterodyning at the scatter center. 

With the systems of Figure 4.8 and 4.9, the Doppler shifted radiation is mixed with 
unshifted radiation. This has been termed the reference/scatter (RS) mode. The amount 
of the Doppler shift is given by Equation 4.1 for both cases, but there is an important 
difference in the actual physics of the two cases. The system of Figure 4.8 has two real 
wavefronts entering and creating real fringes in the sample volume. Anything that 
happens to the radiation after leaving the sample volume will happen identically to both 
beaus, since they are overlapped. In this sense, the heterodyning has already taken 
place at the sample volume and the system can be treated as an incoherent system 
thereafter. 
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Figure 4.9 Reference scatter laser velocimeter heterodyning at receiver. 

Contrast this with the reciprocal system of Figure 4.9. No real fringes exist in 
the sample volume. The alignment and coherence of the two beams is extremely critical 
all the way to the photodetector where heterodyning occurs. 

Figures 4.2 and 4.5 illustrate the types of systems that heterodyne two different 
Doppler shifted waves to produce a beat frequency. If the two different portions of the 
wavefront have Doppler shifts equal to 

f D r 0/2 ” • ( k 9l - & o) 

f D2 - W2, . (k s2 - k 0 ) (A.5) 

then the associated intensity variation of the mixed waves is modulated by the difference 
frequency 


V 0/2 - • (*.1 - ' k „ 2 ) (4 - 6 > 

which is proportional to the velocity component lying along k . - k y direction, which 
is independent of the illuminating beam direction. This is the so-called dual-scatter 
(DS) mode, also referred to as the Doppler difference mode. Note that the system of 
Figure 4.2 is input aligned, while that of Figure 4.5 is output aligned. The DS system, 
with real fringe or input aligned optics, has become, by far, the most widely used 

configuration for a variety of reasons to be discussed. 

4.1.2.3 Signal Analysis and Interpretation 30 *"^ - Since the elementary form of the 
Doppler signal does not usually exist in practice, we must expand the discussion to cover 
signal characteristics that are commonly encountered. Most of the complicating effects 
can be described in terms of the fringe model. 

The point SC, passing through the sample volume of a DS system, scatters light in 
ail directions with an intensity form illustrated in Figures 4.4 and 4.10, that depends 
upon where it traverses the sample volume. A path through the center of the sample 
volume produces a sinusoidal signal enveloped in a Gaussian. The Fourier transform of 
the signal contains two major components, first, a lower frequency part of the spectrum 
that makes up the Gaussian part of the signal, that is commonly called the signal 

"pedestal" or sometimes (incorrectly) the DC part of the signal and, second, the 

frequency of the sinusoid which is the Doppler part of the signal. It is important that 
these two parts of the signal be separable, and this can be designed into a system. 

The signal form can be altered by many practical effects presenting the following 
questions: 

(1) How is the signal affected by 

(a) finite size SC? 

(b) more than one SC at a time? 

(c) different velocity SC at one time? 

(2) When does the fringe model break down? 

(3) What are other broadening effects? 

If the SC has a finite size, then the scattered light signal can never be zero, 
since some part of the particle must be illuminated at all times, even if another part of 
the particle is in a dark fringe. This reduces signal contrast or signal "visibility". 
The visibility parameter has actually been used as a particle sizing parameter, but it 
has not proven reliable and has been abandoned for this purpose by most investigators.3i 

A second effect of finite particle size is the failure of the fringe model to 
predict the signal character for large angles between the two illuminating beams. Finite 
sized SCs exhibit a large forward scattering component, so that the light scattered from 
two wide-angled beams will differ in every viewing plane except the one that bisects the 
angle between the two beams. This will lead to a second reduction in signal contrast 
that is not predicted from the fringe model. 
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Signal Fourier Transform 



Sources of Signal 


Point S.C. through sample 
volume center 


1. S.C. size comparable to fringe space 

2. Unequal beam intensities 

3. Wide angle between beams plus large S.C. 


1. Electronic noise 






1. Multiple S.C. 


Figure A.10 LDV signal types - dual scatter system. 

Inherent in the discussion so far Is the requirement that the SCs do actually 
follow the flow. Large SCs may not be capable of following the flow. The subject of the 
seeding of the flows has become so dominant in reliable LDV that this will be treated 
separately. 

More than one SC In the sample volume affects the two different LDVs, reference 
scatter and dual scatter, differently. This is further compounded by the existence of 
different size SCs, some of which follow the flow and some of which do not. Finally, the 
actual fluid velocity Itself can change in space and time within the sample volume. All 
of these result in Doppler frequency broadening, and usually in a decrease in the signal* 
to-noise ratio. Typical signal types for reference scatter systems are illustrated in 
Figure 4.11. 

Isolating broadening effects, to determine which ones are caused by the instrument 
and which ones are caused by the flow, is one of the major problems of the LDV signal 
interpreter. For example, optical system broadening occurs when the fringes are not 
parallel and equally spaced in the sample volume. This can be caused by poor optics or 
by system misalignment. Electronic broadening is caused by the instrument incapability 
to respond to the required frequency. These are usually characterized by first using the 
instrument on a known, fixed velocity. 

The signal-to-noise ratio is an important parameter to characterize an LDV. This 
can be estimated by assuming that the major noise source is the photodetector dark 
current noise. In general, other sources of noise, such as laser noise and ambient light 
and electronic noise, must be considered. The average Intensity of light striking the 
SCs Is given by 

I 0 - 8P/,Dj (4.7) 

where 

P is the laser power 

Df is the diameter of the laser beam at focus 

The amount of scattered radiation is equal to 

Ig - »(9)I 0 N sc /K Z (4.8) 



where 






109 


Signal 


Fourier Transform 





Single point SC through sample 
volume center 


i. Multiple SC in sample volume 


1. Electronic noise 


Figure 4.11 LDV signal types - reference scatter systems 
o(e) is the scattering cross section 
N sc is the number of SCs 
R is the distance to the collector. 

Brayton** has calculated the S/N for the two types of LDVs as the following: 

(S/S), s - 

. ! h,(9)p K 8C [ie)‘ (4 ' 9) 

4hcBx 

(S/N)pg ■ 


where 


ge is the scattered light collection angle 
n is the detector quantum efficiency 
B is the system bandwidth 
h is Planck’s constant 
c is the light speed 
X is the wavelength 

Especially noteworthy is the S/N improvement with cclle .ting aperture which favors 
the DS system. Note also that the RS system S/N improves with the number of SCa. 
Finally, it must be noted the S/N decreases with velocity because the required system 
> bandwidth, B, increases with velocitv. 

j By their very nature, LDVs produce discontinuous, digital-like data because one 

j velocity measurement is produced per particle transit at the very most. These data 

points must be properly treated statistically to arrive at meaningful interpretations, 
f For example, in a highly turbulent flow, more high-speed data samples will be collected 

? simply because they arrive at the sample volume more often, so one might statistically 

f bias the data to an incorrectly high mean velocity. This bias error is easy to correct 

I and the time averaged velocity can be calculated by 



A. 
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<V> - E V.t ./T (4.11) 

1 1 1 

where Is the velocity of the I th particle, t^ Is the transit time of the I th particle 
and T la the average transit time. On the other hand, the reduced S/N for high-speed SCs 
can bias the data towards low velocity. 

If a particle sl 2 e distribution exists and if all of the particles do not track the 
flow, the problem of interpretation of turbulence data is virtually impossible. 

4.1.3 Optical Geometries 

4.1.3.1 Systems - The ootical requirement of a DS system is that a sample volume of 
desired dimens ions be produced by focusing and crossing two laser beams at a desired 
point in space, followed by collection of the maximum amount of scattered light from the 
volume and from nowhere else. This can be reduced to the following optical functions 
(refer to Figure 4.12 and 4.3). 



Figure 4.12 Laser velocimeter geometry. 

(1) Referring to Chapter 1, determine the F/number required to produce the (transmitter 

defined) sample volume dimension and location, then expand the laser beam to the 

required input diameter. This is called beam conditioning. 

(2) Produce two parallel, collimated beams of light, separated bv the desired amount 
(see Figure 4.12). 

(3) Focus and cross the two beams at their focus, by passing through a high quality 
lens or mirror or by some other means. 

(4) Collect the light scattered from the sample volume with an efficient receiving 
system that is stopped to further define the geometrical sample volume. 

A wide variety of beam splitter elements have been used to produce the two parallel 
beams of light. Figure 4.13 illustrates some of the most commonly used elements (see 
also Figure 4.20). With elements like 4.13a-b, the angle between the two output beams 

can be adjusted by rotating one of the components. Although sometimes useful, this can 

add to the system alignment problem. The remaining beam splitters automatically 
transform a single input beam into two parallel beams (if the components are accurately 
made). Note that, of the remaining components, only c, g, and h have equal pathlengths 
for the two beams. This can be a critical point for lasers with a short coherence 
length, and pathlength compensation may be provided by additional components. For 
example, (b), (d) and (e) can be used in combination to produce two collimated beams of 
equal pathlength, as is shown in Figure 4.14. 

Other types of beam splitters produce two waves that diverge from a point, and 
these waves can be reconverged to produce a detection volume that contains Interference 
fringes. The two beams can be produced by an amplitude or phase grating, or by a 
holographic optical element (see Figure 4.15). This can be understood simply as the 
imaging of the grating into a detection volume. If the grating moves, then the fringes 
likewise move, providing a useful mechanism for determining the direction of the velocity 
of the scattering center. An SC that moves against the direction of the fringe motion is 
associated with a higher frequency than one that moves in the direction of the fringe 
mo*’* on. 


Aeonsto-optical cells can be used to spilt the beam into two components, each 
having slightly different frequency. When the two beams are reconverged, the fringes 
will move as described in Chapter 1. This can also be understood as the Imaging of an 
initially moving set of fringes. 






Ill 



(a) Prisa Pair (b) Beam Splitter Mirror (c) Double Plate 



(d) Parallel Plate 



<g) Grating 



(e) Parallel Plate (f) Prlaa Pair 


Hoe Hoe 



(h> Prism Pair 


Figure 4.13 Beam splitting Into two parallel beams* 







FiRure 4.14 Pathlength compensation 




Figure 4.15 Beam splitting with gratings. 


Photcdetectors 



Figure 4.16 Beam splitting with a Wollaston Prism 

Another beam splitter in this class is the Wollaston prism, which produces two 
beams that are polarized normal to each other; so, when the two beams are reconverged, 
interference does not occur. The procedure is illustrated in Figure 4.16. Polarizers 
can be placed with the polarization direction at 45 degrees with respect to the two 
crossing beams so that interference can occur and virtual fringes are produced. If two 
different polarizers are used with a 90 degree difference in their respective 
polarization directions, then the two resulting beams will produce signals that are 
180 degrees out of phase with each other. These two can be mixed electronically to 
remove the pedestal from the signal. 
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4. 1.3. 2 Two Component Dual Scatter Systems - Figure 4.17 illustrates the basic two 
component DS system. One can interpret the process as the production of two sets of 
fringes in Che detection volume. The DS system can take on a second form illustrated in 
Figure 4.18. In Figure 4.18, the scattered light is divided with a beam stop into two 
sets of beams, each having an associated set of virtual fringes. In Figure 4.17, the 
fringes are real. 

The only added problem, in going from one to two components of velocity, is the 
separation of the two components upon detection. The separation can be done by polariza¬ 
tion, by color, by electronic means, or even by coherence discrimination, but in recent 
years most Investigators have applied color separation. Argon lasers are Ideally suited 
for this process, since two dominant wavelengths (488 nm and 514 nra) are of nearly equal 
power in available lasers. In Figure 4.17, beams 1 and 2 are of two different colors, 
while beam 0 is a mixture of the two. 

In a two component system, the orientation of the fringes, with respect to each 
other and with respect to the flow direction, can be optimized for most cases. If the 
flow direction is approximately known, then the fringes should be oriented symmetrically 
about this direction. Neither fringe set should be placed with fringes parallel to the 
velocity direction, for then the SC will not cross any of the fringes. The optimum angle 
between the fringes is rarely 90 degrees, and is more likely to be about 60 degrees to 
get the most accurate data for both components of velocity. 

4.1.3.3 Three Component Dual Scatter Systems - Three component systems present a number 
of unique and difficult conditions for the ’designer. The problem is to produce three 
independent sets of fringes in a sample volume, and then to separate the scatterd light 
signal into components that can be associated with the correct velocity component. 
Practical difficulties arise from the lack of a good laser candidate that has three 
different wavelengths of equal intensity, and from the need to orient the fringes so that 
the SC crosses a sufficient number of fringes in each of the three sets of fringes. 

The problem of producing three fringe sets of equal intensity has generally been 
solved, either by using two lasers (for example, an argon-ion laser and a krypton laser) 
to produce three wavelengths 33 ; or, by tuning down the power in the two principal wave¬ 
lengths of argon lasers and choosing a third, say at 470 nm (see Figure 4.19). Either 
method leads to a rather difficult and expensive configuration. Pfeifer 34 (1985) has 
presented a novel approach to this problem, in which the argon laser is used partly to 
pump a dye laser for one of the components, while using the two principal wavelengths for 
the other two components. 
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Aq — 514.5nm 
A b = 488nm 
A v = 476.5nm 


L,: 8 W all lines 

L,; 3 W on K = 476.5nm 


Figure 4.19 Scheme of ONERA 3D IV 


Lj : 15 W argon laser operated at 8 W all lines 
1*2 15 W argon laser operated at 3 W of a • 476,5 nm 

Dj: dichrolc plate transmitt ing only the green line 
D 2 : dichrotc plate transmitting only the blue line 
C: Bragg cell 
BS; beam splitter 
A. afocal telescope 
G: green (\ - 514,5 ntn) 

B: blue (a * 488 nm) 

V: violet (a * 476,5 nm) 


4.1.3.4 Moving Fringes * In any of the geometries described above, a bias frequency can 
be added to the Doppler signal. In terms of the fringe model, this can be described as 
producing a system with moving fringes, so that for a stationary SC in the sample volume, 
the scattered light will be modulated at a frequency determined by the rate at which the 
fringes move over the SC. Fringes that are produced by interfering beams of two 
different optical frequencies will move at a rate that Increases with the frequency 
difference. Any point in the fringe pattern will vary sinusoidally at the frequency 
difference of the two waves. Moving fringes can be produced by several different 
methods. A frequency shift can be added to either or both of the two waves with a Bragg 
cell. If the fringes are produced by using a grating, then the fringes can be moved by 
moving the grating. In practice, this is implemented by using a radial grating on a 
rotating disk. 

4.1.3.5 Fiber Optics - Fiber optics can be used to conduct the light from the laser to 
the transmitter.I? the transmitting fiber is a single mode fiber, then the end of the 
fiber is treated as a point source, just as it would be if it had been produced by 
focusin« the direct beam from the laser. It is then collimated and split in any of the 
ways described ^bove. Alternately, light being conducted in a fiber can be split, as 
shown, by a fiber coupler into two beams that can subsequently be collimated into 
parallel beams. (Figure 4.20) 

Single mode fibers can conduct no more than about 100 milliwatts of power. Above 
this power, nonlinear optical interaction with the material can cause prohibitive 
losses. If greater power lasers are used, then the light must be transmitted through 
coherent fiber bundles. The output beam is no longer a simple diverging wavefront, but 
contains speckles. This somewhat complicates the sample volume, first by enlarging it 
and by producing patches of fringes that are not in phase with each other. 

The scattered light can be collected by a lens that focuses the light into a fiber 
conducting it to a remote location into a photodetector. When incoherent detection or 
heterodyning at the SC is involved, this fiber can be a much larger fiber since coherence 
is no longer required (Figure 4.21). 


. ? 
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FiRure 4.21 Fiber-optic receiver 

4.1.3.6 Reference Scatter and Other Direct Doppler Systems - Direct Doppler systems rely 
on the measurement of the Soppier shift of the scattered llRht (as opposed to the 
difference In Doppler shift used In dual scatter systems). Only one beam need he 
projected Into the flow, and the geometrical sample volume, which lies alonR this 
transmitted beam, is defined by the overlap of the imaRe of the recelvinR aperture 
projected across the transmitted beam (Illustrated in Figure 4.22). The properties of 
such a system can actually be approximated by ualnR the frlnRe model of 4.1.1.2, wherein 
the virtual fringes are produced by an imaRlnary beam of llRht, originating at the 
receiver aperture and projecting back through the system optics to cross the transmitted 
beam from the laser (i.e., the reciprocal system). 


PNOTOOETCCTOR 



Figure 4.22 The sample volume of a reference scatter LDV. 

To determine the amount of the Doppler shift of the scattered light, either an 
optical spectrometer/Interferometer can be used directly or a heterodyning technique can 
be used to beat down the frequency of the signal to a frequency that Is measurable 
electronically. 

The heterodyning procedure requires that a reference wave be mixed with the 
scattered light. An optical system that is similar to the DS systems already described 
can be used to achieve this as shown in Figure 4.12. The beam splitter is designed to 
produce a reference wave that is typically 5 to 50 times as bright as the scattered wave, 
and, for this confIguration to perform well, the reference wave must pass through the 
flow field largely unscattered and unmodulated. This arrangement is particularly 
convenient In terms of alignment, since the reference wave, which passes through the 
sample volume, Is automatically aligned with the scattered wave. This is sometimes 
referred to as a self-aligned or input aligned system. 

In highly seeded or highly turbulent flows the configuration of Figure 4.12 may 
not be applicable, since the reference wave can contain a large quantity of scattered 
light and phase noise that can reduce the signal-to-nolse ratio. Therefore, In some 
cases the signal can be optimised by passing the reference wave around the flow field. 
Precise alignment of the reference wave with the scattered wave in such a system, is 
extremely difficult. As described In Section 4.1.1, It Is Important that the reference 
wave and scattered wave overlap precisely over the collecting aperture. This requires 
the F-nuifcers of the two waves to be the same. This Is a simple design problem for the 
system of Figure 4.12, but Is more difficult when the reference and scattered waves 
emerge from two different points. 
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When a spectrometer is used, the only requirement of the receiving system is that 
the optics be focused to collect light from the focal region of the transmitting system. 
The crossover region of the transmitter and receiver optics defines the geometrical 
sample volume. Some of the possible geometries are illustrated in Section 4.1.2. 

4.1.4 Electronic Signal Processors 

4.1.4.1 08Cilloscope8 - The simplest form of signal processor is an oscilloscope that is 
set to trigger at some threshold. If the signals generated by different SC are similar, 
then a sine wave with a Gaussian envelope will be observed, and the frequency can be 
determined simply by measuring the scope screen. If the signals from different SC are 
not similar, then a storage scope must be used. The scope is set to trigger on the 
smallest possible burst signal, allowing the examination of individual "Doppler" 
bursts. The bursts will have the appearance shown in Figure 4.10. 

Using an oscilloscope is a good wav to evaluate the quality of the overall system, 
but it is not a practical method to collect a useful amount of data. The signals should 
exhibit high contrast and low noise. The measured amplitude will help in setting the 
threshold of other types of electronic processors. Some processors actually use the 
scope trigger to trigger the processor that analyzes the Doppler bursts. Low contrast 
signals are a symptom of poor system alignment, incorrect beam intensity ratios, too 
large SC, or insufficient electronic frequency response, or combinations of these. 

4.1.4.2 Spectrum Analyzers - The first sienal processors to be used in LDV were spectrum 
analyzers. Even though they are simple to use, their limitations severely restrict their 
usefulness. The signal from the photodetector is fed directly to the analyzer, which 
sweeps a variable electronic filter over a preset frequency range, and presents an 
amplitude versus frequency plot on a cathode ray tube at rates ranging from seconds down 
to milliseconds. When used in LDV, the sweep rate must typically be quite slow to allow 
enough Doppler bursts to enter the system at each frequency. The system outDut is 
broadened electronically by the signal envelope, by noise, and bv phase fluctuations 
caused by more than one SC in the sample volume at one time. Therefore, it is difficult 
to apply the Instrument to anything other than laminar, constant velocity flows. 
Spectrum analyzers are useful to determine the frequency range in which other electronics 
must operate. 

4.1.4.3 Frequency Trackers - If the Doppler signal were merely a frequency modulated 
sinusoid with constant amplitude, conventional electronic frequency to voltage conversion 
would constitute an acceptable method of signal processing. However, except in some 
special instances, a conventional frequency tracker is incapable of dealing with large 
amplitude fluctuations, noise, and especially signal dropout. Therefore, the signal must 
usually be conditioned before it is analyzed, and signal dropout must be considered. 

Frequency trackers inherently require large amounts of seeding to provide a nearly 
continuous signal. They find their greatest use in the LDV of liquid flows where it is 
easy to provide the required seeding. They have virtually vanished from use in gaseous 
flows. LDV in gaseous flows is basically a binary process, while frequency tracking in 
the usual sense is an analogue process. Because of the limited use, they will not be 
discussed further and the reader is referred to the references for further 
information. 3 s-3 7 

4.1.4.4 Electronic Counters - Commercial counters have evolved to a relatively uniform 
set of electronic functions. 3 6 These are illustrated in Figure 4.23. The sienal 
entering the processor is a sinusoid with a Gaussian envelope plus a vast range of 
noise. The signa* Is first high pass filtered to strip away the pedestal, leaving a 
signal that oscillates about zero. Next, the signal is bandpass filtered to remove noise 
that lies outside the frequency range of interest. The signal-to-noise ratio is 
inversely proportional to the bandpass, so it pays to select Lhls as closely as 
possible. Validity checks are done at this point to insure that the signal has the 
correct appearance. For example, a condition for the signal to be accepted for 
processing might be that, after a maximum or minimum is observed, the signal must cross 
zero before another maximum or minmum is observed. This is called a "three-level 
check". The signal is further amplified and clipped so that it becomes a binary signal. 

In some systems, additional checks are done at this point to insure that the 
frequency of the signal does not change during a Doppler burst. This is done bv counting 
the time required for a set number of cycles occurring early in the burst, and comparing 
it with the time equired for the entire burst. When divided by the total number of 
cycles, this number should be the same to within a small percent. Otherwise, the signal 
should be rejected as invalid. 

4.1.4.5 Simultaneity Test - Multiple component measurements often require that valid 
signals be observed simultaneously from the different components for the data to he 
useful. Therefore, coincidence circuits reject signals that do not occur simultaneously 
In all of the components. Almost without fail, this severely reduces the data rate. 
Fairly simple electronic circuitry is used for this function. When a threshold is 
exceeded In either channel, a pulse Is sent to an AND gate in the other channel. The AND 
gate remains open as long as the pulse remains, so no signal is passed unless both and 
gates are open. A preset time window can be produced in this wav allowing the acceptance 
of the two velocities. By increasing the pulse duration, the data rate can be increased 
at the expense of precise simultaneity. The optimum pulse duration is a function of the 
turbulence frequency in the flow and, to be meaningful, the measurements must be made to 
within times that are small compared to the characteristic times for the flow. 
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Signal Rejected by "Three Level" Filter 




Signal Rejected by Period Test 


Figure 4.23 Signal processing sequence. 

A. 1.4.6 Transient Recorders - Today's digital transient recorders can digitize an analog 
signal at rates in excess of 200 MHz, so it is possible to store the entire Doppler burst 
in digital form.39“40 Then the stored signal can be processed by an ordinary computer. 
The memory of such a system will be filled quickly, limiting the use of the recorder to 
short times. 

A variety of methods have been devised to extend the use of such systems to longer 
times. Data compression schemes can insure that the recorder accepts data only when a 
Doppler burst is present. Buffer memories are also used to dump the recorded data at 
slower rates as the recorder memory is filled. This provides a sampling of the flow 
velocity with dead times lasting a few milliseconds for each data dump. 

Use of a transient recorder is especially important for two specific cases, one 
when the signal is extremely poor and requires additional processing to extract usable 
signals, and a second when the quantity of data is so small that every signal is needed 
to constitute enough data. Having the entire signal for processing at a slower rate 
permits the use of powerful filtering procedures such as correlation. Portions of the 
signal can be correlated with itself to improve the signal-to-noise ratio. Another 
procedure is to Fourier transform the signal to determine its strongest spectral 
content. Once the information signal is identified bandpass filtering can be applied 
around the frequency of the desired signal to further remove the noise in the superfluous 
frequency band. 

A. 1,4.7 Correlators - The previous section described how correlation is used to improve 
the signal-to-noise ratio of a stored signal. 41 High-speed digital correlators have been 
developed to perform this function on the live Doppler bursts. The fastest such machines 
can correlate signals of frequency in excess of 100 MHz. If the Doppler bursts are all 
of the same frequency as in a laminar flow, then nearly perfect Doppler bursts can be 
constructed from an extremely noisy signal, and the correlation time will last as long as 
the Doppler burst. If the frequency of the Doppler bursts varies as in a turbulent flow, 
Chen the correlation time will provide a measure of the turbulence frequency of the 
flow. Correlators are commonly used in front of burst counter processors. 

A. 1.5 Applications 

Hundreds of interesting and Impressive applications of LDV could be cited; however, 
the objective here is to illustrate the principles, the state of the art, and to describe 
some recent developments that have not yet reached textbooks on the subject. Therefore, 
the approach taken is to pick a few examples from the various laboratories visited during 
preparation for this book that illustrate the methods and system variations especially 
well. Examples are given for use in wind tunnels, shock tubes, turboraachinery, and in 
the atmosphere. 
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4.1.5.1 Wind Tunnels - ONERA (France) : Excellent examples of elaborate LDV systems are 
those that have been developed by ONERA in France 42- * 4 5 for use in the Modane and Le Fauga 
wind tunnel facilities (Boutier et al. 1983). The S3 Mondane facility la a blow-down 
transonic tunnel with a run duration of less than 200 seconds. The velocimeter is a two 
component forward scatter system using an Argon laser at 9 watts, blue and green fringes 
with spacing of 24.5 micrometers, oriented at 45 degrees with respect to the tunnel axis, 
and modulated with Bragg cells to 5 and 6.25 MHz respectively. The forward scattered 
light is collected by a Cassegrainean telescope that has a 1 .2 meter working distance at 
f/6. 


Figure 4.24 is a diagram of the systems, showing the extensive use of computers to 
command and track the independent traversing of the transmitter and receiver units in 
three dimensions. This system can take measurements at one point and then move to a 
different point, for a second measurement, in a few seconds. Incense seeding is done to 
increase data rate and to insure that the particles follow the flow. The system has been 
used successfully to study boundary layers, wakes, and flow profiles. 



Figure 4.24 General diagram of a two-component LDV system. 

In continuous tunnels, like the S2 Modane and FI Le Fauga closed circuit facili¬ 
ties, where more time is available for collecting data, backscatter systems can be 
used. Such a system is shown in Figure 4.25. At higher velocities, even here, forward 
scatter systems have been required. In pressurized systems such as the Le Fauga 
facility, where it is required to place the transmitter in the pressurized plenum, 
another problem occurs. At higher pressures, the Brewster angle in the laser changes due 
to the increase in the refractive index inside the laser cavity. This problem has been 
solved here by placing the laser in a sealed box. 

The most recently developed system of ONERA Is the three-component system diagramed 
in Figure 4.19. The system uses three beam pairs at different wavelengths. These are 
extracted from two different argon lasers. Again, the transmitter and receiver are 
traversed under computer control. This is one of only a very few three-component systems 
in operation. 

LDV systems that are similar to the ones in use at ONERA are in use at other 
research centers. A number of these are described in more detail in NASA Conference 
Publication 2243 (1982). 4 5 Some examples are noted in the following: 

NASA Ames Research Center (USA) : Two- component LDV systems are used in various 
tunnels such as the 13 by 26 meter subsonic, the 2.5 by 3.3. meter transonic, and the 
2 meter supersonic tunnels. Backscatter systems are primarily used. A long-range LDV 
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(a) 2-D backscatter LDV of ONERA. 



(b) LDV application in the FI Le Fauga Facility of ONERA 

Figure 4.25 Backscatter LDV in a transonic wind tunnel. 

system has been developed for the 13 by 40 meter subsonic wind tunnel, which is the 
largest such test facility in the free world. 

NASA Langley Research Center (USA) : Two-component backscatter systems for the 5- 
raeter transonic tunnel use a two-dimensional traverse and a zooming system for the third 
dimension. A two-component backscatter system is used in the 4 by 7 meter subsonic 
tunnel. This is a somewhat unique system in several respects. It employs what is prob¬ 
ably the largest two-dimensional traverse for LDV, covering a sampling volume of 1 by 2 
meters and zooming for the third dimension. The traverse alone weighs in excess of 
2700 Kg. A tilt and pan system also allows rapid scanning of the transmitted beam in 
addition to the traversing action. A 12-watt argon laser is used. These two systems are 
permanently dedicated installations. Other smaller systems are employed in tunnels 
ranging up to hypersonic Mach numbers. 

DFVLR (Germany) : A two-component backscatter system is used in the 3 by 3 meter 
subsonic runnel and in smaller facilit ies. 
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NASA Lewis Research Center (USA) : A two-component backscatter LDV has been used in 
the 277 3 meter transonic tunnel. This system uses a two-dimensional traversing 
system with zooming optics and a 15-watt argon laser. The system is fully enclosed in a 
water cooled, vibration isolation housing. 

Naval Surface Weapons Center (USA) : A two-component backscatter LDV is used in low 
velocity flows in a 1 bv K3 meter tunnel; a two-component forward-scatter LDV has been 
used up to Mach 5; and a three-component system has been used at Mach 3.** 7 "5i 

U.S. Air Force. AEDC (USA) : Two and three-component LDV systems have been used in 
a variety of tunnels ranging in test section diameter up to 5 meters and in Mach numbers 
up to Mach 3.52 Systems at AEDC employ what is probably the most advanced signal 
processing system currently in use. This processor, which has been called a "real time 
Fourier transform signal processor digitally records the entire Doppler signal, analyzes 
Its spectral content and stores this information at rates in excess of 100 KHz. This is 
essentially a specialized digital recording system that has eliminated the problem of 
limited memory described above. 

U.S. Air Force, WPAFB (USA) : A two-component backscatter LDV system has been used 
in a number of wind tunnels up to 1 meter test sections at Mach numbers up to 6.53 *$<♦ 

1SL (Germany/France) : Two-component forward and backscatter LDV systems have been 
developed for use in a variety of tunnels. A different approach to the problem of 
obtaining three wavelengths for three-component systems has been developed by Pfeifer 
(1985). In this system, a 10 watt argon laser is used to provide two of the wavelengths, 
and also to pump a 600 nm line from an R6G dye laser to provide the third wavelength. 
Also, a unique traversing corner reflector system makes traversing faster and easier than 
traversing the whole LDV system.55-60 Most large LDV systems require some adjustment in 
the output beams to insure that they cross and focus at the same point in space. A 
recently deployed ISL system employs a segmented output lens, with adjustable segments 
that provide a convenient method for making this adjustment. 

A.1.5.2 Shock Tubes - Shock tubes present a difficult applications environment because 
of the high velocities and short run times. If the velocity is required at more than one 
point, then more than one sample volume must be produced. Such facilities are usually 
high pressure facilities requiring remote operation of the system. 

A system for shock tubes and outdoor shock phenomena, developed by Spectron 
Development Laboratories 33 , has been used in the 3 meter and 7 meter CERF shock tubes at 
Kirtland Air Force Base and on the White Sands Test Range. This forward scatter system 
employs fiber-optics to conduct the light from a A-watt argon laser to four sample 
volumes that are located between two blast wings. The collected light is conducted by 
fiber-optics to the photomultiplier tubes, which are located up to 300 meters away from 
the sample volume. In such transient applications, one cannot afford to reject any of 
the data. The entire signal from the PM tubes is recorded on transient recorders, and 
signal is gleaned from noise using the procedures described in Section A.I.A.6. 

Sraeets 23 ” 30 has applied an LDV using direct Doppler frequency measurement with the 
Michelaon spectrometer (described in Section A.1.2.1) in a variety of flow situations, 
including shock tubes. Figure A.26a and b illustrates the setup. The transmitting and 
receiving optics is extremely simple for this type of system, augmenting Its use in 
severe environments. The velocity components under observation can be discerned by 
examining the directions of the virtual fringes. The fringes are produced by the 
transmitted laser beam and a second beam that is projected toward the receiver. So, for 
the case illustrated, the fringes for the s velocity component are perpendicular to the 
flow axis. 

The innate high-frequency response is illustrated in the data samples shown which 
also illustrates that the system is most readily used to measure velocity change, making 
it an excellent candidate for turbulence measurement. 

A.1.5.3 Long Range Systems - LDV systems have been used at ranges up to several hundreds 
of meters in wind sensing. Some of the earliest work in this application employed argon 
lasers 6l “*5 and CO 2 lasers.53 The basic configuration for a long-range fringe system is 
illustrated in Figure 4.27. A Cassegrainean telescope is used to project and cross two 
laser beams. The backscattered radiation is collected either by the same telescope or by 
a second receiver system that sits off-axis from the transmitter. Earlier systems 
employing standard signal processing methods required large SCs to work at long range. 
More recently, the successful application of correlators has significantly improved the 
applIcatIon.® 6 » 67 Durst et al5 4> 5S using a rotating grating beam splitter, a separate 
transmitter and receiver, a 7-watt argon laser, and a correlator in the signal processor, 
have established that such a system is applicable to a range of about 500 meters. 

Atmospheric perturbation of the beams is quite serious at visible wavelengths at 
long range. This provides incentive to go to longer wavelengths. Also OOy and YAG 
lasers are available at much higher powers. Perhaps of even greater practical Importance 
is that the 10.6 micrometer wavelength of COo eyesafe at much higher powers. Such 
systems are operated in the direct Doppler mode and measure the velocity component along 
the projected beam. The systems, having been operated at ranges up to several kilom¬ 
eters, hold great promise for use in aircraft for the detection of windahear and clear 
air turbulence. 
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Figure 4.26a Refractive index variation caused by a shock and detonation 
wave in an explorive alcohol/air mixture. 

4.1.5.4 Fiber-Optical Systems - LDV systems employing fiber-optics are rapidly gaining 
in use. An optical transmitter and receiver can be located near a measurement volume, 
while the laser and electronics are located many meters away, without much penalty in 
laser power. Equally important is the use of fiber-optics to miniaturize part of the LV 
system. 

Two systems that illustrate these features are presented. The first is the so- 
called LDV endoscope of Reithmueller 67 , illustrated in Figure 4.28. The light from a 2- 
watt argon laser is focused onto the end of a 100 micron diameter multimode fiber that 
conducts it to the measurement head. At the measurement head, the light diverges to a 
collimating lens that produces a 4 mm diameter beam. This beam Is then split into two 
parallel beams in the usual way (Figure 4.13). A focusing lens then focuses and crosses 
the two beams. The backscattered light Is collected by the same lens, which then focuses 
it onto a receiving fiber. The fiber then conducts it to a photodetector. 

In this system design, fiber-optics is exploited principally to produce a small, 
lightweight measuring head that can be traversed independent of the laser and the 
photodetectors. These elements, which are usually not only large and bulky but also 
subject to misalignment and electrical noise sensitivity, can then be located in a well- 
controlled room that is located in a convenient place. The Von Karmen Institute has at 
least four of these systems in routine operation in a variety of wind tunnels and other 
test facilities. 
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Figure 4.26b Details of optical head for aerodynamic measurement. 

Durst has even further exploited the properties of fiber-optics that allow the 
miniaturization of a laser transmitting and receiving optical head. 68 Figure 4.29 
illustrates the design concepts employed. Figure 4.29a Is a one-dimensional system, in 
which the laser beam is first split into two parallel beams that are frequency shifted by 
Bragg cells, then focused into two separate fibers. The beams are then transmitted to a 
miniaturized element that recollimates the two beams before finally focuslnR and crossing 
them at the sample volume. The backscattered light is collected by the same collimator, 
which focuses the light onto a return fiber that ultimately carries it to the photo¬ 
detector. 

Figure 4.29b Illustrates how the same concepts can be used to produce a two 
component miniaturized system. This system uses an argon laser to provide the two colors 
for the two components. These are collected and carried to the photodetector by a single 
fiber. At the photodetector, the light is split into the two components by a dichrolc 
splitter. These systems have been constructed in almost pencil-sized optical heads, 
small enough to be inserted into the head of an internal combustion engine. 

4.5.5.5 Turbomachinery and Internal Combustion Engines - The ability to produce ruRged 
LDV systems has expanded their use In tKe harsh environments presented by Internal 
combustion applications. One example is provided by the fiber-optic LDV system of 
Durst 68 described in the last section. A major problem for this application has always 
been the problem of providing optical access to the inside of the engines. This LDV is 
sufficiently small that it could be inserted into a small spark plug port. Swirl 
velocities were measured at several points in the cylinder, ranging up to about 10 meters 
per second. 
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Figure 4.27 A one velocity component long-range dual-scatter backscatter LDV system. 

Optical access methods also Include ualng glass or quartz cylinders and Installing 
windows Into the cylinders at various points. Extensive study of the windowing tech¬ 
niques for such engines has been made at the U.S. National Combustion Laboratory at the 
Livermore Laboratories.*3 Sehock 79 has developed a unique approach to the problem of 
producing a good sample volume through the otherwise aberrating cylindrical walls of a 
quartz cylinder. Using a holographic optical element designed to focus and cross two 
laser beams through the cylindrical walls without aberration, he was able to remove the 
astigmatism In the beams, producing a diffraction limited focus. The principle Is 
Illustrated In Figure 4.30. The HOE Is produced by forming a hologram of a pair of laser 
beams that have diverged from the point through the cylindrical walls. When the hologram 
is Illuminated again, the two beams have the necessary negative astigmatism to exactly 
correct that Introduced by the walls as the two beams pass through. 

Investigators at Imperial College In England have made extensive measurements In an 
Internal combustion engine model that was constructed of plexiglass to allow optical 
access. 7 ! 

LDV systems have also been used extensively In turbomachinery. 77 * 77 

4.1.5.6 Spatial Correlation Measurements - One of the most useful measurements made by 
hot-wire systems Is that of spatial correlation to characterize turbulence. To produce 
the equivalent measurement with an LDV requires that two systems be assembled In close 
and variable proximity. The complexity of accomplishing this has discouraged it for 
years, and such systems are relatively rare. Two somewhat unique systems for performing 
spatial correlation measurements are described here. 

In an extremely efficient system, Phelfer" and Schafer 79 have used the concept 
illustrated In Figure 4.31. A spherical mirror was used to collect the transmitted, 
unscattered light from the sampled volume, and then to retmage It back Into the flow at a 
slightly displaced location. In this way, two sample volumes with adjustable spacing 
were produced with only one transmitter. The two receivers are Chen placed on either 
side of the flow. 

A more elaborate, two-component, two-point LDV produced at ONERA Is Illustrated In 
Figure 4.32. This system employs two colors to produce the two components, and two 
independent sample volumes are positioned by computer controlled traverses. 

4.1.5.7 Flow Seeding - The subject of proper seeding of the flow has been acknowledged 
as one of the most critical problem areas In LDV. Fundamental to all applications of LDV 
Is the assumption that the flow can be characterized by the motions of finite sized 
particles, that either reside In the flow naturally or are Introduced Into it artifi¬ 
cially. Apparently, much of the published data taken with LDV Is In error because of the 
breakdown of this assumption. The problem Is considered of sufficient Importance that 
entire conferences have been devoted to the subject (see, e.g., NASA Conference Publica¬ 
tion 2393). 79 Nichols" has summarized the effects of particle dynamic effects on LDV 
measurements in the following way: 

(1) Interpretation of LDV measurements requires the capability of analytically 

examining particle dynamics effects. 

(2) Broad SC size distributions produce artificial turbulence and should be avoided. 
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Figure 4.28 LDV Endoscope of Reithmueller. The installation in a cooling tower test 
shows how the laser and photomultiplier can be located in a control room 
with the light being transmitted to the optics head through fiber-optics. 
The optics head can be easily positioned now because of its small size and 
weight. 

(3) Maximum acceptable SC size for typical AEDC LDV mean velocity measurements is less 
than 0.5 microns. 

(4) Absolute maximum acceptable SC size is probably driven by turbulence, not mean-flow 
requirements. 

(5) Correction for mean-flow SC lag may be feasible for monodisperse seed particles. 

Figure 4.33 presented in Nichol’s paper illustrates the response of various size 
particles to a shock wave, showing that even the smallest particles required a consider¬ 
able relaxation distance to attain the gas velocity after passing through the shock. 

Different approaches to understanding the seeding problem include drag computation, 
measurement of monodisperse particle velocity response to various known flow conditions, 
and response of particles to shock waves. 

The approaches to dealing with the problem experimentally include seeding with 
known particle distributions or preferably with monodisperse particles, and also 
attempting to reject data arising from particles that are too large to accurately follow 
the flow. In the early days of LDV, attempts were focused on making measurements without 
artificial seeding. Today, it is almost universally accepted that starting with a clean 
flow and artificially seeding with a proper material is preferred. 
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Figure 4.29b Two-component system. 

Although the detailed discussion of seeders is beyond the scope of this treatment, 
some of these will be mentioned. The reader is referred to the above mentioned NASA 
publication for a detailed discussion of seeders. Crosswy®* has provided a discussion 
and characterization of various seeders. Figure 4.33 summarizes some of the seeding 
methods. One of the most effective is the so-called vaporization condensation seeder 
using dioctyl phthalate (DOP). This produces nearly monodlsperse spheres of submicron 
diameters. Included in the figure are size distributions of other materials. 

Another popular seeder uses Dow Coming polystyrene spheres, which can be obtained 
in monodlsperse solutions with sizes as small as 0.3 micrometers In diameter. The 
solution is nebulized into droplets that mostly contain only one particle, and then 
passed through a dryer that removes the liquid by evaporation before passing the 
particles into the flow. This is practical where small amounts of seeding suffice. 
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Figure 4.30 The use of holographic optical elements in LDV to correct aberrations. 
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Figure 4.31 Optical arrangement for cross-correlation measurements. 
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Figure 4.32 Mechanical and optical arrangement of the laser veloclmeter for spatial 
correlations at A17 CEPr. L: argon laser; D: bean dividing system (cf photo 
n w 13); B: blue way; BR: blue receiver; vg: blue probe volume; G: green wav; 
GR: green receiver; vq: green probe volume; J: let exhaust (L, D, B. vg, BR) 
moves along the jet axis as for (G, vq, GR) allowing various distances, 
vqvb• But vg mey be located in various positions inside the jet flow. 

Ideally, an LDV should Incorporate a measure of particle size to validate the 
measurement. Hess 02 has developed such a system. The system adds particle measurement 
to the LDV measurement in the following way. Surrounding the normal LDV sample volume is 
a larger beam of different color or polarization (Figure 4.34b). When a signal is 
detected from the LDV volume, it Insures that the particle has passed through the 
approximate center of the large beam and, therefore, is exposed to a known Intensity of 
light. Measuring the intensity scattered from the large beam then provides a scattering 
cross section of the particle and, thus, its size. 

4.2 T ransit (Two-Spot) Anemometry (Veloclmetry) 

4.2.1 Theoretical Basis 10 

4.2.1.1 Concept - In simplest terms, a laser beam is split into two components that are 
focused to two separate points in space, and scatter centers are timed as they transit 
the two points (Figure 4.35). Even though this is conceptually simpler than LDV, LTA 
systems require a more complicated optical design to glean all of the advantages of the 
technique. Injecting all of the laser energy into two extremely small spots exposes 
scatter centers to the greatest possible intensity, enabling the detection of smaller 
scatter centers and producing the smallest possible sample volume. A smaller sample 
volume leads the way to better background noise rejection, and enables measurement closer 
to light scattering surfaces. As will become evident, the technique is more accurate for 
both speed and angle measurement. Even so, the method suffers other limitations and it 
cannot always replace LDV. 

Judicious choices of detection and signal-processing methods make possible further 
enhancement of the signal and reduction of noise. The method is directly adaptable to 
the use of photon limited detection and represents a limit in small SC detection. So, 
ultimately, much more Is involved than simply timing the transit of SC across the two 
points. Recently, the concept has been Implemented with laser diodes. 

The concept of velocity measurement, with a single-point system to produce even 
greater intensity and smaller sample volume, has been explored without much success. 83,8 ‘* 
Some of the signal enhancement methods that make two-spot systems work so well cannot be 
employed with a single spot. 

4.2.1.2 Sample Volume - The geometrical sample volume Is determined by both the 
transmitter and the receiver. The transmitter produces two focal regions of Gaussian 
Intensity profile normal to the beam propagation direction and approximately cylindrical 
in shape. The length of the sample volume In the propagation direction is determined by 
the receiver. This is illustrated in Figure 4.36. The receiver image overlaps the 
transmitted focal volume in a way that is determined by the stops used in the receiver. 
This will be discussed in greater detail below. 

4.2.1.3 Nature of the Signal - The greatest S/N can be attained if the intensities 
scattered ~Trom the two spots” - are monitored separately, and this Is the case to be 
discussed here. With ample illumination, the scattered light signal is Gaussian so that 
its peak denotes when the SC passed the center of the sample volume. There are three 
possible events: 1) the SC is detected as it passes through the first spot only, 2) the 
SC is detected as it passes through the second spot only, and 3) the SC is detected as it 
passes through both spots. 

If the signals are observed with an oscilloscope, with the first detector trig¬ 
gering the scope for each detection of an SC and the second detector voltage acting as 
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Figure 4.33a Seeding particles - use In LV. 
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Figure 4.33e Fluidized bed eeeder. 
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Figure 4.33f Alumina Polishing Power 

Nominal size -1.0 micron 
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Figure 4.33g Alumina PolishinR Powder 
Nominal Size - 0.3 micron 
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Figure 4.33h 
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Figure 4.33i Cigarette Smoke. 
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Figure 4.331 Incense smoke 
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Figure 4.33k DOP in the Laskin nozzle. Figure 4.331 DOP in the Colliaon 

nebulizer. 



Figure 4.33m DOP without dilution flow Figure 33n DOP with dilution nitrogen 
in the vaporization/con- In the vaporiation/ 

den8ation seeder. condensation feeder. 
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Fl«ure 33o Ethanol and DOP in the ratio 99/1 
in the vaporization seeder. 
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Transmitced Beam 



Figure 4.36 Illumination and receiving volume. 

the vertical input to the scope, a grouping of signals will be observed at the time S/V, 
where S is the spot separation and V is the mean flow velocity. The group will be narrow 
for laminar and broad for turbulent flow. The signal will peak when the two spots are 
lined up with the flow. Other signals resulting from case 1) or 2) above are spread 
randomly in time and, therefore, do not add up like those for case 3). 

A photon limited signal is made of random spikes that do not exhibit a Gaussian 
profile. The probability is greatest that photons are collected when the SC passes the 
peak intensity, but they can occur anywhere in the sample volume. 

4.2.1.4 Enhancing the Signal with Filters and Correlators - Since the scattered 
intensity Is Gaussian, it is optimum to electronically filter it with a Gaussian 
filter. The filter must be chosen with a passband that is compatible with the fastest 
signal. 

If the signals from the two detectors are correlated, the randomly timed signals 
representing cases 1) and 2) above will exhibit little or no correlation, while those 
from 3) will correlate strongly at time S/V and will be greatly enhanced. Consequently, 
a signal that may not be observable on an oscilloscope may be observable after 
correlation. 

4.2.1.5 Directional Sensing and Turbulence Measurement - The flow direction can be 
determined by rotating the two spots until the correlated signal shows a peak at the time 
S/V. If the flow is laminar, this can be determined precisely. If the flow direction is 
not constant, the peak will appear much broader. If the flow direction is varying 
rapidly over wide angles, then the number of occurrences of cases 1) and 2) above 
increases and the number of occurrences of case 3) decreases, resulting in a poor 
signal. Therefore, the method is not easily usable for highly turbulent flow. 

4.2.2 Optical System 

A typical optical system is illustrated in Figure 4.37. In the following sections, 
we shall discuss the selection of optical components and the various options that are 
possible. The design objective is to concentrate as much of the laser light as possible 
into two smallest possible spots, to collect scattered light from these two spots 
Independently, to reject light from everywhere else, and to provide for rotation of the 
spots. 

4.2.2.1 Beam Splitting - Several types of amplitude division beam splitters are in use, 
the two most suitable for LTA being the Rochon and Wollaston polarizing beam splitters. 
Either component produces two collimated beams of crossed polarization emerging with an 
angle between them. Other beam splitters are in use and offer advantages of adjusta¬ 
bility, but this can cause mechanical inconvenience or unreliability. 

The collimated beam from the laser is divided into two collimated beams that may be 
1 brought to two focal spots in the plane of a diffraction limited lens. This lens 

[ performs a Fourier transform function, the splitting face of the prism being the sane 

optical distance behind the back principal plane of the lens as Is the focus of the two 
spots in front of the front principal plane. In fact, this distance should be slightly 
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Figure 4.37 Significant optical components. 

modified to produce parallel focal zones in the test space. Converging beams at the 
first focus permit the relay lenses, whose longitudinal magnification would otherwise 
create diverging focal regions in the test space, to render these parallel. 

The imaged spots are relayed Into the test space by back-to-back diffraction 
limited collimator lenses. The spot size is determined by the diffraction limited image 
of the Gaussian input laser beam after magnification by the three collimator lenses. The 
separation Is controlled by the divergence of the splitter prism, the focal length of the 
focusing lens and the magnification of the two output lenses. The spot diameter to 
separation ratio Is, thus, easily designed to be an optimum value. 

If the spot separation Is 8, the beam splitter divergence is e , the Input laser 
beam diameter is D, the half width at the Gaussian waist is r, the focusing lens and 
Internal and external lenses of focal lengths f1, f2, and f3 respectively, we may write 

s - (4a) 

and 1 3 2 

r - 2/vUff/D) (f 3 /f 2 )/w (4b) 

giving a spot separation to diameter ratio of 

s/2r - wfiD/4x (4c) 

4.2.2.2 Stops - The accurate design of stops Is extremely critical for LTA systems. 
Their function is to reject light from all sources except that scattered from moving 
particles In the sampling volume. Although a simple pinhole at the conjugate image of 
each sampling region would seem adequate, several potential sources of flare must be 
considered. Multiple reflection from within the optical system, even with well-coated 
optics, are a source of scattered light. These may be reduced by double stopping where 
the image stops are themselves reimaged on two more stops. This is effective in reducing 
second and higher order scatter sources, and has advantages not only in improving the 
rejection of flare generated from within the optical system, but also from multiple 
scatter in any test situation and beam walkoff sources in windows. Typical rejection 
factors of twelve orders of magnitude have been achieved. 

In addition to the field stops discussed so far, two other types of stops are 
necessary to the system. An aperture stop is automatically provided by the outer edge of 
the final lens and by the edge of the turning mirror that separates outgoing from 
Incoming light. It is necessary to extend the stopped area of the turning mirror to 
prevent diffracted flare from its edges. This is conveniently done by two more stops, 
one of which restricts the input beam to an area slightly smaller than the turning 
mirror, and one of which occludes an area slightly larger than the turning mirror from 
•r behind. 



Flare may be seen by one photomultiplier that has arisen from light in the other 
channel. This may be prevented by stopping part of the received area. This results in a 
reduction in optical efficiency of the machine, and this stop is only necessary when wall 
proximity is important. This is executed as an extension to the stop immediately 
adjacent to the turning mirror. 

4.2.2.3 Spot Intensity and Aberrations - To achieve the maximum benefit of the LTA 
principle requires that the system be nearly diffraction limited. Because of the small 
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field, the primary aberration la spherical. Projected spot quality and receiver 
efficiency are both Impaired by normal and more by oblique incidence through a thick 
window. Typically, up to 5 cm thick windows with less than 10 degree Incidence do not 
cause serious deterloration with f/4 systems. Curved windows, like large angles of 
Incidence, Introduce astigmatism. 

The photon flux through the sampling volume is 

2 

- P/hv*r photons/sec/area (4d) 
where P Is the laser power, r Is the radius, and hv is the energy per optical photon. 

The photon flux, scattered by a particle of radius, a, during Its transit of the 
peak, is F q(a) »a*F,q(a) . where q Is the scattering efficiency. Hence, the total 
number of photons scattered during the transit, G e , will be the flux mutlplled by the 
transit time, 2r/V, hence 

G g - 2, « 2 Pq/hvrV (4e) 

This Is necessarily approximate, because no allowance is made for which part of the 
Illuminated region Is intercepted by the particle, nor the profile function, nor the 
collection and detection efficiency for scattered photons. 

4.2.2.4 Spot Rotation - Rotating the orientation of the two soots is required to 
determine flow direction. This can be done by 1) rotating the beam splitter and the 
8tops, 2) using an image rotator prior to the output lens, and 3) rotating the whole 
machine. 

The last method Is simple but mechanically difficult. The first two lead to a more 
compact design. The first has more moving parts with opportunity for less than perfect 
alignment. The second Is efficient and rather expensive. It Is usual to arrange for 
computer control of rotation with corrective feedback of the angular position attained. 
Rotational precisions of 0.1 degree are easily attainable. 

4.2.2.5 Detection - The detection of a transit event Is straightforward, but, as the 
system becomes photon limited, other factors come Into play. In an experiment without 
additional background light, where only light scattered by Isolated particles Is 
detected, two to four photons per spot are required for an LTA. This Could actually be 
reduced to one per spot If the spots are far enough apart and there are no other photon 
sources. 

In roost experiments, the peak photon rate is almost always too high for quantum 
resolution, and analog filtering represents an excellent method of discriminating events 
from nonevents, where an event is associated with a group of photons from a single 
particle. A separate filter, that makes an estimate of time of occurrence of the highest 
Instantaneous photon flux, may generate a digital signal to be correlated with a similar 
signal from the second spot. 

4.2.3 Signal Processing and Interpretation 85 , 88 

4.2.3.1 Broadening - Broadening of the data from an LTA is a manifestation of the 
Inability to know precisely when the SC crossed the centers of the two spots. Broadening 
decreases as the spot diameter to separation decreases. Signal-to-noise ratio Is 
critical here. When a photon limited signal occurs, a positional error for any one event 
can be almost as high as the spot diameter. 

4.2.3.2 Data Rate - Data rate is proportional to the swept volume. This is the cross 
sectional area normal to the flow multiplied by the velocity. Natural aerosol distribu¬ 
tions exhibit a number density approximately proportional to the inverse fourth power of 
radius, so. In an unseeded flow, the data rate can be vastly Increased by extending the 
detection range to smaller SCs. 

4.2.3.3 Time Measurement - Direct timing and correlation are used to measure the time 
intervals^ In sparse particle conditions, where the transit time from one spot to the 
other Is short compared with the mean interval between particles, there appears at first 
to be no significant advantage of one method over the other. Where particles occur more 
frequently, the efficiency of timing falls rapidly compared with correlation. In all 
cases, it is much more satisfactory to use two photomultipliers. With one PMT, timing 
can only stop on the next consecutive event, and this yields a background that must be 
compensated and Is difficult to model. Autocorrelation yields a flat and subtractable 
background that is four times greater than that derived by cross-correlation of the 
signals from two photomultipliers. 

4.2.3.4 Data Analysis - The signal from the two photodetectors are cross-correlated to 
produce a transit time probability curve. Because the background level is constant, it 
can be subtracted from the curve. The resulting transit time probability histogram must 
be converted to a velocity probability distribution. A velocity bias must be removed in 
the conversion. As with all LV systems, the effects of SC bias must be considered. The 
two-spot system measures magnitude and angle over an angular range governed by the 
spacing and diameters of the spots, typically the order of one degree. A curve of the 
type illustrated in Figure 4.38 is generated for each angular setting of the two spots. 
For a full angular sweep, a curve such as that shown In Figure 4.39 results. 
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Figure 4.38 (a) Transit tine correlograms, (b) Velocity probability curve 

from (a) and (c) Bias Introduced by multichannel analysis. 



Figure 4.39 Velocity/probability plot from 
laser anemometer. 


4.2.4 Comparison of LTA to LDV 

Applying the equations presented above and those developed earlier for LDV, It is 
straightforward to compare the properties of the two system types. For particles of 
radlua, a. and sample volume radii of r D and r» for an LDV and an LTA respectively, 
assume that the same laser power is used for eacn system. Some comparative parameters 
for the systems are shown In Table 4.1.96 


Table 4.1 

RATIOS OF LTA VERSUS LDV FOR VARIOUS PROPERTIES 
OF EQUIVALENT SYSTEMS 


Ratio LTA/LDV 


Number of photons scattered In transit 


Peak photon rate detected 




2 


Sample volume length 


Data rate for monodlsperslon 


Date rate for n(a) a 

(aj,,a T - minimum detectable radius) 



Data rate for N(a) (with allowance 
q(a) and retrieval efficiency) 


k 3/4 r T 
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5/4 
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A typical ratio of probe voluae radii for high-speed flow r D /r T la ten to twenty. 
Using this with the quantities In the above table, a number of generalisations can be 
made. 

(1) The nuaiber of photons scattered for a given SC Is larger by over an order of 
aagnltude for the LTA. Therefore, the LTA can detect saaller particles. This 
tends to Increase the data rate, although not much since the scattering cross 
section Increases as the fourth power of the SC radius. 

(2) The saaple voluae for an LDV la such larger than for an LTA. This can make the 
data rate for the LDV higher than that for the LTA unless the LTA is able to detect 
many aore saall particles than the LDV. When the flow la seeded with particles that 
can be detected by both systems, then the data rate for the LDV will be higher. 
The smaller saaple volume length of the LTA allows measurements to be made closer 
to a surface than with an LDV. Typical ratio of probe volume radii for high-speed 
flow r])/r T Is 10 to 20. 

(3) The LTA Is not easily used In highly turbulent or vortical flow. 

(4) The LTA can be used In higher speed flows, since the spot separation can be 
adjusted to Improve the accuracy of the transit time measurement. For the same 
reason, the speed measurement can be made more accurately. 

(5) In nearly laminar flow, the flow direction can be measured more accurately with the 
LTA. 

4.2.3 Applications 

LTAs have been applied successfully in a variety of flow fields. In the following 
sections, some examples are given to Illustrate the features of the technique. LTAs are 
commercially available from at least two companies. 

4.2.5.1 Turbomachinery - The mapping of flow velocity between the compressor blades of 
turbomachines 1 0 •ss- >s has been an Important application of LTA. This Is an extremelv 
demanding application, because the measurement must be made In a hackscatter mode, 
through small windows that quickly become dirty, close to both the windows and the 
blades, and with little or no seeding of the flow. All of these are conditions In which 
the LTA excels. 

Scientists at Rolls-Royce 77 pioneered this application. Their system, Illustrated 
In Figure 4.40, employs a two-watt argon laser focusing Its energy Into two 16-mlcroraeter 
diameter spots separated by 160 micrometers. The optical assembly Is contained In a box 
that Is traversed by a three-dimensional traversing system. The measurements In rotating 
blade rows are made by gating one of the two photomultiplier tubes, so that It accepts 
data only when the sampling volume Is In one selected part of the blade row. The gate Is 
typically 5 or 10 microseconds long and Is at a known time from a triggering signal 
produced by a blade. Alternately, It can be set to occur on every passage to obtain the 
average velocity at the same position In all passages. Typical velocities measured In 
this application range from a few m/sec to hundreds of m/sec. 

Detailed velocity maps have been produced between the blades, resulting In the type 
of data illustrated In Figure 4.41. 

4.2.5.2 Wind Tunnels - LTA Is useful In wind tunnels, especially where only small SC are 
present and where the scattered light Is too small for LDV. Also, in applications where 
measurement must be made close to models or windows, LTA excels. LTA Is an excellent 
choice for measuring flow angularity. Even so, LDV has largely dominated wind tunnel 
applications because of versatility and ability to characterize turbulence. 

4.2.6 Two-Line Anemometers 

To overcome one of the principal limitations of LTA, researchers have developed 
systems that employ two lines Instead of two spots." 5 The optics for such a system are 
Illustrated In Figure 4.42. The unique properties of several optical components provide 
an efficient system. The two major lines of an argon laser are used to produce the two 

lines. The two lines are produced by separating the two colors angularly with a prism. 

Converting the two spots into two lines can be done by introducing astigmatism into the 
optical system. In the Illustration, this Is accomplished by passing the diverging light 
through tilted plates. The lines In this system are 20 by 100 micrometers In 

dimension. The collected, scattered light is separated into Its two-color components by 

a dtchrolc reflector. 

The two lines can also be generated by using cylindrical optics In the projecting 
system. Although somewhat more complicated, this produces a better sample volume. The 
astigmatic focus Is complicated by the existence of a circle of least confusion, and a 
second focal region that contains a pair of lines oriented in a direction normal to the 
first. 

The two-line LTA offers some of the advantages of both LDV and LTA. The system can 
handle flows of greater turbulence, but the benefits of rotating the two spots are lost. 
These kinds of systems are relatively new and their full benefits have not been widely 
demonstrated. 
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Figure 4.41 Laser anemometer results at 701 blade height 
in a research fan. 
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CHAPTER 5 

LASER SPECTROSCOPY AND SCATTERING METHODS 
5.0 INTRODUCTION AND BACKGROUND 

Light (photon) interaction with molecular matter is classified as either "elastic" 
(wherein the photons do not lose or gain energy) or "inelastic" (wherein the photons do 
exchange energy with the internal energy states of the matter). The classifications are 
illustrated ir> Figure 5.1. 


INTERACTION BETWEEN LIGHT AND MATTER 


ELASTIC INELASTIC 

LIGHT IS SCATTERED INSTANTANEOUSLY 
WITHOUT EXCHANGING ENERGY WITH THE INTERNAL 
STATES OF THE GAS MOLECULES 


RAYLEIGH SCATTERING HIE SCATTERING 


MODULATION OF THE INCIDENT LASER BEAM 
BY THE INTERNAL ENERGY STRUCTURE 
OF SPECIFIC MOLECULES 
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a. Overall Classification 
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b. Raman Interaction 


RAMAN ACRONYMS 

COORS Conventional Old-Fashioned Ordinary 

Raman Scattering 

CARS Coherent Anti-Stokes Raman Scattering 

IRS Inverse Raman Scattering 

CSRS Coherent Stokes Raman Spectroscopy 

(pronounced scissors) 

RIKES Raman Induced Kerr Effect Spectroscopy 

HORSES Higher Order Stokes Effect Scattering 

PARS Photo-Acoustic Raman Scattering 

CARE Coherent Anti-Stokes Raman Elllpsoaetry 
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ENERGY LEVEL DIAGRAMS 



Emission Absorption 



Fluorescence Raman Stokes 



Raman Anti-Stokes CARS 


c. Energy Level Diagrams for the Various Inelastic Processes. 


Figure 5.1. Continued 

The elastic scattering process is essentially instantaneous. One can think of the 
photon striking the molecule and rebounding with no energy loss like a hard sphere 
collision. Electromagnetic theory models the interaction as the acceleration of the 
bound electrons in the molecules by the changing electric field of the impinging 
radiation and the electrons radiate at the same frequency as the field which causes their 
acceleration. Even though no energy exchange occurs, the interaction can be used for 
diagnostics because changes are effected in the distribution of the radiation in space as 
well as the polarization state of the radiation. Also, if the molecule is in motion, the 
frequency of the radiation is Doppler shifted by the scattering center by an amount that 
is directly proportional to its velocity. 

Each molecule has a scattering cross section that determines the amount of light 
scattered from impinging radiation as well as how the scattered radiation is distributed 
in space. The cross section depends on the electrical properties and the size of the 
scatter centers. The total amount of radiation scattered from the Impinging wave is 
directly proportional to the number of molecules and thus, the number density. 
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The diagnostics process can be further enhanced by illuminating the volume of 
Interest with well defined distributions of light, such as one or more spots or plane 
sheets. Both transmitting the light to the region of interest as well as collecting the 
scattered light in a carefully restricted way with specially designed apertures and 
lenses can improve the signal to noise ratio and help to spatially resolve the place of 
measurement. 

As scatter centers pass from one point to another in a well defined light field, 
the time and space dependence of the scattered light can be used to determine velocity 
and turbulence, as described in Chapter 4. 

A photon can Interact with matter inelastically in the following five different 

ways (Figure 5.1c). 


1. Absorption - The photon is absorbed entirely through an energy exchange process in 
which 1 ts energy Is absorbed, forcing the molecule into a higher energy state. 
This can be an electronic energy state or a molecular rotational or vibrational 
state or a combination of all three. The change of energy of the molecule equals 
the energy of the photon. 


Subsequently, the molecule can decay spontaneously to a lower energy state while 
emitting a photon whose energy equals the change In energy of the molecule. 

Spontaneous emission occurs after about 10-* seconds for most molecules. Some 

states are much longer lived, lasting even seconds. These are called the 

fluorescent states. 

2. Stimulated Emission - The striking photon stimulates (or drives) an excited 

molecule down to a lower energy state resulting in a decrease in energy equal to 
the photon energy. This transpires three orders of magnitude faster than 

spontaneous emission. Then the striking photon and the emitted photon will have 
the same frequency and leave the molecule in phase with each other. 

3. Raman-Stokes Scattering - A photon strikes a molecule, adds part of its energy, 
and drives it to a higher level. To conserve energy, a photon must be emitted to 
carry away the remaining energy. The scattered photon has lower energy than the 
striking photon (and, therefore, lower frequency). The probability (or cross 
section) for this type of interaction is about three orders of magnitude lower 
than that for elastic scattering. 

4. Raman-Anti-Stokes Scattering - An already excited molecule is driven to a lower 
level by a striking photon causing an energy decrease in the molecule that (unlike 
the stimulated emission described above) is not related to the energy of the 
striking photon. To conserve energy the emitted photon then has greater energy 
(higher frequency). The cross section for this interaction is even less than that 
for the Stokes interaction. 

5. Stimulated Raman Scattering - The photon strikes a molecule that has been excited 
to a virtual state and stimulates (drives) it to a lower state with an energy 
decrease equal to the energy of the photon. This can take place only in 
conjunction with another photon that is undergoing Raman scattering. The photons 
that leave the interaction have a distinct phase relationship with respect to each 
other. 


f 





This chapter covers: 

1. elastic scattering from molecules (Rayleigh Scattering) to measure density, 
temperature and velocity, and 

2. inelastic scattering (absorption/emission, fluorescence, and Raman spectroscopy) 
to measure temperature, species, velocity, and density. 

5.1 Rayleigh Scattering 

The measurable parameters using Rayleigh scattering are temperature, density,and 
velocity. 

5.1.1 Principle of Operation 

Figure 5.2 Illustrates a typical Rayleigh scattering instrument for the 
measurement of the above parameters. A pulsed laser is focused to produce a high 
intensity illumination of a small volume in the gas. A large aperture, high resolution 
imaging system further restricts the collection of scattered light to a small volume in 
space. The total amount of scattered light is proportional to the number of molecules 
illuminated and, therefore, is proportional to gas number density. 

The linewidth of the scattered illumination is broadened by the Doppler shift 
impressed on the scattered light by the motion of the molecules and, therefore, is 
proportional to (temperature) 1 /*. The center frequency of the scattered light is Doppler 
shifted by the bulk motion of the gas by an amount that is proportional to the gas 
velocity. 
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Figure 5.2. Rayleigh scattering diagnostics. 

The Rayleigh scattering cross section for ambient air is large compared to the 
cross section for Inelastic scattering, but it is still so small that high power lasers 
must be used to produce a measurable amount of scattered light. Consequently, a pulsed 
laser must be used. The cross section varies directly with the minus fourth power of the 
wavelength, so much gain is to be had by using short wavelength lasers. This technique 
has flourished especially with the emergence of ultraviolet excimer lasers that can 
provide short wavelength light at intensities never before available. 

Perhaps the greatest attraction of all for Rayleigh scattering is the prospect of 
measuring temperature, density, and molecular velocity all simultaneously (Refs. 1 and 
2). The essential feature of Rayleigh scattering is that the scattered signal is not 
shifted in frequency with respect to the incident beam (except for a possible Doppler 
shift). 


5.1.2 Analytical Description 

The intensity of Rayleigh scattered light, I, from an incident beam of Intensity, 
I, can be written as: 


1 r - I L oNta (5.1) 

where N is the gas number density, n, the solid angle of collection, t the probe volume 
length in the laser beam direction, and <j the total scattering cross section. If the gas 
constituents are known, the total scattering cross section may be calculated from the 
scattering cross section of the individual constituents such that 
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H a i 


(5.2) 


where x^ is the mass fraction of constituent i and oi is its scattering cross section. 
Individual Rayleigh scattering cross sections can be accurately predicted using Mie 
theory, and are functions of the wavelength, complex index of refraction and the 
depolarization factor. 
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Equations 5.1 and 5.2 can be used to determine the number density from the 
measured Id and the temperature can be deduced from the spectrum of the scattered 
light. The spectral distribution is formed by the Doppler shifts of the light scattered 
from individual molecules in thermal motion. Namely, 

I R L 

I x - - 77 - exp (- '. 2 * 2 ) (5.3) 


where 

L2 - M c2 / 8 kT \2 sin 2 (e/2) (5.4) 

and c is the speed of light, k is Boltzmann's constant, M is the molecular weight, T the 
translational temperature and « - x - X 0 is the shift from the incident wavelength. 
Take, for example, nitrogen molecules, a laser wavelength of 488 nm, and backscattering 
(8 - 180 degrees), then the full width half maximum will be related to the temperature 
through 

FWHM - 0.0235 /T/1000(°K) A (5.5) 

The ability of the technique to measure temperature depends on the laser line 
width. Typical laser linewidth at 300 nm is of the order of 0.001 A, which is far below 
the FWHM calculated from Equation 5.5. 

5.2 Absorption/Emission Spectroscopy 

The measurable parameters using c’^sorption/emission spectroscopy are temperature, 
pressure, species concentration and vedrxty. 

5.2.1 Principle of Operation 

Emission spectroscopy analyzes the spontaneous radiation spectrum from a sample 
volume as compared to absorption spectroscopy which measures the attenuation as a 
function of wavelength of radiation that is passed through the medium from a light probe 
(Figure 5.3). The temperature is inferred from the shape of tft« absorption or emission 
profile. The concentration or pressure is inferred from the amount of emission or 
absorption, and the species concentration is inferred from the overall absorption or 
emission spectrum. 

TUNABLE 



DETECTOR 

ARRAY 

Figure 5.3. Species identification and temperature measurement by absorption 
spectroscopy 

In the case of absorption, the source may be either continuous, discrete, or 
ultra-narrow bandwidth in the form of a laser. This discussion will concentrate on 
absorption spectroscopy using laser light sources. With a laser source the spectral 
purity can be extremely high and the spectral linewidth can be less than the individual 
lines in the absorption spectrum. Therefore, one laser line can be tuned over an 
absorption line to measure the broadening of that line to infer temperature. 
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The measurement is an integrated effect and like interferometry, the point 
function must be determined by solving the integral. If the distribution of the measured 
parameter is constant, the solution is simple. If the distribution is axisymmetric the 
solution is the so-called Abel inversion. If the distribution is three-dimensional, then 
a tomographic analysis such as described in Chapter 2 is required to produce the point 
measurement. 


5.2.2 Analytical Description 


The intensity of a beam of light passing through a medium along s changes in 
accordance with the following equation: 


di /ds ® -k (i 
v v v 

is the intensity in watts/steradian/cm*, 


.) 


(5.6) 


where t v is the intensity in watts/steradian/cm*, B is the blackbody function, and k is 
the absorption coefficient. The solution to this equation contains two terms, v one 
describing the absorption of the incident radiation and the second describing the 
radiation added along the path from molecular emission. 
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(5.7) 


The second term depends on the blackbody function times the absorption coefficient 
integrated along the pathlength attenuated by the intervening gas between the observer 
and the emitting volume. The absorption coefficient for a single transition, k in units 
of inverse length is given by 


k 

v 



4 1 f(v) 


(5.8) 


where A is the Einstein A coefficient, f the normalized lineshape function, N the state 
number density, and g the degeneracy of the energy level. 

The Einstein coefficient indicates the absolute probability of a transition from 
one energy state to another and is determined by the quantum mechanics of the two states. 

The lineshape function depends upon several mechanisms. For relatively low 
pressures, Doppler broadening dominates and the shape parameter is a Gaussian with 
temperature as the primary variable. If pressure is high enough, the broadening 
mechanism becomes Lorentzian in shape with the collision frequency the dominant factor. 
Collision frequency depends on both pressure and temperature. 


The Doppler width is given by 
b d = 

and the collisional width is given by 

b _ 
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where b c is the measured width at a standard temperature and pressure. 

If Doppler and collisional mechanisms are of the same order, then the Voight 
(mixed lineshape) function must be used (see, for example, Refs. 3 and 4). 


Absorption spectroscopy has not vet been widely used in flow diagnostics; however, 
the availability of usable hardware and the need for diagnostics of unseeded flows is 
likely to change this in the near future. The strong oxygen absorption band in the 170- 
190 nm region can be used to measure temperature with an excimer laser as a light 
source. The technique, in principle, is rather simple, covers a wide dynamic range and 
requires only small windows. 


5.3 Laser Induced Fluorescence 


The measurable parameters using laser induced fluorescence are number density of 
Che fluorescing species and temperature (from which pressure can be derived from the 
equation of state) and velocity. 

5.3.1 Principle of Operation 

A tunable laser is tuned to resonantly excite a particular transition of a species 
in the test volume. The species may be naturally occurring or perhaps a more convenient 
species introduced into the flow as a seeding material. Fluorescence is a spontaneous 
radiation that arises from the decay of the excited molecules. A typical experimental 
arrangement for laser induced fluorescence measurement is shown in Figure 5.4. The 
fluorescence signal can be used in the following ways: 

1. The population of the excited state, which is proportional to the amount of 
radiation can be related to the total number density of the species. 
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2. By measuring the population of two states near the ground state, the population of 
the ground state and the temperature can be determined from the Boltsmann 
distribution. 

3. The narrow linewidth laser can be tuned over the absorption line and the 
fluorescence spectrum compared with that of a reference stationary cell (see 
Reference 5). The frequency shift of the center frequency is the result of the 
Doppler shift caused by moving molecules. The broadening of the spectrum can 
provide temperature and pressure if the broadening parameters are known. The 
requirement of a reference cell can be obviated by using counter-propagating laser 
beams. This provides a Doppler upshift from one of the beams and a Doppler 
downshift from the other, and the average of the two is the equivalent of the 
stationary velocity. 



Figure 5 4 Planer laser induced fluorescence (PLIF) measurement of gas density. 

5.3.2 Analytical Description 

An approach for analyzing the fluorescence signal is described by Hiller, 
Reference 6. The laser frequency is tuned to the linear region of the fluorescence 
intensity versus wavelength curve. The fluorescence from two opposing beams is Doppler 
shifted up and down by the velocity of the gas, and the amount of the intensity change 
corresponding to the Doppler shift is (see Figure 5.5) 



s<*> 



Figure 5.5. Fluorescence intensity versus frequency for static absorption line (heavy 
line) and for blue-shifted (laser beam in same direction as velocity 
component) and red-shifted cases. The curves are drawn for one specific 
velocity u. The laser frequency is v Q . The approximations assume that S. 
and are located on the linear part of the line. (After Hanson, Ref. 7) 
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The frequency shift is 



The derivative must be determined from the calculated or measured specular distribution 
function. 

In the work reported in Reference 6, three components of velocity were measured by 
probing the gas from three directions with four laser sheets. The measurement was done 
in 250 ms. 

Recent implementations of this method employ laser sheets so that the velocity 
distribution over an entire plane can be measured. The sheet of light is imaged onto a 
photodiode array camera after having been intensified by an image intensifier. The laser 
beams from the different directions can be recorded successively with a full matrix of 
intensities being recorded for each. Then Equation 5.12 is solved for Doppler shift for 
each point in the sheet. 

5.3.3 Special Considerations 

The limiting factors in the accuracy of the velocity measurement are: 


1. Lifetime Broadening 
2# Collision Broadening 

3. Thermal Broadening 

4. Laser Broadening 

5. Turbulence Broadening 


The typical errors cited by researchers is about 5 m/s. This technique is more 
applicable to high speed flows, and has been demonstrated in subsonic jets at about 
50 m/s (Ref. 7) and for hypersonic flows up to about 1700 m/s (Ref. 8). 

To be used at low speeds, the gas pressure must also be small to produce narrow 
lines from which the Doppler shift can be extracted. The highest practical pressure is 
about one atmosphere (Ref. 9). Other limitations are caused by the imperfections or dirt 
on the test cell windows that will affect the fluorescence intensity. Spectral analysis 
of the fluorescence circumvents this problem at the expense of a longer processing time. 

5.4 Raman Spectroscopy 

The measurable parameters are species concentration, vibrational and/or rotational 
temperature of the particular species, and flow velocity. Translational temperature can 
be estimated with considerably less accuracy from the Doppler broadening of the scattered 
light. The technique is typically useful for concentrations of about 1 part per million. 

The variations of Raman spectroscopy are many (see Figure 5.1) and cannot all be 
discussed in this work. However, the most important are Conventional Raman (COORS), 
inverse Raman (IRS) and Coherent Raman Spectroscopy (CARS and CSRS) and a brief 
discussion is Included. 

5.4,1 Principle of Operation 

About one out of every thousand times a photon strikes a molecule it loses a 
fraction of its energy to the molecule and alters the internal energy states of the 
molecule, and the emerging photon has a longer wavelength (typically tens of angstroms) 
from the original photon. This is called Raman-Stokes scattering. The observed 
frequency shift depends upon which energy states were involved in the interaction, so 
what is observed is actually a spectrum of frequencies that are characteristic of the 
molecule. Every molecule has a unique Raman spectrum. This provides a fingerprint of 
the molecule and is, therefore, a diagnostics for the type of species present in the 
sample. 


How the intensities are distributed between the frequencies in the Raman spectrum 
depends on the original distribution of energy states in the sample, therefore, providing 
a measure of the temperature associated with the particular state, such as vibrational or 
rotational temperature. 

In a related but less probable interaction a photon colliding with a molecule can 
emerge from the interaction with greater energy (and higher frequency) by leaving the 
molecule in a lower energy state. The scattering process is called Raman anti-Stokes. 
The same fingerprinting characteristic applies, and the ratio of the Stokes to anti- 
Stokes scattering cross section provides a measure of the initial distribution between 
energy states and, therefore, can be used to measure temperature. 

In any of the cases, the flow velocity of the molecules can be determined through 
the Doppler shift of the scattered radiation. A calibration cell Is used to provide the 
spectrum with no Doppler shift and that is compared with the spectrum of the light 
scattered from the volume under measurement. 
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5.4,2 Technique Description 

Raman spectroscopy is relatively simple experimentally. A laser beam is 
concentrated in a region of interest and the scattered light is collected and analyzed 
(see Figure 5.6). The largest amount of light emerging from the sample volume arises 
from Rayleigh or Mie scattering, or fluorescence. However, the Raman scattering occurs 
at a different wavelength and can be separated from the composite by using narrow-band 
filters, or a monochromator. Then the intensity of the light is measured wth a 
photodetector or recorded on film. 



Figure 5.6. Species identification by Raman Scattering. 

The complete spectral analysis of the signal requires a monochromator or 
spectrograph, but in special cases the measurement is much simpler. For example, if the 
constituents are known, then their concentration can be measured simply by measuring the 
amount of light existing in certain narrow spectral bands. Alternately, the distribution 
of a species in space can be observed by illuminating the volume with a laser sheet, then 
recording through a narrow line filter centered about one of the Raman spectral 
components of the species. This is called a Raraano graph. The signal Is extremely weak, 
so a combination of image intensifiers, multiple exposures, and high power lasers is 
required. 

When the laser line approaches resonance with a state in one of the molecular 
species, Raman scattering is enhanced. The transition and, indeed, the distinction 
between resonance Raman and resonance fluorescence is not completely understood, but 
Raman scattering is defined as that process that is almost instantaneous, involving 
virtual states while fluorescence is a relatively time-consuming process (typically 
microseconds) involving real states. 

The structure of the Raman spectrum can be extremely complex, even for simple 
molecules, because all of the allowable transitions in the molecule can become 
involved. The impinging photon can drive the molecule into any of the allowed final 
states as long as energy is conserved, (unlike the process of resonance absorption). In 
quantum mechanics, the allowable changes in quantum number in an energy transition are 
zero and plus or minus one for vibration and zero or plus or minus two for rotation. So 
for a given vibrational state transition, for example, a change of vibrational quantum 
number zero to one, the spectrum contains three branches. The excitation to higher 
rotational quantum numbers produces what is called the S branch, while excitation to the 
lower quantum numbers produces the 0 branch, and no change in rotational number produces 
the Q branch. 
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5.5 Coherent Raman Spectroscopy 

The measurable parameters using CARS, CSRS, and IRS are the same as with Raman 
spectroscopy, that Is, species, and rotational/vibrattonal temperature and flow velocity. 

5.5.1 Principle of Operation 

The difficulty in using conventional Raman spectroscopy arises from the small 
Raman cross section and also from the homogeneous scattering in all directions, making 
the collection of scattered light inefficient. Coherent Raman solves this problem by 
stimulating the process with a tunable probe laser which forces the molecular transition 
to occur with a fixed phase relationship to the pump and stimulus beams. Figure 5.7 
illustrates a CARS experiment. 



(b) Physical Process 



(a) Experimental Set Up 


Mono chromat o r 
or 

Detector 


Figure 5.7. Description of a CARS experiment 

1 . The Impinging photon from a pump laser raises the energy state of the molecule to 
a virtual level, which is hyp above the initial state. 

2. The stimulus or probe laser is tuned to stimulate a downward transition to a state 
below the virtual state by an amount hvg . 

3. A second photon strikes the molecule adding energy hv p to the molecule raising it 
to a 81ill higher virtual state. 

4. The molecule decays to the initial state emitting a photon of energy 2hv p - hv g . 

5. The emitted photon leaves in a preferred direction, defined by the directions of 

the input beams. This is a result of a constant phase relationship between the 
photons. 

6. The probe laser can be tuned over the entire Raman spectrum allowing a plot of the 
spectrum with a resolution that depends only on the linewidth of the probe laser. 

7. In CSRS the process involves the stimulated Stokes process. 

S. In IRS the energy removed from the probe beam by the stimulation process is 

measured. The process is called "inverse” since a decrease in intensity is 
observed. 
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5.6 Application Examples 

5.6.1 Rayleigh Scattering for Measurement of Density in the Shuttle Reentry Path 

Rayleigh scattering has not yet been applied routinely in aerodynamics. An 
application planned by McKlnsey (Ref. 10) for measurement of the density along the flight 
path and the flow field near the space shuttle during its reentry provides a good example 
of the potential of the technique. The planned system is illustrated in Figure 5.8. The 
recommended system will employ an ArF laser operating at the 193 nanometer wavelength. A 
pulsed output energy of less than 1 mlllljoule was shown to be sufficient to provide an 
adequate return signal. (Commercially available lasers now exceed 200 mlllljoules at 



Figure 5.8. Rayleigh scattering setup for space shuttle application. 

pulse repetition rates in excess of 200 Hz.) The Cassegrain telescope has a collecting 
aperture of 23 cm. The photomultiplier is Cs Te on fused silica having a quantum 
efficiency of 15% and a maximum gain of 2 x 106. A 22 nm narrow-band filter removes 
ambient light. 

The system will provide a Rayleigh return from oxygen and nitrogen molecules in a 
column that extends out to about one kilometer from the shuttle. The integration effect 
allowing contribution from all of the molecules in the long column is one of the keys to 
the success of the application. The cross section of air used for the computations was 
about 3.5 x 10"2* cm2. (After McKenzie. Ref. 10) 

Computations of the signal-to-noise ratio show that such a system is capable of 
providing direct measurements of ambient density with an uncertainty of less than 1% over 
an altitude range from 50 to 90 km with a spatial resolution of 1 km. 

McKlnsey has also produced a conceptual design of an instrument that could 
spatially resolve the density profile Inside the shock wave of the shuttle and could 
actually locate the shock. Such a system must time resolve the return signal and 
associate the Intensity of the return signal at a known time with a distance from the 
shuttle. This is a much more demanding task both on the laser and on the detection 
system, since one must detect scattered power and not total integrated energy as in the 
above application. Therefore, a much more powerful laser is required and a detection 
system that has a frequency response in excess of 1 GHz for a resolution of about .3 
meters. 


5.6.2 Laser Induced Fluorescence Measurements of Oxygen in Supersonic Flows 

Temperature measurements have been made in supersonic flows below 300 K using 
laser induced fluorescence on NO seeded into a nitrogen flow (Ref. 11) and iodine seeded 
into an air flow (Ref. 12). Velocity and pressure fields in supersonic nozzle flows have 
been determined using the Doppler shift of molecular absorption lines for iodine seeded 
into a nitrogen flow (Ref. 13). The use of seeding materials in flows is rarely 
desirable or easy to maintain. Unfortunately, some of the best seeding materials for 
fluorescence are highly toxic and unsafe to use. 

An extremely attractive candidate for LIF measurements would be molecular oxygen 
because of its natural presence in most flows of interest. Reference 14 describes a 
study in which oxygen was used as the fluorescing material. Available excimer lasers 
operating far in the ultraviolet make such measurements possible. The availability of 
the narrow-band laser source also makes LIF Doppler-shift velocity measurements 
possible. If a component of the velocity lies in the direction of the exciting 
radiation, the molecules absorb at a frequency that is Doppler-shifted from the laser 
frequency. Since the place at which the laser falls on the absorption lineshape function 







varies with velocity, the fluorescence intensity will be be a function of velocity 
(Ref. 15). (See Figure 5.5) 

The experimental setup is similar to that shown in Figure 5.4. The laser, tuned 
to the oxygen absorption transition at 193.055 nm, was the excitation source. This wave¬ 
length was chosen because it experiences negligible absorption at room temperature (so 
little loss is experienced in getting the light to the measurement volume). Absorption 
increases by an order of magnitude between 300 and 1000 K. The output of the laser was 
nominally 100 mj in a 13 nsec pulse at a rate of 11 Hz. A 38 cm focal length cylindrical 
lens and 100 cm focal length collimator formed the light sheet. The sheet was about 100 
micrometers thick at the sample region. The lenses were of UV-grade fused silica. 

The sensor was placed at 90 degrees to the sheet and viewed through a KBr high- 
pass filter. The sensor comprised an ITT Model 4144 dual-plate image intensifier with an 
S-20 cathode. Mounted to the image intensifier was a Reticon Model MC520 100 x 100 
element array camera using fiber-optic bundle coupling. The data from the image array is 
then fed to a computer which processes and displays the data in two-dimensional formats. 

The resulting oxygen fluorescence images demonstrated velocity, temperature, and 
pressure sensitivity. Separating the parameters velocity, temperature, and pressure is 
not a trivial process. The recommended procedure emerging from this study is to use an 
absorption line that is insensitive to temperature to measure pressure. The collection 
is done at 90 degrees where there is no Doppler shift to isolate the measurement from 
velocity dependence. To measure temperature, a line that is highly dependent on 
temperature can be selected. Finally, velocity measurement must be done at an angle 
different from 90 degrees with respect to the illumination. 

In reacting flows where the oxygen concentration changes, the measurement can be 
used to determine oxygen number density. Where multiple parameters are varying, the 
possibility of using two exciting wavelengths with broadband collection or a single 
wavelength with two detection channels was suggested. 

5.6.3 Wind Tunnel Applications of Coherent Raman Scattering 

A limited number of applications of Raman scattering in wind tunnels have been 
demonstrated to date but some of these are extremely promising. Raman Doppler 

velocimetry (RDV) is important since the particle lag problem in standard LDV is 

solved. Also attractive is the prospect of measuring temperature, density, and velocity 

with one instrument. RDV possesses significant advantages over Rayleigh Doppler 

velocimetry although the technique ia much more complicated. 

RDV has been demonstrated in a variety of forms. For example, CSRS (Ref. 15), 
CARS (Ref. 16), IRS (Ref. 17), stimulated Raman gain spectroscopy (SRGS) (Ref. 18) have 
all been demonstrated in wind tunnel environments. High resolution CARS and CSRS have 
been proposed as a method for measuring turbulence (Ref. 19). 

Exemplar applications of Raman techniques have been demonstrated in the Unitary 
wind tunnel at NASA Langley Research Center, illustrating an Important utility in high 
speed flows (Ref. 20, 21). Specifically, IRS was used to measure temperature, density, 
and velocity. Figure 5.9 shows the somewhat unique experimental setup, which 
incorporates some finesse necessary for successful operation in a wind tunnel 
environment. The following are key features of the setup: 

1. Molecules in the sample volume are illuminated by focused beams at two different 

wavelengths. One, the probe beam, is produced by a narrow-band (< 10 MHz 

linewidth) CW argon laser operating at a power of 0.5 watts at 514.5 nm 
wavelength. The other (the pump beam) is produced by an argon-laser-pumped dye 
laser whose output is amplified by a pulsed dye amplifier producing a 5 nsec, 7 mj 
pulse at 584.6 nm wavelength. The pulsed dye amplifier is pumped by a Q-switched, 
frequency-doubled YAG laser. 

2. After the beams have passed through the sample volume and to the other side of the 
tunnel, the pump beam is removed by a dichroic filter and dumped. 

3. The probe beam is returned by a retroreflector to focus again (from the opposite 
side) on the molecules in the sample volume. 

4. The detector senses light changes in the probe beam caused by two different types 
of interaction. One of the interactions is with the two beams propagating in the 
same direction (when the photons enter the sample volume for the first pass). The 
second Interaction is with the two beams counterpropagating (when the probe beam 
returns to the sample volume from the retroreflector). 

5. The intensity change in the probe beam is measured while tuning it over the Raman 
spectrum. This intensity decreases most when the energy of the probe beam photons 
are properly matched with the energies of the pump beam photons to allow a 
resonance with the final energy state of the molecule. The probe beam drives the 
molecular energy up to a virtual energy state and the pump beam stimulates it down 
to an allowed state. The emitted photon has energy a which is the difference in 
the probe and pump photon energies. This is called the IRS resonance. 
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6. The Doppler shift effects caused by the actions of the nolecules is significantly 
different for the two types of interactions (copropagation and 
counterpropagation), and these produce two distinct peaks in the measured 
spectrun. The measurement of the separation of the two peaks provides an accurate 
measure of velocity. The measurement of the difference in broadening of the two 
peaks provides a neasure of tenperature and pressure. 


MTKQftCTE* 



Figure 5.9. Schematic of the optical configuration for vibration-free Raman Doppler 
velocimetry in a wind tunnel. (After Exton, Ref. 21) 


For molecules moving with a velocity V, an illuminating wavefront Is Doppler 
shifted by (k - k ) • V (see Chapter 4), where k Q is the illumiating wavefront 
and k is the observed, scattered wavefront and k - »/2*. Both the pump beam, k , and 
the probe beam, (^ Dro be^' are " 8eenM by the molecule as Doppler shifted wavefronts. 
Therefore, the waveiengch at which a resonance is observed when tuning the wavelength of 
the probe beam depends on the velocity of the molecule. Without Doppler shift the 
resonance occurs when u> Drobe ■ u ♦ A where a is the associated Raman frequency. The 
Doppler shifted resonance occurs p wnen 


“pump “ (“probe + (^ptunp - k probe ) • V (5.13) 

The difference between the resonant wavelengths for the copropagating beams interaction 
and counterpropagatlng beams interaction can be determined by applying Equation 5.13 to 
be 
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where 9 is the angle between the probe beam direction and the velocity vector. In this 
manner, velocity can be measured without resorting to an absolute frequency 
measurement. Likewise, pressure and temperature can be determined from frequency 
difference measurements. 


Figure 5.10 is an actual plot, produced by Exton, et al for a Mach 2.5 free-stream 
flow, of the power loss from the probe laser during the 5 nsec pump beam periods as the 
wavelength is tuned for the copropagating and the counterpropagatlng interactions. The 
scan time required was about two minutes with a time constant of about 0.3 seconds. The 
velocity of the flow is computed from the peak separation. For this case, the frequency 
difference between the peaks was 1.00 GHz corresponding to a velocity of 602 m/sec 
(Mach 2.37). 



Figure 5.10. Measured spectral data (J - 10) for a Mach 2.50 free-stream flow in the 
Unitary wind tunnel (Re - 6.56 x 10*). The scan time was 2 min. with a 
time constant of -0.3 sec. (After Exton, Ref. 21) 
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The pressure and temperature are computed by measuring the two linebreadths. The 
two linebreadths were ,33 GHz for the copropagating interaction and 1.98 GHz for the 
counterpropagating interaction. A family of theoretical curves of the Voight 
linebreadths (temperature and pressure broadening) can then be used to derive temperature 
and pressure from the two measured linebreadths. The corresponding pressure and 
temperature for this case were 35 torr and 161 K respectively. 

The authors show that the precision in this measurement was about 17% for 
temperature and pressure and 5% for velocity. With an ideal signal the precision of this 
method would be about 2% for temperature and pressure and 1/2% for velocity. The 
technique is shown analytically to be applicable in the hypersonic regime at pressure 
down to 10 torr. Below that pressure, the signal begins to deteriorate below detectivity 
without other processing. With other signal enhancement procedures, this technique can 
be expected to cover virtually the entire current interest range in high speed 
aerodynamics and will doubtless prove to be an extremely important diagnostic tool for 
modern aerodynamics. 

Coherent Raman techniques have seen extensive application in combustion 
diagnostics. Eckbreth and associates (see additional reading) have refined CARS 
techniques to routine use in combustors and jet engine development, and other combustion 
laboratories have used CARS extensively. The complexity of the hardware as illustrated 
in Figure 9 discourages its use in large-scale wind tunnel tests although a wider 
applicability can be anticipated. 

Raman methods provide the unique capability to monitor several species 
simultaneously at a point. This can be contrasted to absorption and fluorescence which 
monitor a single species for each wavelength, therefore, requiring wavelength tuning to 
cover several species. 
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FIGURE Al—REPRESENT ATION OF PLANE WAVE FRONTS 



FIGURE A2. REPRESENTATION OF PLANE WAVEFRONTS JF WAVELENGTH 
SLIGHTLY LARGER THAN THAT OF FIGURE Al. 
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FIGURE A4 


SPHERICAL WAVE9 OF 9LIGHTLY LARGER 
WAVELENGTH THAN IN FIGURE A3. 
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